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FOREWORD 


The  research  work  reported  herein  was  conducted  by  the  McDonnell  Aircraft 
Company (MCAIR) ,  McDonnell  Douglas  Corpora'' ion,  St.  Louis,  Mo.  for  the  Metals 
and  Ceramics  Division,  Air  Force  Materials  Laboratory,  Air  Fa  a-  Systems 
Command,  Wright-Patt^rson  Air  Force  Base,  Ohio,  under  contract  F33615-72-C-1882. 
This  contract  was  initiated  under  Project  7351,  "Metallic  Materials",  Task 
735106,  "Behavior  of  Metals",  with  Mr.  R.  C.  Donat,  AFML/LLN,  acting  as  project 
engineer. 

The  study  was  conducted  by  the  Structural  Research  Department  of  MCAIR 
under  the  guidance  of  Messrs,  E.  D.  Bouchard  and  M.  L.  Rand  and  under  the 
direct  supervision  of  Mr.  L.  F.  Impellizzeri,  as  Principal  Investigator,  with 
assistance  from  Messrs.  A.  E.  Siegel  and  R.  A.  McGinr4*?.  Other  Structural 
Research  Department  personnel  making  valuable  contributions  to  the  program  were 
Messrs.  T.  L.  Benton  and  F.  R.  Foster.  The  period  covered  by  this  effort  was 
1  May  1972  to  1  May  1973.  The  report  was  submitted  by  the  authors  on.  30  May 
1973. 

This  report  has  been  reviewed  and  is  approved. 
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ABSTRACT 


Structural  reliability  analysis  procedures  were  evaluated  for  estimating 
the  variability  in  fighter  aircraft  fatigue  performance.  The  expected  magni¬ 
tude  of  this  variability  was  determined  based  on  an  investigation  of  scatter 
in  fatigue  test  results  for  aluminum  structures.  This  study  revealed  that 
(a)  the  scatter  for  spectrum  fatigue  tests  is  considerably  less  than  for  con¬ 
stant  amplitude  fatigue  tests,  (b)  the  Weibull  probability  distribution 
provides  a  better  fit  of  the  spectrum  fatigue  data  than  the  log-normal,  and 
(c)  the  Weibull  shape  parameter  a  is  5.27.  The  a  value  was  determined  from 
1060  spectrum  test  results  of  which  243  were  full-scale  airplane  and  airplane 
component  tests.  These  included  the  F-3H  Demon  wing  and  horizontal  tail,  the 
F-4  Phantom  II  wing  box  beam,  the  Lockheed  wing  test  panel,  the  F-9F  Panther 
wing,  the  Navy  Lab  box  beam,  the  P-51  Mustang  wing,  the  C-46  transport  wing, 
and  the  British  Piston  Provost  wing.  In  addition  to  these  studies  of  experi¬ 
mental  data,  theoretical  analyses  were  performed  yielding  the  mathematical 
probability  distribution  for  a  Weibull  based  scatter  factor  which  is 

S  “(R/l-R)1^01 

where  S  is  the  scatter  factor  and  R  is  reliability  or  the  probability  of  no 
failure.  VGH  and  load  factor  counting  accelerometer  data  from  the  F-4 
fighter  airplane  were  utilized  to  correlate  that  airplane’s  laboratory  and 
service  fatigue  experience.  Probable  minimum  service  liv^  considering  the 
F-4  fleet  size  and  individual  airplane  usage  were  computed  based  on  the 
Weibull  based  scatter  factor  and  order  statistics.  The  combined  effect  of 
fatigue  test  scatter  and  usage  severity  scatter  was  derived  utilizing  a 
joint  scatter  factor  concept.  Three  fatigue  critical  locations  on  the  F-4 
airplane  were  considered  to  be  amenable  to  analysis  using  the  methods  of 
this  report.  The  correlation  was  excellent  for  one  of  these,  but  not  for  the 
other  two.  Fabrication  variations  in  a  redundant:  load  path  joint  and  outer 
wing  buffeting  were  considered  the  probable  cause  for  the  less  than  favorable 
correlation. 
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1.  ABBREVIATIONS 

axp  =>  exponential  function 
log  =  common  logarithm 

BLUE-3  best  linear  unbiased  estimator  or  estimate 

2.  SYMBOLS 

a  =  Weibull  distribution  shape  or  scatter-controlling  parameter 

6  =  Weibull  distribution  scale  parameter  or  characteristic  life 

f(  ].  =  probability  densit1'  function  of  the  parameter  within  the 

parenthesis 

F(  )  «*  cumulative  probability  distribution  function  of  the  parameter 

within  the  parenthesis 

y  =  log-normal  distribution  scale  parameter  or  mean  life 

n  =  number  of  failure  observations  in  the  test  sample 

R  =  reliability  of  a  randomly  chosen  fleet  member  or  the  probability 
of  no  failure  (probability  that  the  ratio  of  life  in  the  labora¬ 
tory  to  life  in  service  is  less  than  the  scatter  factor) 

S  =  scatter  factor 

a  =  log-normal  standard  deviation 


x  ■  life  of  individual  laboratory  specimen 

y  ■  life  of  a  randomly  selected  specimen  or  aircraft 

t  «  usage  severity  scatter  ■  actual  counts/estimate^  counts 
”(  )  ■  Gamma  function  of  argument  within  the  parenthesis 

3.  ACCENT  MARKS  DENOTING  ESTIMATION 

(a)  Bar,  as  in  a  or  a,  denotes  the  sample  point  estimate  for  one  group 
of  specimens 

(b)  Double  bar,  as  in  a  or  a  denotes  the  sample  estimate  for  more  than 
one  group  of  specimens. 
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1.  INTRODUCTION 


1.1  General 

The  Metals  and  Ceramics  Division  of  the  Air  Force  Materials  Laboratory, 

Air  Force  Systems  Command,  has  sponsored  research  efforts  on  application  of 
reliability  methods  to  the  estimation  of  probable  aircraft  structural  fatigue 
performance.  Results  of  these  efforts  are  reported  in  References  (1)  and  (2). 
The  reliability  analysis  approach  has  been  evaluated  for  cargo  (C-130)  and 
tanker  (KC-135)  type  aircraft.  The  purpose  of  the  present  program  has  been 
to  extend  reliability  analysis  methods  to  fighter  aircraft.  This  objective 
was  achieved  through  utilization  of  fatigue  test  experience  and  service 
repeated  load  data  which  has  been  generated  in  conjunction  with  the  heDonnell 
Phantom  II  (F-4)  airplane. 

1.2  Traditional  Approach  to  Structural  Reliability 

Structural  fatigue  integrity  is  currently  designed  into  present  day  air¬ 
craft  utilizing  a  rather  well  organized  plan.  Spectrum  fatigue  tests  are 
performed  in  the  laboratory  on  the  full  scale  structure  with  the  target  goal 
being  some  multiple  of  the  expected  service  life.  This  is  the  traditional 
scatter  factor  approach  which  has  been  used  for  quite  a  few  years.  There 
have  been  aircraft  accidents  during  this  period,  but  there  have  been  many 
causes  othe  than  airframe  structural  deficiencies;  e.g.,  pilot  error,  power 
plant  failure,  electrical  system  malfunction.  The  percentage  of  accidents 
resulting  from  airframe  fatigue  failure  has  been  small.  However,  fatigue 
failures  have  led  to  rather  large  structural  maintenance  expenditures  over 
the  years.  In  general,  these  have  been  caused  by  unexpectedly  severe  usage, 
an  unusual  design  detail,  or  a  material  deficiency;  but  the  simple  scatter 
factor  approach  may  have  been  partly  responsible. 

1.3  Additional  Considerations 

There  are  a  number  of  questions  that  the  traditional  approach  to  struc¬ 
tural  reliability  does  not  take  into  account.  For  example:  Should  the  same 
scatter  factor  be  used  for  fighters,  bombers,  transports,  commercial  airliners, 
for  different  materials,  for  different  structural  details  and  arrangements, 
etc.?  Should  not  the  number  of  laboratory  test  articles  affect  the  required 
scatter  factor?  Should  the  expected  fleet  size  affect  the  required  scatter 
factor?  General  aircraft  attrition  due  to  factors  other  than  fatigue  is 


another  Important  consideration.  The  fleet  size  is  reduced  by  attrition.  The 
consideration  of  damage  tolerance  could  very  well  lead  to  different  allowable 
scatter  factors  for  different  components  on  the  same  airplane.  The  underlying 
problem  is  the  need  for  development  of  a  method  to  compute  the  required  scatter 
factor  for  any  preselected  reliability  level  and  the  associated  time  to  first 
failure . 

In  terms  of  structural  reliability,  there  are  so  many  unknowns  that  the 
traditional  approach  certainly  offers  some  advantages.  For  example,  it  is 
not  possible  to  accurately  predict  the  usage  that  an  airplane  will  experience 
in  some  future  time  period.  However,  if  certain  variables  of  reliability 
can  be  quantified  without  adding  undue  complexity  to  an  airplane  design,  it 
should  be  done.  A  list  of  pertinent  considerations  would  include: 

(a)  basic  fatigue  scatter  for  material  in  question  -  defined  by  a 
distribution  function  and  certain  scaling  parameters,  e.g. ,  mean 
and  standard  deviation. 

(b)  number  of  full  scale  laboratory  test  articles  -  these  would  define 
the  probability  distribution  sample  mean. 

(c)  expected  size  of  total  fleet. 

(d)  expected  service  usage  severity  scatter  from  airplane  to  airplane 
within  the  fleet. 

(e)  attrition  considerations  -  loss  of  aircraft  due  to  reasons  not 
associated  with  structural  failures,  or  retirement  of  aircraft 
due  to  obsolescence. 

(f)  damage  tolerance  of  primary  structural  components. 

(g)  degree  of  simulation  of  the  service  loading  including  environmental 
effects  by  laboratory  testing. 

(h)  level  of  confidence  in  the  predicted  repeated  load  usage  severity. 

(i)  type  of  airplane,  e.g.,  fighter,  bomber,  commercial  transport,  etc. 

The  first  five  considerations  are  strictly  concerned  with  the  statistical 
aspects  of  the  overall  structural  reliability  problem.  In  other  words  (a) 
through  (e)  provide  an  estimate  of  the  failure  probability  associated  with  any 
given  scatter  factor  and  any  number  of  airplanes  out  of  the  total  fleet. 

Considerations  (f)  through  (k)  have  to  do  with  design,  maintenance,  dollar, 
and  safety  trade-offs.  What  magnitude  of  failure  probability  is  acceptable? 

What  are  the  potential  consequence  of  a  failure  in  particular  structural 
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components?  How  much  structural  weight  can  be  tolerated  to  minimize  the 
failure  probability?  Some  compromise  must  be  reached  between  a  high  per¬ 
formance  design  on  the  one  hand  and  one  hundred  percent  safety  of  flight 
and  zero  maintenance  expenditures  on  the  other.  In  reaching  this  compromise, 
it  is  important  to  consider  damage  tolerance  and  the  effectiveness  of  perio¬ 
dic  inspection  in  detecting  fatigue  cracks  before  they  can  reach  catastrophic 
proportions.  The  compromise  is  also  affected  by  the  type  of  aircraft.  Some 
difference  in  philosophy  would  normally  be  expected  between  military  and 
commercial  aircraft.  Further  differences  might  be  expected  between  fighters, 
bombers,  and  cargo  type  aircraft.  On  each  new  airplane,  decisions  regarding 
these  areas  must  be  made  within  the  context  of  available  funding,  schedule, 
contingency  planning,  and  national  priorities. 

1.4  Tailored  Scatter  Factor 

Reliability  requirements  could  be  more  effectively  established  if  some 
quantitative  measure  of  failure  probability  was  calculable  for  various  design 
trade-offs.  If  the  design  scatter  factor  is  increased  or  decreased  for  a 
given  component,  how  will  this  affect  failure  rates  in  a  fleet  of  aircraft? 

How  will  maintenance  expenditures  be  affected?  A  change  in  the  design  scatter 
factor  will  obviously  affect  structural  weight.  The  influence  of  structural 
weight  on  cost  and  performance  must  then  be  considered.  The  basic  payoff  to 
be  derived  from  using  mathematical  probability  methods  In  structural  fatigue 
design  is  to  provide  for  the  calculation  of  a  "tailored  scatter  factor"  that 
would  fit  the  particular  set  of  circumstances  existing  during  an  airplane's 
design  stage.  Aircraft  program  managers  would  be  able  to  select  a  magnitude 
of  scatter  factor  associated  with  whatever  reliability  level  is  deemed 
appropriate. 

1 . 5  Utilization  of  Phantom  II  Data 

The  HCAIR  Phantom  II  (F-4)  aircraft  was  considered  ideally  suited  for  an 
investigation  of  structural  reliability  design  methods  for  fighter  aircraft, 
More  than  4,000  F-4's  have  been  delivered  and  are  being  used  by  the  U.S.  Air 
Force,  Navy,  and  Marines.  The  airplane  is  currently  being  manufactured  in 
various  models  for  the  U.  S.  military  services  and  also  foreign  governments. 

A  comprehensive  and  detailed  flight  loads  monitoring  program  has  been  it?  effect 
continuously  on  the  F-4  since  its  initial  delivery  to  the  Navy.  This  program 
has  provided  data  from  counting  accelerometers  and  VGH  recorders  for  training 
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as  well  as  combat  operations  in  Southeast  Asia.  The  VGH  instruments  provide 
an  analog  recording  of  load  factor,  speed,  and  altitude.  Gross  weights  corres¬ 
ponding  to  the  VGH  data  are  provided  from  individual  flight  records.  The 
counting  accelerometers  record  load  factor  occurrences  which  equal  or  exceed 
four  predetermined  acceleration  levels.  More  than  2,000,000  flight  hours  of 
load  factor  exceedance  data  have  been  accumulated  in  the  counting  accelerometer 
program.  The  reduction  of  these  data  has  been  completely  automated.  A  typi¬ 
cal  computer  print-out  showing  the  load  factor  usage  history  on  an  individual 
airplane  is  shown  on  Figure  1.  The  laboratory  fatigue  test  program  on  the 
F-U  has  also  been  extensive;  it  has  included  five  complete  wing-center  fuse¬ 
lage  test  articles,  two  half-wings,  numerous  smaller  component  tests,  and 
hundreds  of  element  tests.  A  brief  summary  of  the  F— h  full-scale  testing  is 
shown  in  Figure  2. 

1.6  Overall  Structural  Reliability  Program  Outline 

The  subject  program  for  evaluating  structural  reliability  analysis  methods 
for  fighters  was  based  largely  on  F-h  experience  and  was  conducted  in  four  phases 
Phase  I:  Survey  of  Fatigue  Test  Scatter 

Phase  II:  Evaluation  of  Usage  Severity  Scatter 

Phase  III:  Documentation  of  Laboratory  Test  Results 

Phase  IV:  Correlation  of  Laboratory  and  Service  Failure  Experience 

The  general  approach  included  a  literature  Burvey  to  determine  basic 
fatigue  scatter  trends  for  fighter  type  repeated  load  spectra.  ThiB  was  per¬ 
formed  in  Phase  I.  A  significant  part  of  this  effort  was  the  study  of  in-house 
test  results.  MCAIR  has  performed  about  1000  aluminum  fatigue  testa  on  simple 
elements  and  components.  About  two-thirds  of  these  were  conducted  using 
fighter  type  repeated  load  spectra. 

Study  of  scatter  in  severity  of  repeated  load  history  from  airplane  to 
airplane  in  a  given  fleet  was  an  important  part  of  the  program.  This  repre¬ 
sented  Phase  II  of  the  total  effort.  In  recent  years,  numerous  military 
aircraft  have  been  equipped  with  load  factor  exceedance  counters.  Millions  of 
flight  hours  have  been  logged  with  these  instruments.  The  analysis  of  these 
data  provides  an  excellent  ins i gilt  into  loads  scatter  trends.  The  F-U  air¬ 
plane  loads  monitoring  program  is  the  most  extensive  of  its  kind  ever  under¬ 
taken.  The  F-U  has  logged  more  than  5 >350*000  flight  hours  in  numerous  theaters 
of  operation.  As  mentioned  earlier,  load  factor  exceedance  data  were  available 
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Figure  1 

Typical  Computer  Print  Out  Showing  Load  Factor  History  of  an  Individual  Airplano 
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Figure  2 

Summary  of  the  F-4  Full  Scale  Laboratory  Fatigue  Tests 
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from  counting  accelerometers  for  more  than  2,000,000  flight  hours.  These 
were  supplemented  with  about  9,000  hours  of  VGH  data  which  were  also  studied 
in  the  second  phase  of  the  subject  program.  The  extent  of  the  available  data 
is  summarized  in  Figure  3.  Note  that  about  650,000  flight  hours  of  counting 
accelerometer  data  were  from  combat  operations. 

Comparison  of  laboratory  and  service  failures  was  the  obvious  next  step. 
Tnis  was  accomplished  in  Phases  III  and  IV.  The  study  included  a  detai  led 
investigation  of  test  results  on  the  five  full-scale  F-h  test  articles  and 
the  two  half-wings  in  a  number  of  key  areas.  The  time  to  failure  for  each  key 
area  on  each  teBt  article  was  determined.  A  careful  examination  was  made  of 
any  possible  differences  in  laboratory  loading  and  in-flight  loading.  The 
study  of  VGH  data  was  used  to  assist  in  this  examination.  The  final  stage 
included  a  comparison  of  the  service  life  predictions  and  actual  service 
failures.  The  predictions  were  based  on  the  laboratory  testing,  counter  and 
VGH  data,  and  o^atter  considerations  evaluated  in  Phases  I  and  II. 


VGH  Data 

1965  -  3000  hr  Training  Data 

1966  -  600  hr  Training  Data 
1P39  -  1800  hr  SEA  Combat  Data 
1970  -  3600  hr  SEA  Combat  Data 

Counting  Accelerometer  Data 

F-4  Fighter  Aircraft:  RF*4  Reconnaissance  Aircraft: 

2,394,000  Total  hr  402,000  Total  hr 

517,000  Combat  hr  |  124,000  Combat  hr 

1,877,000  Training  hr  278,000  Training  hr 

Total  Flight  Time,  all  F-4  Aircraft  -  5,350,000  hr 

Figure  3 

Summary  of  Available  F-4  VGH  and  Counting  Accelerometer  Data 
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2.  PHASE  1  -  SURVEY  OF  FATIGUE  TEST  SCATTER 

2.1  General 

The  nature  of  fatigue  is  such  that  there  is  very  little  that  "an  be 
determined  in  a  purely  analytic  sense.  Whatever  is  done  on  a  theoretical 
basis  usually  allows  simple!  interpolation  or  modest  extrapolation  of  experi¬ 
mental  information.  Fatigue  scatter  falls  in  this  same  category.  The  theory 
of  probability  la  applicable  but  there  are  many  mathematical  functions  that 
might  be  used  and  each  one  has  certain  constants  that  need  to  be  defined. 

An  extensive  detailed  study  of  fatigue  teat  results  is  required  to  determine 
which  mathematical  function  is  most  appropriate  and  what  constants  are  neces¬ 
sary  to  provide  the  best  fit. 

There  are  two  mathematical  functions  in  particular  that  have  been  widely 
used  to  describe  fatigue  scatter;  the  log-normal  and  the  Weibull,  and  pre¬ 
dominately  the  two  parameter  version  of  each.  In  recent  years,  the  Weibull 
has  gained  favor  partly  due  to  a  physical  argument  concerning  its  risk  func¬ 
tion  compared  to  that  of  the  log-normal.  The  risk  function  for  fatigue  is 
defined  as  the  probability  of  failure  on  the  (N  +  X)th  cycle  given  that  the 
structure  has  sustained  N  cycles.  The  Weibull  has  a  continuously  increasing 
risk  function  while  the  log-normal  risk  function  first  increases  and  then  de¬ 
creases  with  increasing  cycles.  Since  fatigue  is  considered  to  be  a  veanmt 
process,  an  ever  Increasing  risk  function  seems  more  reasonable.  A  plausible 
physical  argument  can  also  be  made  in  favor  of  the  log-normal.  It  is 
generally  agreed  that  fatigue  originates  at  some  kind  of  discontinuity  in  the 
metal  mic restructure.  It  is  highly  improbable  that  auy  given  specimen  would 
be  totally  homogeneous  and  without  discontinuities  in  its  microstructure. 

The  log-normal  risk  function  does  not  start  to  decrease  until  the  probability 
level  reaches  0.999999  which  is  associated  with  a  very  long  life.  It  is  highly 
improbable  for  a  specimen  to  last  that  long,  but  if  it  does,  it  could  be  an 
indication  that  the  specimen  is  without  discontinuities  and  therefore  the  risk 
function  should  start  to  decrease. 

The  discussion  in  the  preceding  paragraph  is  not  presented  as  an  argu¬ 
ment  to  prove  that  the  log-normal  distribution  is  correct.  It  is  included 
simply  to  indicate  that  experimental  evidence  and  not  theoretical  analysis 
must  be  the  primary  determining  factor. 

Once  a  distribution  function  is  established,  the  next  step  is  to 
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determine  the  constants  of  the  distribution.  Generally,  two  constants  are 
sufficient.  One  of  these  defines  the  magnitude  of  scatter  and  the  other 
some  kind  of  average  life.  Average  life  is  obviously  peculiar  to  each  type 
of  airplane;  but  fatigue  scatter  is  generally  considered  similar  on  different 
airplanes,  at  least  if  they  are  made  of  the  same  material  and  with  the  same 
type  of  construction  and  loaded  in  generally  the  same  manner.  This  leads  to 
the  conclusion  that  the  constant  defining  the  magnitude  of  scatter  can  be 
determined  by  reviewing  past  test  results.  This  was  a  basic  underlying 
assumption  in  the  literature  survey  and  subsequent  study  of  fatigue  scatter 
trends  reported  in  this  Phase  I  portion  of  the  program. 

2.2  Literature  Survey  of  Fatigue  Test  Data 

As  a  result  of  an  extensive  literature  survey,  more  than  2,400  groups 
of  fatigue  test  data  made  up  of  over  12,000  specimens  were  collected.  This 
search  was  limited  to  the  2000  and  7000  series  aluminum  alloys,  which  in¬ 
clude  2020,  2024,  2124,  7075,  7079,  7175,  and  7178.  It  included  test  data 
from  MCAIR  as  well  as  from  numerous  other  aerospace  companies  and  government 
agencies.  In  particular,  the  MCAIR  fatigue  test  results  are  from  F-15,  F-4, 
JXI-10,  and  earlier  in-house  teat  programs.  Specimens  from  sheet,  plate,  bar 
extrusion,  and  forging  are  included  in  configurations  ranging  from  simple 
notched  elements  with  and  without  holes  to  small  components,  as  well  as  full- 
scale  structures.  Test  loading  conditions  include  axial  loading,  both  ten¬ 
sile  and  compressive,  and  flexural  loading.  The  greater  portion  of  the  test 
loading  data,  however,  was  axial.  The  fatigue  tests  consisted  of  both 
constant  amplitude  and  spectrum  loading.  References  of  both  the  literature 
search  and  the  MCAIR  data  are  listed  in  Appendix  I.  Only  those  references 
in  which  the  data  were  analysed  and  summarised  are  listed.  The  following 
variables  of  the  collected  fatigue  data  were  indexed  so  that  comparisons 
could  be  made  among  different  groups  of  specimens;  type  of  structure,  type 
of  specimen,  material,  type  of  test  machine,  and  type  of  loading.  The  actual 
test  results  analysed  in  detail  in  this  survey  are  presented  along  with  the 
appropriate  index  numbers  in  Appendix  It.  In  addition,  the  appropriate 
statistical  parameters  for  the  log-normal  and  Veibull  distributions  for  each 
group  of  specimens  are  also  presented  in  Appendix  IX. 
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2.2.1  Scatter  Trends  for  Constant  Amplitude  and  Spectrum  Loading  - 

Reference  (1)  presented  an  extensive  literature  survey  and  analysis  of  fatigue 

test  data.  A  summary  of  these  data  is  listed  below: 

Number  of  Specimens  6000* 

Number  of  Groups  1250 

Average  Standard  Deviation  0.168 

The  numbers  above  are  for  aluminum  alloy  specimens  of  varying  configuration 
subjected  to  either  spectrum  or  constant  amplitude  loading.  These  teat  re¬ 
sults  from  Reference  (1)  were  subdivided  into  groups  according  to  loading  as 


follows: 

Constant  Amplitude  Spectrum 
_ Loading _  Loading 

Number  of  Specimens  5000*  1000* 

Number  of  Groups  1096  154 

Average  Standard  Deviation  0.180  0.083 


The  difference  between  the  constant  amplitude  and  spectrum  loading  standard 
deviations  is  quite  significant.  Xt  is  possible  of  course  that  part  of  the 
difference  could  be  a  result  of  other  variations  besides  loading  in  the  two 
groups  such  as  a  comparison  of  perhaps  more  open  hole  specimens  in  one  group 
to  more  complex  lap  Joint  specimens  in  the  other.  So  in  order  to  more  care¬ 
fully  evaluate  this  difference  between  spectrum  and  constant  amplitude  data, 
70? 5-T6  aluminum  open  hole,  axial  loaded  specimens  from  Reference  (I)  were 
studied  with  the  following  results: 


Constant  Amplitude  Spectrum 

Loading  Loading 


Number  of  Specimens  516  146 
Number  of  Groups  138  26 
Average  Standard  Deviation  0.211  0.125 


The  same  trend  as  in  the  larger  data  set  is  evident  from  the  above  comparison. 
A  possible  explanation  for  the  lower  scatter  of  the  spectrum  data  is  that  only 
the  higher  load  levels  in  the  spectrum  produce  the  major  part  of  the  total 
damage;  that  is  to  say,  the  spectrum  test  results  are  effectively  low  cycle 
fatigue  even  though  the  total  number  of  cycles  in  the  spec t rum  might  be  quite 
large.  Low  cycle  fatigua  generally  ahova  lower  scatter  (smaller  standard 
deviation)  than  high  cycle  fatigue.  This  trend  is  illustrated  by  the  histo¬ 
gram  of  Figure  4.  The  data  in  this  graph,  which  include  970  specimens  in  269 

♦Estimated 
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Footed  Oats  •  7-075  Material  ■  Constant  Amplitude  Loading 


Figure  4 

Qbsarwjd  Variation  of  Fatigue  Scatter  with  Cyclic  Life  for  970  Aluminum 
7075  Constant  Amplitude  Loaded  Specimens  -  Log  -  Normal  Distribution 


groups,  are  for  7075-T6  aluminum  axially  loaded,  edge  notched  or  unloaded 
hole  specimens. 

To  validate  and  further  understand  the  trends  found  from  Reference  (1) , 
MCAIR  fatigue  test  data  were  assembled  as  a  group  for  scatter  trena  studies. 
The  HCAIR  data  were  obtained  from  44  teat  reports  consisting  of  288  groups 
of  947  specimens.  The  HCAIR  fatigue  teat  results  included  7000  and  2000 
series  aluminum  alloys,  various  specimen  configurations,  with  constant  ampli¬ 
tude  or  spectrum  loadings.  Separating  the  data  based  on  the  type  of  loading 
similar  to  the  Reference  (1)  grouping  gives  the  following  comparison: 

Constant  Amplitude  Spectrum 


Loading _  Loading 


Number  of  Specimens  287  656 
Number  of  Groups  93  196 
Average  Life  in  Cycles  140,000  200,000 
Average  Standard  Deviation  .199  .109 


The  spectrum  data  again  have  the  lower  scacter.  It  is  further  hypothesised 
that,  in  addition  to  the  low  cycle  fatigue  consideration  discussed  in  the 
preceding  paragraph,  the  type  of  fatigue  test  machine  also  could  influence 
the  magnitude  of  scatter.  Spectrum  teats  are  in  a  great  majority  of  eases 
run  with  oithor  a  servo  or  solenoid  load  controlling  device,  while  the 
mechanical  shaker  type  of  test  machine  is  more  often  used  for  constant 
amplitude  testing.  The  mechanical  shaker  type  machine  has  less  load  control 
accuracy.  In  order  to  verify  the  hypothesis  of  differences  in  scatter  due 
to  the  type  of  fatigue  test  machine,  the  HCAIR  Constant  amplitude  test  data 
were  analysed.  The  results  are  as  follows: 


Mechanical  Servo 
Shaker  Control 


Number  of  Specimens  73  128 
Number  d£  Groups  27  39 
Average  life  in  Cycles  178,000  206,000 
Average  Standard  Deviation  0.238  0.158 


Although  the  amount  of  fatigue  test  data  is  small,  the  hypothesised  trend 
of  inherently  mors  scatter  with  the  mechanical  shaker  machine  than  the  servo 
control  machine  is  evident. 

2.2.2  No-Load  Transfer  Elegant  Fatigue  Teat  Data  -  Baaed  on  the  re¬ 
sults  discussed  in  Section  2.2.1,  spectrum  fatigue  test  scatter  is  leas  than 
constant  amplitude  teat  scatter  for  two  distinct  reasons: 

(a)  low  cycle  fatigue  effect,  and 

(b)  fatigua  machine  effect. 
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Since  fighter  aircraft  are  subject  to  c  spectrum  of  loads  in  service,  the 
primary  emphasis  in  this  study  of  fatigue  test  scatter  was  directed  toward 
spectrum  data. 

A  literature  survey  was  undertaken  to  obtain  additional  spectrum  fatigue 
test  results.  These  additional  data  were  combined  with  MCAIR's  spectrum  test 
data  and  the  spectrum  test  data  from  Reference  (1).  Only  axially  loaded, 
edge  notched  or  unloaded  hole  specimens  were  included.  Fatigue  test  data 
obtained  from  mechanical  shaker  machines  were  not  included  for  the  reasons 
discussed  in  Section  2.2.1.  A  compilation  of  the  results  of  this  study  are 
given  below: 

No-Load  Transfer  Spectrum  Fatigue  Test  Data 


Aluminum 

Number  of  Specimens 

1176 

Number  of  Groups 

239 

Effective  Average  Life  in  Cycles 

17  s.400 

Average  Standard  Deviation 

.0855 

Average  wieibull  Shape  Parameter 

6.36 

The  effective  average  life  denoted  above  is  defined  as  the  number  of  higher 
load  level  cycles  to  failure.  The  higher  load  levels  are  defined  as  the  top 
two-thirds  of  the  load  levels  of  the  spectrum.  If  there  are  nine  load  levels 
in  a  particular  spectrum,  only  the  cycleB  due  to  the  six  higher  load  levels 
would  be  considered  effective.  This  technique  was  used  to  minimize  the  mis¬ 
leading  effect  of  the  large  number  of  lower  level  load  cycles  in  a  spectrum 
that  produces  only  a  small  percentage  of  the  total  fatigue  damage.  It 
should  be  noted  that  he  average  standard  deviation  is  not  simply  a  numerical 
average  but  the  unbiased  estimate  of  the  true  population  standard  deviation. 
Similarly,  the  WeiLull  number  is  the  weighted  average  of  the  best  linear 
unbiased  estimator  (BLUE)  of  the  Weibull  distribution  shape  parameter  given 
in  Reference  (3).  The  Weibull  shape  parameter  is  a  measure  of  fatigue  test 
scatter  similar  to  the  standard  deviation  except  in  reverse;  a  large  value 
Indicates  less  scatter  while  a  smaller  value  indicates  more  scatter. 

Additional  studies  were  conducted  by  grouping  the  data  in  different 
ways.  The  7075  aluminum  spectrum  fatigue  test  results,  for  example,  were 
subdivided  according  to  stress  concentration  and  material  form.  Neither 
3tress  concentration  nor  material  form  were  shown  to  have  a  significant 
effect  on  fatigue  scatter.  Another  study  on  these  7075' data  included  a 
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histogram  of  standard  deviations  similar  to  Figure  4  except  it  was  for 
spectrum  data  plotted  versus  effective  average  life.  This  histogram  indi¬ 
cated  that  spectrum  fatigue  scatter  is  reasonably  independent  of  average 
life. 

2.2.3  Effects  of  Special  Specimen  Preparation  -  Specimen  finish  in  the 
fatigue  critical  area  was  one  of  the  parameters  investigated  for  effect  on 
fatigue  scatter.  A  relatively  large  group  of  the  no-load  transfer  element 
specimen  data,  generated  at  the  NASA  Research  Center  in  Langley,  West 
Virginia,  is  interesting  in  this  respect.  This  group  from  References  (47), 
(48),  (49),  (57),  (58),  and  (59)  of  Appendix  I  includes  a  total  of  682  edge 
notch  specimens;  the  edge  notch  itself  was  machined  with  extreme  care.  A 
comparison  of  these  data  to  all  others  is  a3  follows: 


Smooth  Finish 
Edge  Notch  Data 

All  Other  No-Load 
Transfer  Data 

Number  of  Specimens 

682 

494 

Number  of  Groups 

120 

119 

Average  Standard  Deviation 

.0747 

.0995 

Average  Weibull  Shape  Parameter 

7.19 

5.41 

The  effect  of  the  finish  appears  to  be  significant.  For  the  purposes  of  the 
NASA  work  (determination  of  the  effect  of  various  fatigue  spectrum  parameters) , 
the  smooth  finish  edge  notch  provided  a  very  useful  specimen  configuration. 
However,  an  extremely  smooth  finish  would  be  too  costly  for  typical  aircraft 
fabrication  procedures.  Therefore,  these  data  are  not  included  in  the  analy¬ 
sis  in  the  remaining  sections  of  this  report. 

2.2.4  Load  Transfer  Element  Fatigue  Test  Data  -  Lap  joint  and  double 
shear  joint  specimen  test  results  have  been  compiled  through  the  literature 
survey.  Pertinent  details  are  as  follows: 

Load  Transfer  Spectrum  Fatigue  Test  Data 


Aluminum 

Number  of  Specimens 

323 

Number  of  Groups 

59 

Average  Standard  Deviation 

.0998 

Average  Weibull  Shape  Parameter 

5.00 

It  is  of  interest  to  note  that  the  scatter  parameters  for  these  load  transfer 
element  test  data  are  similar  to  those  given  in  the  preceding  section  for  no- 
load  transfer  testinb  excluding  the  smooth  finish  edge  notch  data. 

2.2.5  Full  Scale  Structure  Fatigue  Test  Data  -  The  scatter  trends 
summarized  in  the  previous  three  sections  are  based  on  tests  of  element 


specimens  with  unloaded  holes  and  edge  notches  and  simple  joint  specimens. 

For  this  kind  of  element  specimen  and  also  for  full  scale  structure  fatigue 
test,  scatter  is  a  result  of  life  variations  at  a  stress  concentration.  It 
could  be  expected  therefore  that  the  scatter  trends  in  element  tests  and  full 
scale  structure  tests  would  be  similar.  However,  there  are  a  number  of  con¬ 
siderations  that  could  lead  to  differences.  There  is  generally  more  than  one 
single  location  on  a  full  scale  structure  where  the  failure  could  originate. 
For  example,  it  is  usually  the  case  that  there  are  a  number  of  holes  in  a  row 
that  are  equally  susceptible  to  fatigue.  This  kind  of  situation  where  there 
are  multiple  origin  possibilities  tends  to  reduce  scatter.  From  another 
viewpoint,  however,  full  scale  structure  might  be  expected  to  exhibit  more 
scatter  than  element  specimens.  An  airplane  is  a  complex  built-up  structure 
with  redundant  load  paths.  The  possibility  of  load  differences  at  the  criti¬ 
cal  location  from  one  test  article  to  another  is  therefore  greater  than  with 
a  simple  element  specimen.  These  internal  load  differences  would  tend  to 
increase  scatter. 

As  indicated  in  the  previous  paragraph,  it  is  not  a  certainty  that 
scatter  in  full  scale  testing  is  the  same  as  in  element  fatigue  testing.  A 
literature  survey  of  full  scale  airplane  and  component  fatigue  test  data  was 
therefore  undertaken.  A  total  of  243  spectrum  test  results  and  491  constant 
amplitude  test  results  were  found.  These  test  data  are  summarized  below: 

Full  Scale  Spectrum  Fatigue  Test  Data 


Spectrum 

Constant  Amplitudi 

Number  of  Specimens 

243 

4S1 

Number  of  Groups 

82 

143 

Average  Standard  Deviation 

.0985 

.1486 

Average  Weibuil  Shape  Parameter 

5.44 

3.70 

Note  that  these  full  scale  structure  data  show  the  same  trend  as  the  element 
test  results;  viz.,  the  average  standard  deviation  for  the  spectrum  test  re¬ 
sults  is  significant!  '  smaller  than  that  for  the  constant  amplitude  test 
results.  The  full  scale  structure  data  include  both  7075  and  2024  material. 

A  detailed  listing  of  the  full  scale  structure  spectrum  fatigue  test  results 
is  presented  in  Figure  5. 

2.2.6  Pooled  Spectrum  Fatigue  Test  Data  -  It  is  expected  that  the  scatter 
differences  between  7075  and  2024  aluminum  are  probably  small.  This  can  be 
determined  for  the  various  kinds  of  element  specimens  and  full  scale  structure 
discussed  in  the  three  preceding  sections.  The  data  comparisons  for  spectrum 
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Figure  5 

Full  Scale  Structure  Spectrum  Fatigue  Test  Data 


fatigue  test  results  are  as  follows: 

No-Load  Transfer  Element  Fatigue  Test  Data 


7075 

2024 

Total 

Number  of  Specimens 

407 

87 

494 

Number  of  Groups 

104 

15 

119 

Average  Standard  Deviation 

.1019 

.0886 

.0995 

Average  Weibull  Shape  Parameter 

5.29 

5.97 

5.41 

Load  Transfer  Element 

Fatigue  Test  Data 

7075 

2024 

Total 

Number  of  Specimens 

141 

182 

323 

Number  of  Groups 

38 

21 

59 

Average  Standard  Deviation 

.1100 

.0927 

.0998 

Average  Weibull  Shape  Parameter 

4.90 

5.06 

5.00 

Full  Scale  Structure 

Fatigue 

Test  Data 

7075 

2024 

Total 

Number  of  Specimens 

170 

73 

243 

Number  of  Groups 

75 

7 

82 

Average  Standard  Deviation 

.0932 

.1056 

.0985 

Average  Weibull  Shape  Parameter 

5.97 

4.87 

5,44 

Note  that  the  7075  and  2024  values  are  generally  about  the  same.  There  is  a 
slightly  larger  scatter  indication  for  7075  element  fatigue  specimens,  but 
thla  is  offset  by  the  smaller  :ctter  indication  for  7075  full  scare  struc¬ 
ture.  The  total  values  do  not  show  sizeable  scatter  difference  among  no-load 
transfer  element  fatigue  specimens,  load  transfer  element  fatigue  specimens, 
and  full  scale  structure.  Fooling  all  of  these  gives  the  following  averages: 

Aluminum  Spectrum  Fatigue  Test  Data 

Total 


Number  of  Specimens  1060 
Number  of  Groups  260 
Average  Standard  Deviation  .0994 
Average  Weibull  Shape  Parameter  5.27 


2.3  Scatter  in  Scatter 

Figure  5  lists  standard  deviations  and  Weibull  shape  parameters  computed 
from  spectrum  fatigue  test  data  for  full  scale  structure.  It  is  hypothesized 
that  all  of  these  data  are  from  the  same  statistical  population  in  tern.s  of 
scatter.  The  range  in  stendard  deviations  is  from  .0415  to  .1326;  the  range 
in  Weibull  shape  parameters  is  from  10.96  to  3.38.  It  is  further  hypothesized 
that  these  data  are  in  the  same  statistical  population  with  the  element  spec¬ 
trum  fatigue  test  data  discussed  in  previous  sections.  These  element  test  data 
exhibit  even  a  wider  spread  in  computed  standard  deviations  and  Weibull  shape 


parameters.  The  question  might  be  aakad  at  this  point,  "If  all  of  these  data 
are  from  the  same  parent  population,  why  don't  they  all  show  a  similar  scatter 
number?"  The  answer  is  that  the  standard  deviations  and  Weibull  shape  para¬ 
meters  computed  from  the  data  are  sample  values  and  as  such  are  actually  random 
variables.  These  random  variables  have  their  own  theoretical  probability  dis¬ 
tributions  and  should  exhibit  scatter  according  to  those  distributions. 

Figure  6  presents  a  comparison  of  the  theoretical  probability  distribution  to 
the  spectrum  fatigue  data  standard  deviations  for  sample  size  n  *  3.  Figure  7 
presents  a  comparison  of  the  theoretical  probability  distribution  to  the 
spectrum  fatigue  data  Weibull  shape  parameters  also  for  sample  size  n  n  3. 

Note  that  both  graphs  indicate  reasonably  good  correlation  between  the  data 
and  the  theoretical  curves. 


2.4  Theoretical  Scatter  Factor.  Ratio  of  Two  Fatigue  Lives 

Every  discussion  on  scatter  must  Include  some  definition  of  "scatter"  and 
relate  it  to  a  parameters  such  as  the  mean,  the  median,  or  some  characteristic 
value.  In  terms  of  airplane  fatigue  life,  there  is  usually  a  full-scale  lab¬ 
oratory  test  article.  The  life  of  this  test  article  is  the  obvious  choice  for 
a  number  about  which  to  center  the  scatter.  However,  the  full-scale  test 
article  life  is  in  fact  a  random  variable.  If  one  considers  a  fleet  of  air¬ 
craft,  the  laboratory  teat  article  can  be  thought  of  as  one  airplane  selected 
at  random  'am  the  total  fleet.  Then  the  scatter  factor  between  the  labora¬ 
tory  life  and  the  service  life  of  another  airplane  picked  at  random  from  the 
fleet  ia  the  "ratio  of  two"  randomly  selected  variates  from  the  same  popula¬ 
tion.  A  scatter  factor  so  defined  being  the  ratio  of  two  statistically 
independent  random  variables  is  itself  a  random  variable.  Its  probability 
distribution  can  generally  be  derived  given  the  distribution  of  the  parent 
population. 


2.4.1  Ratio  of  Two  Log-Norm*!  Variates  -  The  probability  distribution 
for  the  ratio  of  two  statistically  independent  log-norrai  variates  ia  utilized 
to  evaluate  fatigue  scatter  in  Reference  <4).  The  relationship  between  re¬ 
liability  R  (probability  of  no  failure)  and  scatter  factor  8  ia  shown  to  be 

,^/n+i 

F  n 


S  -  10 


G01 


where  n  is  the  number  of  laboratory  teat  articles  and  t  ia  the  normal  variate 
obtainable  from  a  normal  probability  table  such  as  shown  in  Figure  3.  The 
Area  noted  in  Figure  8  correspond*  to  th*  reliability  R.  The  scatter  factor 
for  th*  above  formula  i*  defined  as  the  ratio  of  the  log  mean  life  of  u 
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Figure  6 

Comparison  of  the  Theoretical  and  Observed  Distributions  of  Estimates 
of  the  log  -  Normal  Standard  Deviation  for  Sample  Size  n  *  3 
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laboratory  test  articles  to  the  life  of  an  Individual  airplane  in  service. 
Because  of  the  economics  of  full  scale  airplane  laboratory  fatigue  testing, 
there  is  usually  only  one  test  article  which  means  n  ■  1.  This  is  the  ratio 
of  two  situation  for  which  the  scatter  factor  Is 

S  -  iotcr>/^ 

Using  this  formula  with  a  »  .0994  and  t  ■  2.327  from  Figure  8  gives  a  scatter 
factor  equal  to  2.12  for  99%  reliability. 

2.4.2  Ratio  of  Two  Weibull  Variates  -  The  detailed  mathematical  deriK— 
tion  for  the  ratio  of  two  Weibull  variates  is  included  in  Appendix  hi.  The 
result  of  the  derivation  giving  the  relationships  between  reliability  R  and 
scatter  factor  S  la  as  follows: 


(l^sV 

jher©  a  again  is  the  number  of  laboratory  test  articles.  The  scatter  factor 
for  the  above  formula  is  defined  as  the  ratio  of  the  sample  characteristic 
life  of  n  laboratory  test  articles  to  the  life  of  an  airplane  in  service. 

When  there  is  only  one  laboratory  teat  article,  n  **  1  and  the  formula  reduces 
to  the  ratio  of  two  case  for  which  the  scatter  factor  is 


1 


Using  this  formula  with  a  -  5.27  gives  a  scatter  factor  equal  to  2.39  for  995 
reliability. 

2.5  Comparisons  of  Data  to  Weibull  and  Log-Normal 

As  mentioned  in  Section  2.1,  experimental  evidence  au^c  be  the  primary 
factor  in  evaluating  what  distribution  function  should  be  used.  The  diffi¬ 
culty  has  always  been  in  obtaining  enough  data  to  allow  a  determination 
among  the  various  probability  distributions.  The  data  studied  herein  include 
no-load  transfer  element  specimens,  load  transfer  element  specimens,  and  full 
scale  aircraft  structure  test  results.  A  total  of  IGbO  data  points  were  com¬ 
pared  and  appeared  to  be  from  the  same  statistical  population  in  terns  of 
scatter.  These  test  results  can  be  plotted  in  different  ways  and  compared 
to  the  Weibull  and  log-normal  probability  distribution  functions. 
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The  1060  data  points  are  from  260  groups  where  the  group  sizes  vary 
from  just  two  specimens  to  as  many  as  twenty.  Figure  9  presents  an  example 
group  of  six.  For  the  purpose  of  comparing  theoretical  probability,  distri¬ 
butions  for  the  ratio  of  two  to  experimental  ratios  of  two,  three  statisti¬ 
cally  independent  data  points  can  be  obtained  from  the  example  group  of  six. 
However,  depending  on  the  way  the  results  are  paired,  a  number  of  different 
sets  of  three  can  be  obtained  as  shown  in  Figure  9.  A  total  of  fifteen  diffe¬ 
rent  random  sorts  were  used  to  pair  the  individual  test  results  for  the  three 
types  of  specimens,  i.e.,  no-load  transfer,  load  transfer,  and  full  scale 
structure.  An  additional  five  random  sorts  were  made  for  all  the  data 
combined.  The  theoretical  curves  for  the  Weibull  and  log-normal  ratio  of  two 
are  compared  to  typical  random  sorts  of  the  data  in  Figures  10  through  13. 

It  aopears  that  the  Weibull  provides  the  better  fit. 
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Set me  Factors  for  Five  Osftawu  Random  Sorts  of  a  Gtotip  of  n  *  6 
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CompiiTitan  of  Theofstic^.  PrabaKUty  Ofstritoiitkmsta  Aluminum  Fatigue  Test  Data 
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Comparison  of  Tbeonrtiosf  Prafoahlffty  Distribution*  to  AJumfnum  Fatigue  Test  Data 


Figure  12 

Comparison  of  Theoretical  Probability  Distributions  to  Aluminum  Fatigue  Test  Data 


Figure  13 

Comparison  of  Theoretical  Probability  Distributions  to  Aluminum  Fatigue  Test  Data 


3.  PHASE  II  -  EVALUATION  OF  USAGE  SEVERITY  SCATTER 

3.1  General 

The  study  of  fighter  airplane  usage  data  is  presented  in  this  section* 

It  includes  identifying  scatter  trends  in  usage  severity,  examining  effects 
of  changes  in  operational  usage,  identifying  scatter  trends  in  total  hour 
accumulation,  and  evaluating  probability  distributions  for  "goodness  of  fit" 
with  compiled  data.  Counting  accelerometer  and  total  flight  hour  data  were 
studied  to  quantitatively  identify  variations  in  usage  from  airplane  to  air¬ 
plane  in  a  given  fleet.  In  addition,  VGH  data  were  analyzed  to  determine 
trends  in  speed,  altitude,  and  gross  weight  usage.  These  data,  the  speed/ 
altitude/ gross  weight  information,  are  utilized  in  conjunction  with  counting 
accelerometer  data  to  compute  service  fatigue  damage  at  various  fatigue  criti¬ 
cal  locations  on  the  airplane. 

3.2  F-4  Training  and  Combat  VGH  Data 

A  total  of  8200  hours  of  F-4  Phantom  II  VGH  data  were  analyzed  to  deter¬ 
mine  trends  fcr  airspeed,  altitude,  and  gross  weight  usage  both  in  combat  and 
training  operations.  Histograms  for  these  parameters  were  prepared  for 
maneuvers  exceeding  3.0g's,  4.6g's,  and  6.6g’s  (levels  consistent  with  VGH 
recording  intervals) .  Average  values  for  each  parameter  for  maneuvers 
exceeding  each  of  the  three  load  factors  levels  were  computed  for  both  combat 
and  training  operations.  In  addition,  the  flight  regimes  in  which  the  major¬ 
ity  of  maneuvers  occur  in  both  training  and  combat  were  determined. 

The  histograms  for  maneuvers  exceeding  4.6g's  for  combat  and  for  tracing 
operations  are  presented  in  Figures  14  through  19.  These  histograms 
illustrate  usage  treads  typical  of  those  noted  at  each  of  the  three  load 
factor  levels  examined.  The  summary  of  average  speedo,  altitudes,  and  gross 
weights  for  combat  and  training  operations  for  maneuvers  exceeding  each  of 
the  three  load  factor  levels  is  presented  in  Figure  20.  The  F-4  usage  trends 
apparent  from  all  of  these  data  are  as  follows: 

(a)  The  majority  of  maneuvers  in  both  combat  and  training  operations 
are  pulled  in  a  limited  Mach  number/altitude  regime.  In  bath  types 
of  operation,  the  majority  of  maneuvers  are  pulled  at  between 
350  and  550  knots  and  below  10,000  feet  (the  combat  average 
altitudes  being  slightly  higher  than  those  for- training). 


30 


28  32  36  40  44  4B 

Walfiht- 10001b 

Figure  17 

Histogram  of  Weight  Usage  During  Maneuvers  for  Training  Operations 
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Figuri  18 

Kutogrmn  of  Altitude  Utag*  Durofi  M«ttuv«n  fo*  Trmaine  Op&ttioa* 


35 


3548  hrs  of  VGH  Training  Data 
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Figure  20 

Summary  of  F4  Weigh  t/Albtuda/ Airspeed  U*aga  fu.  Combat  and  Training  Operations- 
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(b)  The  gross  weights  in  combat  for  maneuvers  exceeding  each  of  the 
three  load  factor  levels  are  higher  than  those  for  maneuvers 
executed  in  training  operations. 

(c)  In  both  training  and  combat,  the  average  speed  at  which  maneuvers 
are  pulled  increases  as  the  load  factor  increases. 

The  above  trends  were  not  unexpected.  The  F-4  has  been  used  primarily 
for  air-to-ground  weapon  drops  for  which  low  speed  and  low  altitude  provide 
more  accurate  deliveries.  Thus,  the  majority  of  maneuvers  should  be  expected 
to  be  executed  in  a  limited  speed/altitude  range.  The  slightly  higher 
average  altitude  in  combat  reflects  ground  fire  avoidance.  The  higher  gross 
weights  for  combat  reflect  the  higher  weapon  payload  required  in  actual  com¬ 
bat  service.  Finally,  load  factor  capability  is  dependent  upon  airspeed. 
Higher  airspeeds  are  required  to  pull  High  load  factor  maneuvers.  Thus, 
higher  level  maneuvers  should  he  expected  to  be  executed  at  correspondingly 
higher  average  airspeeds  than  lower  level  maneuvers. 

In  Phase  IV  of  this  study  service  failures  are  compared  to  life  pre¬ 
dictions  based  on  all  the  available  data.  Average  life  predictions  for  each 
key  area  were  derived  from  counting  accelerometer  information  describing  each 
airplane's  load  factor  us&gu*  and  from  V<S4  data  which  helped  define  stress 
versus  load  factor  relationships.  Load  factor  exceedance  information  was 
eaneutial  for  making,  a  .'.treasonable  estimate  of  m  airplanes  expected  fatigue 
life,  but  an  accurst# ' required  even  sore  information.  Speed, 
altitude,  and  gross  weight  4te  the  important  parameters  chat  influence 
internal  structural  leads.  The  ¥S8  information  noted  above  thus  provided 
the  prissaty  data  necessary  to  define  the  relationships  between  stress  and 
maneuvering  load /fnctor.  For  example,  Figure  21  shows  bow  stress  in  the  F-4 
wing  main  torque  box  lover  skin  varies  with  spued  ana  altitude.  From  the  VGH 
trend®  it  may  be  noted  that  the  majority  of  maneuvers  are  pulled  in  a  flight 
regime  in  which  wing  bending  assent  per  g  is  about  85S  of  that  at  the  critical 
design  condition. 

3.5  F-4  'Counting  Accelerometer  Data 

3.3*1  Keshcd  for  Evaluating  I'sago  Scatter  -  The  examination  of  scatter 
trends  in  load  factor  counting  accelerometer  data  was  a  completely  automated 
operation.  A  computer  program  developed  for  this  purpose  calculated  scatter 
iu  load  factor  count  accumulation  versus  effective  flight  hours  for  all  air¬ 
planes  reporting  counting  accelerometer  data.  The  scatter  computed  represents 


40.000  Pounds  Gross  Weight 
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Figure  21 

Typical  F4  Wing.  Station  Comparison  of  Flight  Condition  end  Laboratory  Test  Bending  Moments  Par  g 


the  ratio  between  the  load  factor  counts  accumulated  by  a  given  airplane  and 
the  load  factor  counts  which  would  have  been  accumulated  by  that  airplane  if 
it  had  been  operating  at  the  fleet  average  accumulation  rate. 

Counting  accelerometer  data  for  F-4  airplanes  are  stored  on  magnetic 
tape,  the  computer  program  uaed  was  designed  to  read  the  data  directly  off 
the  storage  tapes.  For  each  airplane  reporting  counting  accelerometer  data,, 
the  reporting  dates,  comets  accumulated,  and  aircraft  flight  hours  were  read 
directly  from  the  storage  tapes.  Fleet  average  daily  count  per  hour  rates 
were  then  computed  using  the  data  from  all  aircraft  in  the  fleet.  Following 
this,  scatter  for  a  given  load  factor  exceedance  level  was  computed  at  each 
of  the  reporting  dates  listed  for  each  airplane.  These  scatter  values  were 
determined  by  taking  the  ratio  of  the  actual  counts  accumulated  through  the 
reporting  date  divided  by  the  estimated  accumulated  counts  (estimated  counts 
were  represented  by  the  summation  of  the  products  of  the  flight  hours  accumu¬ 
lated  on  given  days  by  the  airplane  multiplied  by  the  applicable  fleet  average 
daily  count  per  hour  rates).  Finally,  the  usage  scatter  data  were  scanned 
and  evaluated,  and  information  necessary  to  construct  histograms  shoving  the 
distribution  of  scatter  at  various  numbers  of  effective  flight  hours  was 
produced,  effective  hours  (estimated  counts  divided  by  the  fleet  overall 
average  counts  r«r  hour)  were  used  in  this  study  so  that  usage  scatter  from 
aircraft  operating  in  different  time  periods  could  be  compared  directly.  The 
significance  of  effective  hours  is  discussed  in  the  following  suction . 

3.5.2  Effective  Flight  flours  -  F-4  operational  usage  has  varied  signi¬ 
ficantly  since  the  airplane  was  introduced  into  service.  As  discussed  la 
Section  3.3.5,  monthly  load  factor  count  accumulation  rates  liave  shown  sub¬ 
stantial  fluctuations  over  the  years.  An  circraft  entering  service  early 
in  the  F-4  program  would  be  expected  to  have  accumulated  considerably  fever 
coitnta  in  the  same  number  of  hours  than  would  an  airplane  entering  service 
in  the  last  few  years.  Therefore,  if  aircraft  with  the  same  actual  number  of 
hours  were  directly  compared ,  the  estimated  number  of  counts  would  in  general 
vary  from  one  airplane  to  another.  The  estimated  number  of  counts  is  repre¬ 
sentative  of  tiie  expected  number  of  counts  or  the  average,  hence,  a  compari¬ 
son  based  on  actual  number  of  hours  would  consist  oi  airplanes  with  different 
averages  in  the  same  group,  and  scatter  trends  based  on  such  a  comparison 
would  not  really  be  meaningful.  On  the  other  hand,  observed  scatter  trends 
would  have  more  significance  if  based  on  comparisons  of  airplanes  with  the 
same  average  or  estimated  number  of  counts. 
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Effective  hours  ia  equal  to  estimated  number  of  counts  divided  by  the 
fleet  overall  average  counts  per  hour.  Then  comparing  airplanes  baaed  on  the 
number  of  effective  hours  they  have  accumulated  would  place  airplanes  in  the 
same  group  that  have  the  same  estimated  number  of  counts  or,  in  other  words, 
the  same  average.  The  main  difference  between  actual  hours  and  effective 
hours  is  that  aircraft  which  flew  predominantly  early  in  the  history  of  the 
F-4  program  have  fewer  effective  hours  than  actual  hours  since  they  flew  at 
a  time  when  count  accumulation  rates  were  relatively  low;  whereas,  aircraft 
entering  service  in  the  more  recent  history  have  more  effective  hours  than 
actual  hours  because  they  were  flying  in  a  period  in  which  usage  was  more 
severe  than  the  overall  average. 

3.3,3  Usfrge  Scatter  Trends  Versus  Hours  -  Computer  runs  were  made  to 
evaluate  4g  and  tc  usage  scatter  exhibited  by  Air  Force  and  by  Navy /Marine 
airplanes.  The  d*'ca  shewing  the  distribution  of  scatter  at  selected  flight 
hour  intervals  was  generated  during  those  computer  runs.  This  information 
is  presented  in  Figures  22  through  25.  The  outputs  show  the  number  of  air¬ 
planes  with  scatter  values  in  given  ranges  versus  effective  hours.  It  should 
be  noted  that  in  each  200  hour  interval,  only  one  scatter  value  is  counted 
for  each  airplane.  The  scatter  value  included  la  the  one  at  the  number  of 
hours  closest  to  the  midpoint  of  the  interval.  Thus,  the  column  labeled  "401 
to  600  hours’*  yields  the  distribution  of  the  number  of  airplaaaa  with  various 
scatter  values  at  500  hours. 

Figures  22  through  25  indicate  that  the  dispersion  of  usage  scatter 
points  decreases  significantly  with  increasing  hours.  Its  is  trend  may  he 
explained  intuitively  on  the  boats  that  the  lunger  an  airplane  ia  flying,  the 
note  likely  that  ft  will  have  been  subjected  to  a  variety  of  usages  and  its 
repeated  loading  will  have  '’averaged  out,**  The  trend  is  also  in  qualitative 
agreement  with  probability  theory  which  states  that  the  ratio  of  the  standard 
deviation  to  the  wean  fare  number  of  occurrence's  in  a  given  tire  interval  is 
inversely  proportional  to  the  square  root  of  the  time  interval  length.  Al¬ 
though  the  trend  toward  decreasing  scatter  agrees  qualitatively  wife  this 
theory,  sa  examination  of  the  data  has  indicated  that  the  data  and  the  theory 
are  not  in  exact  quantitative  agreeaauc.- 

0s«f o  severity  scatter  aaeoa  and  variances  versus  effective  hours  for 
Air  Force  and  Navy  Marius  4g  and  maneuver  accumulations  are  presented  in 
Figures  26  through  .29*  These  data  show  how  usage  severity  scatter  variances 
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Figure  26 

4g  Usage  Severity  Scatter  Mean  and  Variance  vs  Effective  Hours  for  Air  Force  Aircraft 
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Figure  27 

4g  Usage  Severity  Scatter  Mean  and  Variance  vs  Effective  Hours  for  Navy/Marine  Aircraft 
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change  with  increasing  tine.  As  may  be  noted,  the  variances  decrease  with 
increasing  time  but  not  in  proportion  to  the  inverse  of  time  as  is  suggested 
by  probability  theory.  In  Reference  <5),  it  is  indicated  that  for  such  a 
relation  to  hold,  the  accumulation  of  flight  maneuvers  by  airplanes  in  a 
fleet  would  have  to  be  a  stationary  random  process.  One  of  the  requirements 
for  such  a  process  is  that  the  mean  frequency  of  occurrence  of  maneuvers  must 
be  constant  with  respect  to  time.  As  has  been  shown,  this  has  not  been  the 
case  with  F-4  service  usage;  operational  usage  has  changed  and  fleet  average 
maneuver  frequencies  of  occurrence  have  fluctuated  greatly.  Usage  severity 
scatter  variances  could  not,  therefore,  be  expected  to  decreese  exactly  in 
proportion  to  the  inverse  of  time.  However,  the  fact  that  scatter  variances 
are  shown  to  decrease  and  not  Increase  with  time  can  be  used.  The  variance 
at  1500  hours  is  calculable  from  the  available  data.  The  question  ia  vhat 
will  be  the  variance  at  a  design  life  of  aay  4000  hours.  Since  variances  are 
shown  not  to  increase  but  decrease,  although  not  inversely  proportional  to 
the  square  root  of  time,  it  would  certainly  be  reasonable  to  as sum  that  the 
variance  at  4000  hours  is  not  larger  than  the  value  at  1500  hours. 

3, 3. 3,1  Graphical  Presentation  of  Usage  Severity  Scatter  -  Plata  of  4g 
and  6g  usage  severity  scatter  versus  flight  hours  were  generated  for  both  Air 
Force  airplanes  and  Navy /Hat ice  airplanes.  The  plots,  presented  in  Figures 
30  through  33,  illustrate  graphically  the  decrease  in  scatter  with  increasing 
hours.  As  aay  be  noted,  in  all  cases,  the  dispersion  of  uaaga  severity 
scatter  is  shown  to  decrease  significantly  with  iueraaaiog  flight  hours. 

As  has  been  indicated  previously,  this  trend  may  be  explained  intuitively 
on  the  basis  that  the  longer  ^n  airplane  is  flying,  the  more  likely  that  it 
will  have  been  subjected  to  a  variety  of  usages  and  Its  repeated  loading  will 
have  averaged  out. 

In  accordance  with  this  theory  of  repeated  load  averaging  over  long 
period*,  large  scatter  factor®  can  be  expected  only  early  in  an  sitplane*® 
life  when  the  number*  of  actual  and  estimated  counts  are  relatively  low. 

There  are  captions  vo  this  rule  however.  The  plot  of  6g  usage  scatter  for 
Navy /Marine  airplane*  ahowu  In  Figure  33  contains  85  data  point*  for  F-4B 
Bureau  No.  150492  plotted  at  between  62  and  2482  flight  hour®.  The  last  of 
these,  a  scatter  factor  of  9.31  at  2482  flight  hours  (see  Figure  33),  ia 
unusually  high  for  that  tuny  flight  hour®  and  i®  obviously  inconsistent  with 
the  data  from  the  other  aircraft  in  the  fleet.  Prior  to  November  1966,  this 
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Figure  3  T 

4g  Usage  Severity  Scatter  vs  Flight  Hours  for  Navy /Marine  F-4  Aircraft 
Reporting  Counting  Accelerometer  Data 

UQ39  Aircraft  958,377  Flight  Hours  of  Data  Through  March  1972,  29,507  Data  Points  Plotted) 


Gg  Usage  Severity  Scatter  vs  Flight  Hours  for  Navy/Marina  F-4  Aircraft 
Reporting  Counting  Accelerometer  Data 

(947  Aircraft,  836,616  Flight  Hours  of  Data  Through  March  1972,  25,104  Data  Points  Plotted) 


airplane  had  accumulated  approximately  1200  flight  hours  and  had  been  engaged 
in  operations  typical  of  an  average  airplane  in  the  fleet.  Since  that  time, 
however,  the  airplane  has  been  engaged  in  development:  of  weapon  delivery  tech¬ 
niques.  Since  transferring  into  such  operations,  the  aircraft's  6g  exceedance 
accumulation  rate  has  continually  increased.  Recently,  the  airplane's  average 
usage  lias  been  12  times  more  severe  than  that  of  a  typical  airplane  in  the 
fleet.  This  represents  an  activity  which  is  not  at  all  comparable  to  typical 
operations  in  the  overall  fleet  of  Navy /Marine  aircraft.  The  flight  hours 
accumulated  by  Bureau  No.  150492  are  thus  considered  to  be  from  unusual  opera¬ 
tions,  and  as  such,  the  data  from  this  usage  are  not  included  in  any  scatter 
evaluations  herein. 

It  is  frequently  found  that  unusual  usage  such  as  that  experienced  by 
Buteau  No.  150492  leads  to  fatigue  problems.  Failure  of  an  airplane  such  as 
this  in  unknown  or  unrecognized  special  operations  would  cast  a  shadow  of 
doubt  au  the  structural  integrity  of  the  total  fleet  of  aircraft.  By  having 
counting  accelerometers  installed  in  each  aircraft,  however,  these  operations 
can  be  recognised  and  special  periodic  inspections  of  the  aircraft  involved 
can  be  conducted.  This  should  preclude  any  major  problems.  In  addition, 
scatter  can  be  controlled,  if  deemed  advisable,  by  transferring  airplanes  in 
and  out  of  unusual  activities. 

3.3.4  Evaluation  of  Probability  Distribution  -  The  statistical  repre¬ 
sentation  of  aircraft  flight  maneuver  accumulation  is  discussed  in  detail  ir 
tofereace  (5).  It  is  shown  that  the  distribution  of  the  numbers  of  maneuvers 


accumulated  by  individual  aircraft  in  a  ,;leet  may  be  described  using  the 


negative  binomial  distribution.  Ussgc  severity  scatter,  as  used  herein,  is 
defined  as  the  ratio  of  the  actual  number  of  load  factor  counts  accumulated 
by  an  airplane  divided  by  the  number  of  counts  which  would  have  been  accumu¬ 
lated  by  that  airplane  had  it  been  operating  at  the  fleet  average  maneuver 
accumulation  rate.  Since  this  ratio  is  a  measure  of  an  individual  airplane's 
accumulation  of  maneuvers,  the  distribution  of  usage  severity  scatter  should 
also  follow  the  negative  binomial  distribution. 

The  negative  binomial  distribution,  defined  by  the  probability  density 
function. 


r(s^m) 

Us+l) 
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where : 


z  *  usage  severity  scatter  ■  actual  counts /estimated  counts 

o  =  mean  value  z 
2 

b  ■  m/  (c  -  m) ,  and 
o  -  standard  deviation  of  z 

has  been  evaluated  for  "goodness  of  fit"  with  the  usage  severity  scatter  data 
generated  for  both  Air  Force  and  Navy /Marine  airplanes.  The  evaluation  has 
shown  definite  correlation  between  the  theoretical  distribution  and  the  data. 
As  an  example,  Figures  34  through  36  show  comparisons  between  the  negative 
binomial  distribution  and  the  4g  scatter  exhibited  by  Air  Force  airplanes 
having  accumulated  100 ,  500,  and  1500  effective  hours.  These  comparisons 
show  that  the  scatter  in  usage  severity  at  any  given  number  of  effective 
hours  can  be  represented  using  the  negative  binomial  distribution. 

3.3.5  Effects  of  Changes  in  Operational  Usage  -  The  definition  of  a 
particular  airplane’s  usage  may  change  completely  from  the  design  mission  due 
to  the  character  of  the  actual  theatre  of  operations.  This  has  been  clearly 
illustrated  by  the  usage  of  the  F-4  airplane  in  service.  As  is  well  known, 
the  F-4  has  been  developed  into  a  highly  diversified  weapons  system.  It  has 
been  utilized  as  an  interceptor,  a  fighter-bomber,  and  in  numerous  ground 
support  operations.  The  differences  between  the  usages  in  these  roles  has 
been  drastic.  The  SEA  conflict  has  likewise  had  a  significant  effect  on 
operational  usage  considerations.  In  addition  to  increased  maneuver  frequency 
resulting  from  deployment  in  conventional  bombing  operations,  the  utilization 
rates  and  the  actual  flight  gross  weights  went  beyond  initial  expectations. 
Instead  of  being  flown  20  tc  30  hours  per  month,  the  rate  for  combat  squadron 
activities  increased  to  60  to  70  hours  per  month.  Also,  as  illustrated  by 
the  VGH  data  presented  in  Section  3.2,  the  average  gross  weight  In  combat 
operations  increased  to  a  level  significantly  above  the  design  gross  weight 
of  37,500  pounds.  None  of  these  changes  could  have  been  predicted  statisti¬ 
cally  during  the  design  phase  of  the  airplane. 

Air  Force  and  Navy /Marine  4g  and  6g  fleet  average  exceedance  accumula¬ 
tion  rates  versus  time  are  shown  in  Figure  37.  These  plots  show  how  F-4 
usage  has  varied.  As  noted  previously,  the  usage  scatter  factors  computed 
herein  represent  the  scatter  between  the  exceedances  accumulated  by  a  given 
airplane  and  the  exceedances  which  would  have  been  accumulated  by  that  air¬ 
plane  if  it  had  been  operating  at  the  fleet  average  exceedance  accumulation 
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Figure  34 

Comparison  of  Air  Force  4g  Usage  Severity  Scatter  tc  the  Negative  Binomial  Distribution  (100  Effective  Hours) 


Figure  35 

Comparison  of  Air  Fores  4g  Usage  Severity  Scatter  to  the  Negative  Binomial  Distribution  (500  Effective  Hours) 


Figure  36 

Comparison  of  Air  Force  4g  Usage  Severity  Scatter  to  the  Negative  Binomial  Distribution  (1500  Effective  Hours) 


10,000 


Figure  37 

Load  Factor  Usage  Variation  with  Time  for  Air  Force  and  Navy/Marine  Aircraft 


rate.  Scatter  factors  were  computed  at  each  reporting  date.  This  factor  is 
the  ratio  of  the  actual  counts  accumulated  through  that  date  divided  by  the 
estimated  accumulated  counts.  Both  the  actual  counts  accumulated  and  the 
estimated  counts  reflect  changes  in  operational  usage.  The  actual  counts, 
having  been  recorded  in  service  operations,  necessarily  include  changes  in 
operational  usage.  The  estimated  counts,  being  computed  from  fleet  average 
daily  count  per  hour  rates,  reflect  changes  in  operational  usage  occurring 
during  the  period  in  which  the  aircraft  was  in  service.  Thus,  if  an  aircraft 
had  been  operating  in  a  period  in  which  load  factor  frequencies  continually 
increased,  both  the  actual  counts  and  the  estimated  counts  would  exhibit 
corresponding  increases.  By  taking  the  ratio  of  these  counts,  the  effects  of 
changes  in  operational  usage  have  been  removed,  and  the  scatter  is  simply  that 
between  a  given  airplane  and  an  average  airplane  in  the  fleet. 

Although  it  can  have  a  significant  effect  on  service  fatigue  life,  a 
change  in  operational  usage  is  not  a  statistical  variable  and  cannot,  there¬ 
fore,  be  approached  on  a  statistical  basis.  In  arriving  at  a  total  scatter 
factor,  the  scatter  in  fatigue  test  results  and  in  usage  can  be  handled 
statistically,  but  it  is  necessary  to  make  a  "best  guess"  as  to  how  the 
airplane  will  be  used  in  service.  Changes  in  operational  usage,  although 
quite  important,  cannot  be  anticipated  or  accurately  estimated  during  the 
aircraft's  design  phase. 

3. 4  Scatter  in  Total  Hour  Accumulation  on  Individual  Aircraft 

In  addition  to  developing  data  and  procedures  for  computing  a  design 
scatter  factor  for  a  preselected  reliability  level,  it  is  also  important  to 
determine  what  should  be  the  der  gn  life  of  an  airplane  in  terms  of  hours  or 
years  of  service  operations.  Data  reflecting  scatter  in  the  total  hour  acc¬ 
umulation  on  individual  F-4  aircraft  have  thus  been  investigated  to  obtain 
information  which  can  be  used  to  aid  in  establishing  an  aircraft's  design 
life.  The  study  performed  and  the  information  obtained  are  described  in  the 
following  paragraphs. 

3.4.1  Method  for  Evaluating  Hour  Accumulation  Scatter  -  In  this  study, 
the  scatter  numbers  computed  represented  the  scatter  between  the  flight  hours 
accumulated  by  a  given  airplane  and  the  flight  hours  which  would  have  been 
accumulated  by  that  airplane  if  it  had  been  operating  at  the  fleet  average 
hour  accumulation  rate.  Total  flight  hour  data  for  F-4  airplanes  are  stored 
on  the  magnetic  tapes  which  contain  the  F-4  counting  accelerometer  data. 


61 


Aircraft  reporting  dates  and  hours  accumulated  were  read  directly  from  the 
storage  tapes.  A  daily  fleet  average  hour  accumulation  rate  was  first  esta¬ 
blished  by  determining  the  total  flight  hours  accumulated  on  any  given  day 
and  dividing  by  the  total  number  of  aircraft  in  the  service  inventory  on  that 
day.  Aircraft  known  to  not  have  flown  were  included  in  the  inventory  count 
in  order  to  account  for  expected  aircraft  down-time.  Scatter  was  then  com¬ 
puted  at  each  reporting  date  listed  for  every  airplane.  Estimated  hours  for 
individual  aircraft  were  computed  by  accumulating  the  fleet  daily  average  hour 
accumulation  rates  applicable  for  periods  between  reporting  dates.  The 
scatter  at  each  reporting  date  was  then  computed  by  taking  the  ratio  of  the 
actual  hours  accumulated  through  that  date  divided  by  the  estimated  accumu¬ 
lated  hours.  Finally,  the  scatter  data  were  scanned  and  evaluated  and  infor¬ 
mation  necessary  to  construct  histograms  showing  the  distribution  of  scatter 
at  selected  intervals  of  years  was  produced. 

3.4.2  Hour  Accumulation  Scatter  Versus  Years  -  Computer  runs  were  made 
to  evaluate  the  hour  accumulation  scatter  exhibited  by  Air  Force  airplanes 
and  by  Navy/Marine  airplanes.  In  the  Air  Force  run,  1,701,964  flight  hours 
of  data  accumulated  by  1232  Air  Force  airplanes  were  examined.  In  the  Navy/ 
Marine  run,  1,372,702  flight  hours  of  data  aceumulted  by  1114  Navy /Marine 
airplanes  were  examined.  As  was  the  case  in  the  usage  severity  scatter  study, 
the  scatter  in  flight  hours  accumulated  was  also  noted  to  decrease  signifi¬ 
cantly  with  increasing  time.  The  trend,  exhibited  in  both  the  Air  Force  and 
Navy/Marine  runs,  indicates  that  flight  hour  usage  also  averages  out  and  that 
aircraft  placed  in  the  fleet  at  the  srme  time  will  tend  to  accumulate  similar 
numbers  of  hours  over  a  period  of  years. 

3.4.3  Histograms  of  Hour  Accumulation  Scatter  at  Selected  Intervals  - 
The  information  necessary  to  construct  histograms  showing  the  distribution  of 
hour  accumulation  scatter  at  selected  year  intervals  was  generated  during  the 
computer  runs.  This  information  for  Air  Force  airplanes  and  for  Navy /Marine 
airplanes  is  presented  in  Figures  38  and  39.  It  should  be  noted  that  in  eacn 
2  year  interval,  only  one  scatter  value  is  counted  for  each  airplane.  The 
scatter  value  included  is  the  one  at  the  number  of  years  closest  to  the  mid¬ 
point  of  the  interval.  The  columns  thus  yield  the  distribution  of  the  number 
of  airplanes  with  various  scatter  values  at  the  number  of  years  at  the  inter¬ 
val  midpoints.  The  histograms  plotted  from  the  data  are  shown  in  Figures 

40  and  41.  Histograms  for  the  data  from  the  Air  Force  run  show  the 
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1,701,964  Total  Hours  of  Data 
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Figure  33 

Flight  Hour  Scatter  Distribution  for  1232  Air  Force  Airplanes 
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1,372,702  local  Hours  of  Data 


Flight  * 

Total  Years  in 

Service 

Hour 

6 

■ 

4.01 

6.01 

12.01 

Scatter 

to 

to 

to 

1KSH 

to 

to 

2.00 

6.00 

8.00 

BLULJISi 

14.00 

0.0  -  0.2 

7 

2 

0.2  -  0.4 

39 

8 

3 

2 

i 

0.4  -  0.6 

68 

40 

31 

19 

9 

0.6  -  0.8 

80 

100 

76 

54 

41 

2 

0.8  -  1.0 

131 

208 

166 

94 

57 

7 

1.0  -  1.2 

211 

204 

202 

99 

38 

3 

1.2  -  1.4 

224 

129 

86 

37 

7 

1.4  -  1.6 

157 

59 

25 

10 

1.6  -  1.8 

95 

26 

5 

1.8  -  2.0 

50 

4 

2 

2.0  -  2,2 

24 

1 

2.2  -  2.4 

2.4  -  2.6 

2,6  -  2.8 

2,8  -  3.0 

5 

3.0  -  3.2 

0 

3.2  -  3.4 

1 

Total  Humber 
of  Aircraft 

1114 

781 

. . 

m 

315 

153 

12 

0 

*  Flight  Hour  Scatter  «  Actual  Flig '.*w  .loura/Eatitsated  Flight  Hours 


Figure  39 

Flight  Hour  Scatter  Diitti  button  tor  1 1 14  Wavy/Mafina  Airplanes 
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Hour  Scatter 


Right  Hour  Scatter  I Aetual/£nl««tod  Plight  Houral 


Figure  40 

Percentage  of  Air  Force  Aircraft  with  Given  Flight  Hour  Scatter 
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Figure  41 

Ptrcgfltagu  of  Navy/Manna  Aircraft  with  Givtn  Fii^it  Hour  Scatter 
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distribution  of  scatter  at  1,  3,  5,  and  7  years.  Histograms  for  the  data 
from  the  Navy /Marine  run  show  the  distribution  of  scatter  at  1,  3,  5,  and  9 
years. 

3.4.4  Effects  of  Operations  Changes  on  Hour  Accumulation  Scatter  -  Air 
Force  and  Navy/Marine  fleet  average  hour  accumulation  rates  versus  time  are 
shown  in  Figure  42.  These  curves  show  how  F-4  hour  accumulation  rates  have 
varied.  As  previously  indicated,  the  hour  accumulation  scatter  values  computed 
herein  represent  the  scatter  between  the  hours  accumulated  by  a  given  airplane 
and  the  hours  which  would  have  been  accumulated  by  that  airplane  if  it  had 
been  operating  at  the  fleet  average  hour  accumulation  rate.  Both  the  actual 
hours  and  the  estimated  hours  used  in  computing  the  scatter  values  include 
operational  usage  change  effects.  The  actual  hours  necessarily  reflect 
changes  in  operational  usage.  Estimated  hours,  being  computed  from  daily 
hour  accumulation  rates,  reflect  changes  in  operational  usage  occurring  dur¬ 
ing  the  period  in  which  the  aircraft  was  in  service.  By  taking  the  ratio,  the 
effects  of  operational  usage  changes  have  been  removed,  and  the  scatter  is 
simply  that  between  a  given  airplane  and  an  average  airplane  operating  in  the 
fleet  during  the  same  period.  Aa  has  been  noted  previously,  this  scatter 
decreases  with  increasing  time.  Since  usage  changes  have  been  removed,  this 
trend  of  decreasing  scatter  indicate#  that  aircraft  placed  in  the  fleet  at 
the  same  time  will  accumulate  similar  numbers  of  hours  during  their  service 
tours.  In  addition  to  this,  however,  since  F-4  flight  hour  accumulation 
rates  indicate  that  monthly  hour  usage  is  relatively  invariant,  the  data  would 
also  indicate  Chat  all  aircraft  in  a  fleet  may  be  expected  to  accumulate  a 
similar  number  of  hours  after  a  like  number  of  years  in  service. 
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4.  PHASE  III  “  DOCUMENTATION  OF  LABORATORY  TEST  RESULTS 

4.1  General 

F-4  full  9cale  laboratory  fatigue  testing  has  been  conducted  in  a  series 
of  test  programs  involving  five  complete  wing-center  fuselage  fatigue  test 
articles  and  two  half-wings.  The  fatigue  characteristics  of  five  fatigue 
critical  key  areas  baaeO  on  this  laboratory  testing  were  documented  in  Phase 
III  of  this  program.  Efforts  in  this  area  included  the  compilation  of  spec¬ 
trum  hours  to  failure  and  the  determination  of  the  crack  growth  characteris¬ 
tics  for  each  key  area.  The  key  areas  studied  were: 

(1)  the  wing  main  torque  box  lower  skin, 

(2)  the  outer  wing  lower  akin, 

(3)  the  F.S.  303  bulkhead, 

(4)  dte  wing  main  torque  box  upper  skin,  and 

(5)  the  lower  longeron  dog  bone  fitting. 

4.2  P-4  Fatigue  Test  Program  History 

P-4  full  scale  fatigue  testing  has  been  conducted  in  a  saries  of  test 
programs  that  have  included  seven  full  scale  teat  articles.  The  test  articles 
that  have  been  utilised  are  depicted  in  Figure  2.  The  test  programs,  which 
were  performed  to  substantiate  the  airframe  fatigue  life  and  to  define  fatique 
critical  areas,  ssr  described  briefly  in  subsequent  paragraphs.  This  test 
history  provides  the  information  necessary  to  identify  the  laboratory  failures 
referred  to  in  this  atudy. 

Block  1  Test  Program  -  The  tear,  article  consisted  of  *  Center  fuselage 
and  cecniete  wing  assembly  representative  of  Che  Ship  1  through  95  configura¬ 
tion.  Test  loads  represented  the  critical  loading  condition  for  a  flight 
gross  weight  of  34,500  pounds  and  veto  applied  in  accordance  with  spectrum  A 
of  MIt-A-6866.  This  resulted  in  the  application  of  a  test  loading  spectrum 
which  was  considerably  more  severe  chan  the  original  design  spectrum.  Failure 
occurred  at  the  end  oi  the  fifth  load  program  when  both  R/H  and  L/H  wings 
failed  simultaneously  at  the  wing  root. .  The  failures  were  precipitated  by 
repeated  buckling  of  the  wing  carry  through  auxiliary  bass  shear  Webs  at  high 
load  levels.  Modifications  to  the  auxiliary  beam  ve re  developed  and  refined 
in  subsequent  jig  testing.  A  retrofit  fix  was  incorporated  in  Ships  1  through 
40  and  redesigned  auxiliary  beams  were  installed  effective  Ship  41  and  up. 

Revision  uf  Test  Loads  Repeated  toad  Spectrun  -  A  flight  load 
measurement  program  indicated  that  the  analytically  determined  loads  used  in 
the  Block  1  test  were  conservative.  The  test  loading  distribution  was  revised 
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to  reflect  the  results  of  flight  teat  measurements,  but  concurrently,  the  mag¬ 
nitude  of  test  loads  was  increased  to  reflect  an  increase  in  gross  weight  from 
34,500  pounds  to  40,000  pounds.  The  combination  of  these  changes  increased 
test  loads  by  a  factor  of  1.039  above  the  Block  1  fatigue  test  leading  and 
rest, ited  in  loads  representative  of  flight  loads  at  critical  design  speed  and 
altitude  (Mach  1.1  at  25,000  ft.).  In  addition  to  the  test  loads  revision, 
the  load  spectrum  was  changed  j  Spectrum  A  of  MIL-A-8866  to  what  has  come 
to  be  known  as  the  F-4  test  sptwtrum  (the  list  of  cycles  applied  in  the  F-4 
test  spectrum,  which  was  utilized  in  all  subsequent  fatigue  testing,  is  pre¬ 
sented  in  Figure  43).  This  combi,  ation  of  changes  resulted  in  a  test  spectrum 
which  was  also  considerably  more  severe  than  the  original  design  spectrum. 

Block  1  R/H  Remnant  Wing  Test  -  The  R/H  remnant  wing  from  the  Block  1 
fatigue  test  article  was  jig  supported  (critical  loads  simulated  outboard  of 
B.L.  80)  and  spectrum  testing  was  continued  using  the  F-4  test  spectrum  and 
loads  (40,000  pounds  gross  weight,  Mach  1.1  at  25,000  ft.).  Failure  occurred 
in  the  wing  main  torque  box  lower  skin  at  B.L.  100  (fatigue  critical  area 
shown  in  Figure  44)  after  4G0  additional  spectrum  hours. 

Block  1  L/H  Remnant  Wing  Test  -  A  lower  wing  skin  reinforcing  strap 
(retrofit  fix  for  Ships  1  through  95)  was  installed  on  the  L/H  remnant  wing 
from  the  Block  1  fatigue  test  article  and  testing  was  continued  as  in  the 
R/H  remnant  wing  test.  Testing  was  discontinued  after  9300  spectrum  hours 
without  a  catastrophic  fatigue  failure. 

Block  6  Test  Program  -  The.  test  article  consisted  of  the  Block  1  test 
fuselage  and  a  new  1  through  95  wing  with  redesigned  auxiliary  beams.  Test 
loads  were  applied  in  accordance  with  the  F~4  test  spectrum  and  loads 
(40,000  pounds  gross  weight,  Mach  1.1  at  25,000  ft.).  Four  hundred  spectrum 
hours  were  applied  prior  to  installation  of  the  lower  wing  skin  reinforcing 
strap  to  simulate  possible  service  history  experienced  by  airplanes  prior  to 
the  reinforcement.  Testing  was  discontinued  after  the  completion  of  4200 
spectrum  hours  without  a  catastrophic  fatigue  failure. 

Block  8  Test  Program  -  The  test  article  consisted  of  a  new  center  fuse¬ 
lage  and  a  wing  representative  of  Ship  No.  96  and  up  (primary  modifications 
consisting  of  Increased  strength  in  lower  skin).  Test  loads  were  applied  in 
accordance  with  the  F-4  test  spectrum  and  loads  (40,000  poundp  gross  weight, 
Mach  1.1  at  25,000  ft,).  Failure  of  both  wings  occurred  simultaneously  in 
the  fatigue  critical  area  at  B.L.  100  (same  location  as  R/H  remnant  wing 
failure  as  shown  in  Figure  44)  after  2700  spectrum  hours. 
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List  of  Cycles  Applied  in  the  F4  Test  Spectrum 
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Figure  44 

Failure  Location  -  R/H  Remnant  Wing  and  Bfofik  8  Test  Articte 


ECP-613  Test.  Program  -  This  test  program  was  initiated  to  demonstrate 
that  the  F-4  test  life  could  be  extended  from  2700  to  6000  spectrum  hours  by 
installing  Taper-Lok  fasteners  in  the  lower  wing  skin  fatigue  critical  area. 
The  test  articles  consisted  of  two  full  scale  aircraft: 

(1)  F-4B  Test  Article  -  Representative  of  early  aircraft  with  7079-T6 
spars,  ribs,  and  7075-T651  bulkheads,  and  with  conventional 
fasteners  installed  in  the  fatigue  critical  area  in  production. 

(2)  F-4J  Test  Article  -  Representative  of  higher  effectivity  aircraft 
with  7075-T73  spars,  ribs  and  bulkheads.  R/H  wing  representative 
of  aircraft  with  conventional  fasteners  installed  in  the  fatigue 
critical  area  in  production  and  the  L/H  wing  representative  of 
aircraft  with  Taper-Lok  fasteners  installed  in  the  critical  area 
in  production. 

Test  loads  were  applied  in  accordance  with  the  F-4  test  spectrum  and  loads 
(40,000  pounds  gross  weight,  Mach  1.1  at  25,000  ft.).  After  approximately 
1700  spectrum  hours,  Taper-Lok  fasteners  were  installed  in  ths  fatigue  criti¬ 
cal  area  of  the  F-4B  test  article  wings  and  in  the  critiJu-.  area  of  the  R/H 
wing  of  the  F-4J  test  article  (1700  spectrum  hours  applied  to  simulate 
possible  service  history  experienced  by  airplanes  prior  to  incorporation  of 
Taper-Lok  retrofit  fix) .  Testing  of  both  test  articles  was  discontinued 
following  the  completion  of  6000  spectrum  hours  without  catastrophic  failure. 

FSCP  46  Extension  -  This  testing  was  conducted  to  determine  whether  the 
Taper-Lok  installation  could  extend  the  test  life  of  the  F-4B  test  article 
used  in  the  ECP-613  test  program  to  8000  spectrum  hours.  Cycling  of  the  test 
article  was  continued  in  the  same  manner  as  in  the  ECP-613  test  program  ex¬ 
cept  cycles  in  the  35%,  45%,  and  55%  load  levels  of  the  F-4  test  spectrum 
were  not  applied.  Cycling  of  the  test  article  was  discontinued  following  the 
completion  of  8000  spectrum  hours  without  catastrophic  failure. 

FSCP  60R1  Extension  -  This  testing  was  conducted  to  determine  whether 
the  test  life  of  the  F-4B  test  article  utilized  in  ECP-613  and  FSCP-46  test¬ 
ing  could  be  extended  to  12000  spectrum  hours.  Cycling  was  continued  in  the 
same  manner  as  in  the  FSCP  46  extension.  Testing  of  the  article  was  discon¬ 
tinued  after  11800  spectrum  hours  when  the  L/H  wing  failed  catastrophically. 
This  failure  was  found  to  have  originated  at  the  pylon  fitting  hole  in  the 
wing  main  torque  box  lower  skin  at  B.L.  132.50. 
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4.3  Compilation  of  Test  Results 

Locations  of  the  five  fatigue  critical  key  areas  and  photographs  of 
typical  damage  detected  in  the  laboratory  are  shown  in  Figures  45  through  53. 
Detailed  lists  of  fatigue  damage  detected  in  each  of  the  five  areas  are  pre¬ 
sented  in  Figures  54  through  59. 

4.4  Electron  Microscopic  Examinations  of  Fatigue  Fracture  Surfaces 

The  markings  on  a  fatigue  fracture  surface  represent  successive  positions 
of  the  crack  front.  Electron  microscopy  provides  the  means  for  correlating 
this  crack  growth  to  the  loading  history.  Crack  growth  curves  generated 
during  microscopic  examinations  thus  define  the  fatigue  characteristics  of 
components  in  terms  of  times  to  crack  initiation,  crack  detection,  and  com¬ 
plete  rupture. 

Fracture  surfaces  from  each  of  the  five  key  areas  were  examined  using  the 
electron  microscope.  Crack  growth  curves  generated  for  the  critical  areas  in 
the  wing  main  torque  box  lower  skin,  in  the  F.S.  303  bulkhead,  and  in  the 
outer  wing  lower  skin,  are  presented  in  Figures  60  through  62.  Attempts  to 
obtain  crack  growth  data  from  fracture  surfaces  from  the  wing  main  torque  box 
upper  skin  and  from  two  lower  longeron  dog  bone  fittings  did  not  prove  to 
be  successful.  These  areas  are  loaded  in  compression,  and  although  cracks 
developed  due  to  residual  tensile  stresses,  subsequent  compressive  cycles 
damaged  the  fracture  surfaces  to  such  an  extent  that  crack  growth  could  not 
be  determined. 
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Figure  45 

32*15062  and  32*15531  Outer  Wing  Lower  Skins 
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Figure  47 

Fuselage  Station  303  Bulkhead 
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Fuselage  Station  333  Bulkhead  Failure  F-4J  Test  Article 
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Figure  61 

Dog  Bone  Fitting  (Lower  Longeron  Spiice  Fitting) 
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Figure  54 

List  of  Failures  Detected  in  the  Wing  Main  Torque  Box  Lower  Skin  at  BL  100  in  Laboratory  Testing 
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Figure  55 

List  of  Failures  Detected  in  the  32-15062  and  the  32-15531  Outer  Wing 
Lower  Skins  in  Laboratory  Testing 
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List  of  Failures  Detected  in  the  FS  303  Bulkhead  in  Laboratory  Testir 
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Figure  58 

List  of  Failures  Detected  in  Winn  Main  Torque  Box  Uooar  Surface 
at  Wing  Root  in  Laboratory  Testing 
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List  of  Failures  Detected  in  the  32-32086  Lower  Longeron  Dog  Bone 
Fitting  in  Laboratory  Testing 
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Figure  60 

Crack  Growth  Curves  for  Key  Area  in  Wing  Main  Torque  Box  Lower  Skin 
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5.  PHASE  IV  -  LAB  AND  SERVICE  FAILURE  CORRELATION 

5.1  General 

The  correlation  of  laboratory  and  service  experience  for  the  fatigue 
critical  key  areas  was  performed  in  Phase  IV  of  this  program.  This  work 
included  compiling  data  on  fatigue  cracks  detected  in  service  aircraft  for 
comparison  with  laboratory  test  results  compiled  in  Phase  III,  evaluating  the 
degree  of  simulation  of  the  service  loading  environment  in  laboratory  testing, 
and  comparing  the  service  experience  to  life  predictions  based  on  all  the 
data  gathered  in  the  first  three  phases.  The  combination  of  fatigue  test 
results  scatter  and  loads  usage  severity  scatter  into  a  joint  scatter  factor 
was  a  major  part  of  this  effort. 

5.2  Compilation  of  Service  Failures 

Lists  showing  the  service  experience  for  each  of  the  five  key  areas  dis¬ 
cussed  in  Phase  III  are  presented  in  Appendix  IV.  These  lists  were  compiled 
based  on  studies  of  MCAIR  in-hou9e  aircraft  inspection  records  and  supple¬ 
mental  information  obtained  from  the  Air  Force  and  the  Navy.  As  may  be  noted, 
the  lists  include  information  on  demonstration  team  aircraft  (Blue  Angels  and 
Thunderbirds)  as  well  as  fleet  airplanes,  and  information  on  aircraft  in 
which  fatigue  cracks  were  and  were  not  detected.  The  lists  of  aircraft  with 
cracks  were  developed  mainly  from  Che  MCAIR  records  of  incidents  of  reported 
i  itigue  cracks.  The  lists  of  aircraft  inspected  in  which  no  cracks  were 
detected  were  developed  mainly  from  the  information  obtained  from  the  ser¬ 
vices.  This  Information,  l.e.,  that  used  in  preparing  the  lists  of  aircraft 
in  which  no  cracks  were  detected,  is  discussed  in  the  following  paragraphs. 

For  the  critical  area  in  Che  wing  main  torque  box  lower  skin,  the  list 
reflecting  crack  free  performance  in  service  operations  was  developed  based 
on  hole  inspection  records  compiled  by  the  Air  Pot  ex  during  performance  of 
F-4  Taper-Lok  retrofit  operations.  Tas>er-Loks  ave  being  installed  in  this 
area  of  the  wing  main  torque  box  lower  skin  as  a  result  of  Che  full 
stale  laboratory  test  program  discussed  in  Section  4.2.  Crack  inspection 
information  sheets  were  obtained  from  the  Air  Force  for  474  airplanes  in 
which  Taper-Lok  fasteners  were  installed  in  the  wing  main  torque  box  lower 
surface  at  B.L.  100  on  a  retrofit  basis.  Ait  examination  of  the  data  sheets 
from  these  aircraft  revealed  462  airplanes  in  which  all  fastener  holes  were 
free  of  crack  indications.  These  are  the  462  aircraft  listed  in  Figure  IV-I 
of  Appendix  IV.  There  were  data  sheets  on  twelve  aircraft  in  which 
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indications  of  some  sort  of  discrepancy  were  noted  for  one  or  more  fastener 
holes.  The  sheets  were  not  clear,  however,  as  to  whether  there  were  fatigue 
cracks  in  the  fastener  holes  or  not.  In  gener?  it  would  appear  that  what 
was  described  was  probably  a  scratch  or  an  indication  of  a  stress  corrosion 
crack  since  the  location  of  the  discrepancy  in  the  hole  did  not  correspond  to 
where  fatigue  cracking  would  be  expected.  Efforts  were  expended  in  an  attempt 
to  substantiate  whether  or  not  the  twelve  aircraft  in  question  contained 
fatigue  cracks  at  retrofit.  Although  no  positive  information  could  be  obtained 
on  the  particular  twelve  aircraft,  conversations  with  cognizant  personnel  at 
the  repair  facility  indicated  that  fatigue  cracks  were  not  detected  on  any 
aircraft  during  the  retrofit  operation.  Since  the  results  of  the  twelve 
inspections  could  not  be  verified,  the  aircraft  were  not  included  in  this 
study  and  data  from  the  aircraft  are  not  presented  herein. 

In  the  case  of  the  critical  area  in  the  lower  surface  of  the  outer  wing, 
the  lists  reflecting  crack  free  performance  were  developed  based  on  an  outer 
wing  inspection  of  Air  Force  airplanes  in  March  1970.  The  results  reported 
to  MCAIR  were  by  bases  and  included  the  number  of  aircraft  checked  and  the 
Bureau  Numbers  of  airplanes  in  which  cracks  were  detected  (the  aircraft  with 
cracks  are  included  in  the  list  of  failures).  MCAIR  records  of  aircraft 
stationed  at  the  various  bases  at  that  time  were  examined  to  determine  the 
airplanes  that  were  inspected  and  found  to  be  crack  free.  These  aircraft  are 
the  ones  listed  in  Figure  IV-4  of  Appendix  XV. 

For  the  critical  area  in  the  F.S.  jO"*  bulkhead,  the  list  of  fleet  air¬ 
craft  indicating  crack  free  performance  in  service  operations  is  composed 
primarily  of  aircraft  inspected  by  the  Navy  (a  small  portion  of  this  list 
and  the  lists  of  failures  and  of  demonstration  team  aircraft  with  no  cracks 
were  developed  from  data  from  a  limited  nuaber  of  bulkhead  inspections  con¬ 
ducted  by  MCAIR),  As  a  result  of  testing,  HCAIR  proposed  certain  modifica¬ 
tions  for  the  F-4  including  a  modification  of  the  F.S.  303  bulkhead 
incorporating  an  eddy  current  examination  prior  to  rework.  This  modification 
of  the  bulkhead  was  approved  by  the  Navy  for  the  rework  of  a  selected  number 
of  F-4B's  and  F-aJ’s.  A  limited  number  of  airplanes  have  been  inducted  for 
repair  and  are  currently  being  reworked.  None  of  these  aircraft  have  been 
found  to  contain  cracks.  These  airplanes  are  the  ones  coeuprising  the  bulk 
of  the  list  in  Figure  IV-7  of  Appendix  IV. 
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5. 3  Correlation  of  Laboratory  and  Service  Loadings 

The  development  of  reliability  procedures  for  estimating  time  to  early 
failure  on  fighter  aircraft  was  the  basic  objective  of  this  program.  The 
material  presented  in  the  following  section  provides  a  quantitative  evaluation 
of  methods  developed  toward  this  end.  In  this  final  phase,  service  life  pre¬ 
dictions  are  compared  to  the  actual  service  failures.  The  life  predictions 
are  based  on  data  gathered  in  the  first  thre*  phases,  i.e.,  laboratory  test 
results,  counter  and  VGH  data,  and  usage  and  fatigue  test  results  scatter 
information.  The  comparison  itself  is  on  the  basis  of  equivalent  laboratory 
test  hours  and  not  in  terras  of  actual  flight  hours  at  failure.  This  type  of 
comparison  attempts  to  account  for  the  differences  between  the  laboratory  and 
service  loading  environments. 

In  laboratory  testing,  all  efforts  are  directed  toward  attaining  the 
best  possible  simulation  of  the  loading  experienced  by  the  aircraft  in  service. 
The  problem  of  cost,  however,  as  it  always  does,  necessitates  that  simplified 
test  approaches  actually  be  utilized.  For  the  F-4  laboratory  test  condition 
chosen,  major  structural  components  (e.g.,  the  lower  wing  skin)  were  loaded 
in  much  the  same  way  that  they  are  loaded  in  service.  The  loading  was  an 
approximation,  however,  and  as  a  result,  other  areas  may  not  have  been  loaded 
in  exactly  the  snae  manner  as  in  service  operations.  The  F-4  outer  wing  is 
subjected  to  buffeting  and  small  amplitude  vibrations  in  service.  This  was 
not  simulated  in  the  laboratory  tooting.  The  lower  longeron  dog  bone  fitting 
and  the  wing  main  torque  box  upper  surface  are  primarily  subjected  to  com¬ 
pressive  loads;  failures  in  these  areas  rod a it  from  the  inducement  of  residual 
tensile  stressor  .  Tire  dog  bone  fitting  and  Upper  surface  failures,  therefore, 
depend  very  much  on  when  high  loads  are  applied. 

For  the  reasons  noted  above,  and  aiso  because  of  changes  in  operational 
usage ,  laboratory  test  hours  and  service  flight  hours  could  not  be  correlated 
on  a  one  to  one  basis  for  any  of  the  key  areas  considered  herein,  Such  a 
flight  hour  correlation  would  not  be  exactly  correct  for  either  the  wing  main 
torque  box  lower  skin  or  the  F.S.  303  bulkhead,  and  would  be  particularly 
inaccurate  for  the  outer  wing  lower  skins,  the  lower  longeron  dog  bone 
fittings,  and  the  wing  main  torque  box  upper  surface  at  the  root.  Fatigue 
analysis,  element  test  results,  flight  loads  measurements,  and  full  scale 
test  results  information  have  thus  been  used  to  correlate  the  laboratory  and 
service  failure  experience.  Fatigue  analysis  procedures  were  formulated  to 
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calculate  fatigue  damage  such  that  En/N  =  1  at  the  time  of  failure.  The  same 
procedures  were  then  employed  to  predict  the  En/N  damage  sustained  at  in¬ 
spection  or  failure  in  service  operations.  These  procedures  utilized  element 
fatigue  test  data,  full  scale  airplane  strain  surveys,  and  flight  loads 
measurements.  The  damage  numbers  were  then  converted  to  equivalent  laboratory 
hour.]  by  multiplying  the  En/N  value  times  the  iboratory  life. 

It  should  be  noted  that  comparisons  were  actually  made  only  for  the 
critical  areas  in  the  wing  main  torque  box  lower  skin,  the  outer  wing  lower 
surface,  and  the  F.S.  303  bulkhead.  The  remaining  two  critical  areas  were 
not  amenable  to  analysis  by  the  methods  of  the  subject  program.  As  noted 
previously,  both  the  upper  wing  skin  and  the  lower  longeron  dog  bone  fitting 
develop  fatigue  cracks  by  the  mechanism  of  an  overload  in  compression  produc¬ 
ing  residual  tension  stresses.  The  service  fatigue  life  is  therefore  very 
much  influenced  by  when  m  overload  occurs.  An  overload  could  be  the  result 
of  a  high  positive  maneuver  or  a  hard  carrier  landing.  The  resulting  life 
scatter  is  expected  to  be  extreme.  The  compilation  of  service  failures  bears 
this  out.  For  example,  the  upper  wing  akin  was  found  cracked  on  one  F-4 
after  only  195  flight  hours;  and  on  another  F-4,  no  cracks  were  detected 
after  700  flight  hours. 

5 . 4  Comparison  of  laboratory  and  Service  Failures 

The  actual  comparison  of  laboratory  and  9ervJ.ee  failures  included  efforts 
in  the  following  areas;  equivalent  laboratory  hours  for  all  service  failures 
were  determined;  scatter  in  fatigue  and  in  usage  seven./  were  combined  and 
the  Joint  probability  of  scatter  in  fatigue  and  usage  severity  was  evaluated; 
service  failures  were  compared  to  predicted  failures  on  the  basis  of  the  ex¬ 
pected  number  of  failures  versus  the  actual  number  of  failures  for  each  key 
area;  and  finally,  actual  times  to  first  failure,  second  failure,  etc.,  were 
compared  to  the  times  at  which  failures  were  predicted.  Details  of  these 
efforts  are  delineated  in  the  following  paragraphs. 

S.4.1  Equivalent  Laboratory  Hours  -  The  equivalent  laboratory  hours 
sustained  at  each  of  the  service  failures  (or  each  service  inspection  in 
which  no  cracks  were  detected)  are  listed  in  Figures  IV-1  through  IV-10  of 
Appendix  IV.  As  indicated  previously,  these  hours  were  determined  by  com¬ 
puting  the  En/N  damage  number  for  a  particular  component  and  then  multiplying 
that  number  by  the  life  demonstrated  by  the  component  in  laboratory  testing. 
The  detailed  methods  are  described  in  the  following  paragraphs. 
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5.3  Correlation  of  Laboratory  and  Service  Loadings 

The  development  of  reliability  procedures  for  estimating  time  to  early 
failure  on  fighter  aircraft  was  the  basic  objective  of  this  program.  The 
material  presented  in  the  following  section  provides  a  quantitative  evaluation 
of  methods  developed  toward  this  end.  In  this  final  phase,  service  life  pre¬ 
dictions  are  compared  to  the  actual  service  failures.  The  life  predictions 
are  based  on  data  gathered  in  the  first  three  phases,  i.e. ,  laboratory  test 
results,  counter  and  VGII  data,  and  usage  and  fatigue  test  results  scatter 
information.  The  comparison  itself  is  on  the  basis  of  equivalent  laboratory 
test  hours  and  not  in  terms  of  actual  flight  hours  at  failure.  This  type  of 
comparison  attempts  to  account  for  the  differences  between  the  laboratory  and 
service  loading  environments. 

In  laboratory  testing,  all  efforts  are  directed  toward  attaining  the 
best  possible  simulation  of  the  loading  experienced  by  the  aircraft  in  service. 
The  problem  of  cost,  howe  'er,  as  it  always  does,  necessitates  that  simplified 
test  approaches  actually  be  utilized.  For  the  F-4  laboratory  test  condition 
chosen,  major  structural  components  (e.g.,  the  lower  wing  skin)  were  loaded 
in  much  the  same  way  that  they  are  loaded  in  service.  The  loading  was  an 
approximation,  however,  and  as  a  result,  other  areas  may  not  have  been  loaded 
in  exactly  the  same  manner  as  in  service  operations.  The  F-4  outer  wing  is 
subjected  to  buffeting  and  small  amplitude  vibrations  in  service.  This  was 
not  simulated  in  the  laboratory  testing.  The  lower  longeron  dog  bone  fitting 
and  the  wing  main  torque  box  upper  surface  are  primarily  subjected  to  com¬ 
pressive  loads;  failures  in  these  areas  result  from  the  inducement  of  residual 
tensile  stresses.  The  dog  bone  fitting  and  upper  surface  failures,  therefore, 
depend  very  much  on  when  high  loads  are  applied. 

For  the  reasons  noted  above,  and  also  because  of  changes  in  operational 
usage,  laboratory  test  houis  and  service  flight  hours  could  not  be  correlated 
on  a  one  to  one  basis  for  any  of  the  key  areas  considered  herein.  Such  a 
flight  hour  correlation  would  not  be  exactly  correct  for  either  the  wing  main 
torque  box  lower  skin  or  the  F.S,  303  bulkhead,  and  would  be  particularly 
inaccurate  for  the  outer  wing  lower  skins,  the  lower  longeron  dog  bone 
fittings,  and  the  wing  main  torque  box  upper  surface  at  the  root.  Fatigue 
analysis,  element  test  results,  flight  loads  measurements,  and  full  scale 
test  results  information  have  thus  been  used  to  correlate  the  laboratory  and 
service  failure  experience.  Fatigue  analysis  procedures  were  formulated  to 
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The  laboratory  fatigue  lives  for  the  key  critical  areas  in  the  wing  main 
torque  box  lcwer  skin,  the  F.S.  303  bulkhead,  and  the  outer  wing  lower  surface 
as  determined  from  the  data  presented  in  Section  4.3  are  shown  in  Figure  63. 
Fatigue  analysis  procedures  were  formulated  which  would  predict  these  labor¬ 
atory  lives.  The  same  procedures  were  then  employed  to  predict  the  damages 
sustained  in  service  operations.  The  actual  computations  of  the  En/N  damage 
numbers  sustained  by  the  servi  ^  aircraft  were  made  using  an  existing  F-4 
fatigue  damage  computation  program.  In  making  the  predictions,  the  program 
utilized  each  airplane's  counting  accelerometer  data,  the  overall  fleet 
average  count  frequency  for  estimating  counts  for  those  periods  when  an  air¬ 
plane's  counter  was  inoperative,  and  the  VGH  data  discussed  in  Section  3.2  for 
correcting  counter  load  factor  data  to  wing  bending  moment  or  to  F.S.  303 
bulkhead  load.  Periods  in  which  the  aircraft  were  engaged  in  combat  or  in 
training  were  tracked  and  the  appropriate  VGH  correction  was  utilized.  The 
resulting  En/N  damage  numbers  were  then  multiplied  by  the  component  labora¬ 
tory  lives  to  establish  the  equivalent  laboratory  hours. 

In  addition  to  the  full  scale  test  results  information  of  Figure  63, 
flight  loads  measurements  and  element  test  results  information  were  used  to 
aid  in  determining  the  equivalent  hours.  Flight  loads  measurements  provided 
the  factors  to  be  applied  to  the  counting  accelerometer  data  to  reflect  the 
differences  between  the  flight  regimes  represented  in  the  laboratory  and  in 
actual  operations.  In  Section  3.2  It  was  noted  that  the  majority  of  maneuvers 
are  pulled  in  a  flight  regime  in  which  wing  bending  moment  per  g  is  about  85% 
of  that  used  in  laboratory  testing.  Flight  loads  measurements  indicated  that 
F.S.  303  bulkhead  loads  follow  the  same  relationship  (bulkhead  loads  were 
found  to  be  approximately  proportional  to  the  wing  root  bending  moment) .  On 
the  other  hand,  flight  loads  measurements  indicated  that  the  outer  wing  is 
loaded  more  severely  in  service  in  relation  to  the  laboratory  loading  than 
are  either  the  bulkhead  or  the  wing  main  torque  box  lower  skin.  This  is 
borne  out  by  the  Figure  64  summary  of  average  flight  hours  to  crack  detection 
tor  the  three  key  areas.  The  summary  shows  outer  wing  lower  skin  fatigue 
cracks  at  considerably  fewer  hours  in  relation  to  the  bulkhead  and  the  wing 
main  torque  box  lower  skin  than  the  laboratory  test  results  would  suggest. 

The  factors  from  flight  loads  measurements  thus  had  to  be  included  to  account 
for  the  differences  between  the  "points  in  the  sky"  represented  in  the  lab¬ 
oratory  and  in  actual  flight  operations.  As  indicated  earlier,  element  test 
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results  were  also  utilized.  Demonstration  team  aircraft  pull  considerably 
more  negative  g  maneuvers  than  do  fleet  airplanes;  and,  within  demonstration 
team  groups,  the  solo  airplanes  are  subjected  to  a  more  severe  negative  load 
spectrum  than  are  the  diamond  airplanes.  Factors  for  the  effects  of  negative 
loads  on  the  damages  sustained  by  the  demonstration  team  airplanes  were  ob¬ 
tained  from  MCAIR  element  tests  utilizing  spectra  derived  from  Blue  Angel 
counting  accelerometer  data.  Different  factors  were  developed  for  solo  usage 
and  for  diamond  usage.  Different  factors  were  also  developed  for  different 
time  periods  to  reflect  variations  in  the  negative  load  factor  usage  which 
have  occurred  over  the  period  in  which  F-4  aircraft  have  been  engaged  in 
demonstration  team  operations.  These  factors  were  then  applied  to  the  cal¬ 
culated  En/N  damage  numbers  to  determine  the  equivalent  laboratory  hours 
exhibited  at  the  failure  or  inspection  of  demonstration  team  airplanes. 

As  indicated  above,  the  equivalent  laboratory  hours  listed  in  Figures 
IV-1  through  IV- 10  are  the  equivalent  hours  computed  as  having  been  accumu¬ 
lated  at  the  time  of  crack  detection  or  at  the  time  of  inspection.  These 
hours  are  the  best  estimates  of  the  equivalent  hours  accumulated  through  that 
time.  Where  counting  accelerometer  data  was  missing  and  damage  had  to  be 
estimated  for  more  than  50%  of  the  flight  hours,  that  fact  is  indicated. 

Also,  in  the  case  of  the  outer  wing  lower  skins,  equivalent  hours  were  com¬ 
puted  only  if  there  was  no  indication  of  switching  of  outer  wings  on  the  air¬ 
plane  in  question.  Wherever  possible,  outer  wing  aerial  numbers  for  the 
aircraft  in  the  lists  of  Figures  IV-4  through  IV-6  were  checked  against  the 
serial  numbers  of  the  outer  wings  installed  on  the  aircraft  when  they  left 
MCAIR.  Where  outer  wings  had  been  switched,  equivalent  hours  were  not 
computed  since  the  counting  accelerometer  data  for  the  airplane  might  not 
bear  any  relation  to  the  usage  experienced  by  the  outer  wings.  The  inspection 
information  of  such  aircraft  have  been  listed  for  reference,  but  these  data 
were  not  used  in  any  further  studies. 

5.4.2  Joint  Probability  of  Scatter  in  Fatigue  and  Usage  Severity  -  The 
combination  of  fatigue  test  results  scatter  and  loads  usage  severity  scatter 
into  a  joint  scatter  factor  was  a  major  part  of  the  effort  involved  in  the 
comparison  of  laboratory  and  service  failures.  As  noted  in  Section  2.5, 
scatter  in  fatigue  data  is  reasonably  approximated  using  the  Weibull  distri¬ 
bution,  and  as  noted  in  Section  3.3.4,  scatter  in  usage  severity  can  be 
described  by  the  negative  binomial  distribution.  Since  they  are  statistically 
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independent,  the  joint  probability  density  function  representative  of  com¬ 
bined  scatter  from  fatigue  and  usage  severity  is  given  by  the  product  of  the 
two  distributions.  To  determine  probability  levels,  this  complex  function 
had  to  be  integrated  which  was  not  feasible  analytically.  A  computer  program 
was  thus  developed  to  do  the  integration  numerically.  The  program  was  de¬ 
signed  to  evaluate  th?  function  for  a  preselected  joint  scatter  factor  and 
then  determine  the  area  under  the  curve.  This  area  represented  the  probabi¬ 
lity  of  the  joint  scatter  factor  for  combined  fatigue  and  usage  severity 
being  less  than  or  equal  to  the  preselected  value.  A  number  of  computer  runs 
were  thus  required  to  completely  establish  the  relationship  between  the 
probability  of  a  fatigue  failure  and  the  joint  scatter  factoi  for  combined 
scatter  in  fatigue  and  usage. 

Curves  showing  the  probability  of  failure  versus  the  joint  scatter  factor 
for  fatigue  and  usage  severity  are  presented  in  Figure  65.  These  curves  were 
derived  using  the  means  and  standard  deviations  listed  for  the  4g  usage 
severity  scatter  in  Figure  26  and  the  Weibull  Shape  Parameter  a  of  5.27 
determined  in  Section  2.2.6.  The  Ag  usage  severity  scatter  means  and 
variances  were  used  because  a  major  portion  of  fatigue  damage  is  caused  by 
maneuvering  at  this  load  factor  level. 

It  should  be  noted  that  a  commonly  used  approach  for  obtaining  a  joint 
scatter  factor  for  fatigue  and  usage  severity  is  to  multiply  the  two  scatter- 
factors  to  obtain  a  total  scatter  factor.  This  approach  results  in  an  overly 
conservative  estimate  of  what  the  true  joint  scatter  factor  should  be.  This 
is  illustrated  by  the  comparison  shown  in  Figure  66. 

5.4.3  Comparison  of  Service  Experience  and  Minimum  Service  Life 
Predictions  -  The  actual  comparison  of  service  and  laboratory  experience  has 
been  made  on  the  basis  of  the  expected  number  of  cracked  aircraft  versus  the 
adtvil  number  and  also  on  the  baaie  of  the  predicted  times  to  crack  detection 
v  irsua  the  actual  times.  The  detailed  calculation  procedures  to  obtain  these 
comparisons  are  presented  irt  the  following  paragraphs. 

The  comparisons  of  the  number  of  cracked  aircraft  versus  the  expected 
number  for  each  of  the  key  fatigue  critical  areas  are  presented  in  Figures 
67  through  71.  In  arriving  at  these  comparisons,  the  expected  number  of 
cracked  aircraft  was  found  by  taking  the  summation  of  the  probability  of 
crack  detection  on  each  airplane  inspected.  As  an  example,  if  one  hundred 
airplanes  would  be  inspected  at  the  point  in  time  when  the  probability  of 
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Figure  66 

Comparison  of  Joint  Probability  Scatter  Factor  to  Product  of  Fatigue  end  Usage  Scatter 
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crack  detection  is  .01,  the  expected  number  of  cracked  aircraft  would  be 
100  x  .01  =  1.  The  curves  of  Figure  65  were  used  for  determining  the  prob¬ 
ability  of  crack  detection  for  each  airplane  inspected.  The  scatter  factor 
to  be  used  for  each  airplane  was  found  by  dividing  the  laboratory  life  (i.e., 
the  number  of  laboratory  test  hours  to  crack  detection)  by  the  equivalent 
laboratory  test  hours  accumulated  by  the  particular  airplane.  Figure  65  was 
then  entered  using  this  scatter  factor  and  the  probability  of  a  crack  being 
detected  was  read.  When  counting  accelerometer  data  was  available  for  more 
than  50%  of  the  flight  hours  accumulated  by  a  given  airplane,  the  curve 
denoted  as  "known  usage"  in  Figure  65  was  utilized.  Where  counter  data  had 
to  be  estimated  for  more  than  50%  of  the  flight  hours  accumulated  by  a  given 
airplane,  the  probability  of  failure  was  determined  usin^  the  curves  denoted 
as  "usage  unknown."  As  would  be  expected,  this  resulted  in  a  greater  prob¬ 
ability  of  failure  for  a  given  scatter  factor  when  the  usage  of  an  airplane 
was  unknown. 

The  predicted  times  to  crack  detection  for  a  fleet  of  aircraft  depends, 
of  course,  on  the  number  of  aircraft  in  the  fleet.  For  purposes  of  correla¬ 
tion  in  t u-d  report,  since  all  airplanes  in  the  fleet  were  not  flying  at 
the  same  rate  of  damage  accumulation,  the  calculation  of  expected  times  to 
crack  detection  required  special  consideration.  For  example,  say  the  labora¬ 
tory  test  life  is  3000  hours  and  that  there  is  a  1000  airplane  fleet.  Further 
assume  that  all  1000  airplanes  have  reached  or  exceeded  500  equivalent  labora¬ 
tory  test  hours,  then  the  expected  number  of  cracked  aircraft  is  1000  times 
the  probability  of  failure  determined  from  Figure  65  at  a  scatter  factor 
equal  to  3000/500  «»  6,  The  calculation  at  500  hours  is  fairly  straightforward 
since  all  1000  airplanes  in  the  fleet  have  reached  or  exceeded  500  hours, 
Uowever,  to  calculate  the  expected  number  of  cracked  aircraft  in  less  than 
or  equal  to  1000  hours  is  somewhat  different  if  it  is  assumed  that  all  the 
aircraft  have  not  reached  1000  hours.  Assume,  for  example,  that  750  of  the 
1000  aircraft  have  reached  or  exceeded  1000  hours.  Then  in  terms  of  those 
750  aircraft,  the  expected  number  that  would  be  found  cracked  is  750  times 
the  probability  of  failure  determined  from  Figure  65  at  a  scatter  factor 
equal  to  3000/1000  »  3.  This  is  not,  however,  the  total  expected  number 
out  of  1000,  The  reason  is  that  the  250  aircraft  that  have  not  reached  1000 
hours  also  contribute  toward  the  probability  of  failure  in  less  than  or  equal 
to  1000  hours.  Their  contribution  would  be  calculated  based  on  the 
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probabilities  of  failure  determined  from  Figure  65  at  the  appropriate  scatter 
factor  for  each  airplane;  and  then  each  of  these  probabilities  would  be  summed 
and  finally  added  to  the  expected  number  for  the  750  aircraft  that  have 
reached  or  exceeded  1000  hours.  This  process  can  be  continued  to  provide  a 
graph  of  expected  number  of  cracked  aircraft  versus  hours.  The  graph  can 
then  be  entered  at  expected  number  equals  one  to  determine  the  predicted  life 
for  the  first  cracked  aircraft,  it  can  be  entered  at  expected  number  equals 
two  to  determine  the  predicted  life  for  the  second  cracked  aircraft,  etc. 

The  predicted  lives,  thusly  determined,  are  50%  probable  type  numbers. 

Using  the  technique  discussed  in  the  preceeding  paragraph,  graphs  of  ex¬ 
pected  number  of  cracked  aircraft  versus  equivalent  laboratory  test  hours 
were  constructed  for  each  of  the  key  fatigue  critical  areas.  These  graphs 
are  presented  in  Figures  72  through  74.  Then  entering  these  graphs  at  expected 
number  equals  one  and  two,  the  predicted  lives  for  the  first  and  second 
cracked  aircraft,  respectively,  were  determined  and  are  shown  in  Figures  75 
through  78.  Also  presented  for  comparison  are  the  actual  service  lives  from 
the  lists  of  Figures  1V-1  through  1V-10  in  Appendix  IV. 

5.4.4  Discussion  of  Lab  and  Service  Experience  Correlation  -  A  review 
of  the  lab  and  service  experience  comparisons  in  Figures  67  through  78  indi¬ 
cate  favorable  correlation  for  the  wing  main  torque  box  lower  skin  but  some¬ 
what  leas  than  favorable  for  the  F.S.  303  bulkhead  and  the  outer  wing  lower 
skins.  Pertinent  considerations  are  discussed  in  the  following  paragraphs. 

Figure  67  indicates  that  six  Navy  Blue  Angel  airplanes  were  inspected 
using  the  eddy  current  technique  in  the  wing  main  torque  box  lower  skin 
fatigue  critical  area.  Two  of  these  were  the  aircraft  deployed  in  the  severe 
solo  operation.  One  had  accumulated  1830  equivalent  laboratory  test  hours 
and  the  other  had  accumulated  2100  hours.  Fatigue  cracks  were  detected  in 
both  aircraft  as  expected  since  cracks  were  detected  in  the  laboratory  at 
1600  hours.  The  remaining  four  airplanes  were  flown  in  the  much  less  severe 
diamond  formation  and  had  accumulated  on  the  order  of  400  hours.  None  of  the 
four  were  found  to  be  cracked.  Figure  75  shows  that  the  predicted  life  for 
the  first  cracked  aircraft  is  1570  hours  as  compared  to  the  actual  value  of 
1830  hours. 

Fip/r*'.'  68  indicates  that  462  F-4  aircraft  were  inspected  using  the  dye 
penetrant  technique  in  the  wing  main  torque  box  lower  skin  fatigue  critical 
area.  None  of  these  were  found  to  be  cracked.  This  is  not  unreasonable 
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500  1000  1500  2000 

Equivalent  Laboratory  Test  Hours 


Figure  72 

Expected  Time  to  First  Cracked  Aircraft  and  Second  Cracked  Aircraft 
for  the  Key  Area  in  the  Wing  Main  Torque  Box  Lower  Skin 
(Demonstration  Team  Airplanes) 


2500 


OP  73-0430 


Figure  73 

Expected  Time  to  First  Cracked  Aircraft,  Second  Cracked  Aircraft,  Etc., 
for  the  Key  Area  in  the  FS  303  Bulkhead 


QP7  3  0439-7 
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Figure  74 

Expected  Time  to  First  Cracked  Aircraft,  Second  Cracked  Aircraft,  Etc.,  for  the  Key  Area 
in  the  Outer  Wing  Aft  Lower  Skin  (32-15062  Block  26  and  Up  Configuration) 

OI*730«3B-8 
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inner  Wing  Lower  Skin  Time  to  Crack  Detection 
(Equivalent  Laboratory  Test  Hours) 

Cracked  Aircraft 

Expacted 

Actual  | 

1570  hrs 

1830  hrs  1 

■ 

>2100  hrs 

2100  hrs  j 

Figure  75 

Comparison  of  Expected  and  Actual  Times  to  First  Cracked  Aircraft  and  Second 
Cracked  Aircraft  for  the  Kay  Area  in  the  Wing  Main  Torque  Box  Lower  Skin 

(Demonstration  Team  Airplanes) 

QP7V043B-0 
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FS  303  Bulkhead 

Time  to  Crack  Detection 
(Equivalent  Laboratory  Test  Hours) 

Cracked  Aircraft 

Expected 

Actual 

First 

1530  hrs 

407  hrs 

Second 

2140  hrs 

410  hrs 

Third 

>  3500  hrs 

672  hrs 

Note;  Unfavorable  correlation  between  predicted  and 
actual  llvet  due  to  fabrication  vartationi.  See 
Section  5,4.4  and  Figure  69. 


Figure  76 

Comparison  of  Expected  and  Actual  Times  to  First  Cracked  Aircraft, 
Second  Crecked  Aircraft,  Etc.,  for  the  Key  Area  in  the  FS  303  Bulkhead 

0*73  043*13 
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Outer  Wing  Aft  Lower  Skin 
Time  to  Crack  Detection 
(Equivalent  Laboratory  Test  Hours) 


Cracked  Aircraft 

Expected 

Actual 

First 

■rfqh 

Second 

l 

Third 

1395  hrs 

Fourth 

2450  hrs 

SEE  1 

Nota:  Unfavorable  correlation  between  predicted  and 
actual  live*  due  to  outer  wing  buffeting.  See 
Section  5.4.4  and  Figure  70. 


Figure  77 

Comparison  of  Expected  and  Actual  Times  to  First  Cracked  Aircraft, 
Second  Cracked  Aircraft,  Etc.,  for  the  Key  Area  in  the  Outer  Wing 
Aft  Lower  Skin  (32-15062  Block  26  and  Up  Configuration) 


Outer  Wing  Forward  Lower  Skin 
Time  to  Crack  Detection 
(Equivalent  Laboratory  Test  Hours) 


Cracked  Aircraft 

Expected 

Actual 

First 

4100  hrs 

666  hrs 

Note;  Unfavorable  correlation  between  predicted  and 
actual  live*  due  to  outer  wing  buffeting.  Saa 
Saction  5.4.4  and  Figure  71. 


Figure  78 

Comparison  of  Expected  and  Actual  Times  to  First  Cracked  Aircraft 
for  the  Key  Area  in  the  Outer  Wing  Forward  Lower  Skin 
(32-15531  Block  26  and  Up  Configuration) 


Q»*»C4W9 


since  cracks  could  not  be  reliably  detected  with  dye  penetrant  in  the  labora¬ 
tory  until  2400  hours,  and  the  majority  of  the  462  airplanes  inspected  in 
service  had  accumulated  less  than  1000  equivalent  laboratory  test  hours.  The 
highest  number  of  hours  on  any  of  the  462  airplanes  was  1500. 

The  F.S.  303  bulkhead  fatigue  cracking  originates  in  fastener  holes  in 
the  bulkhead  flange  as  shown  in  Figure  47.  A  more  detailed  view  of  the 
attachment  of  the  bulkhead  to  the  wing  main  spar  is  shown  in  Figure  79.  Note 
that  the  tension  stresses  in  the  critical  flange  result  from  wing  spar  curva¬ 
ture  tending  to  pull  the  spar  cap  away  from  the  "foot"  of  the  bulkhead  through 
the  two  inboard  fasteners.  The  magnitude  of  tension  stress  depends  signifi¬ 
cantly  on  each  fasteners  torque-up,  the  exact  location  of  the  fasteners,  and 
the  stiffness  of  the  bolt  and  nut  combination.  These  types  of  parameters  can 
vary  from  one  airplane  to  another.  It  is  considered  probable  that  these 
fabrication  variations  are  the  causa  of  the  relatively  poor  correlation 
between  predicted  and  actual  lives  for  the  F.S.  303  bulkhead  shown  in 
Figure  76. 

Within  the  speed-altitude  envelope  where  the  F-4  airplane  executes  the 
majority  of  its  maneuvers,  buffeting  ia  fairly  common  at  high  angles  of 
attack.  For  example,  at  the  airplane’s  design  gross  weight  of  37,500  lbs. 
and  at  Mach  0.7  and  10,000  ft.  altitude,  buffet  onset  ia  about  4  g's.  How¬ 
ever,  the  buffeting  originates  and  primarily  remains  in  the  outer  wing  panel. 
The  resulting  vibratory  loads  superimpose  on  the  basic  maneuver  wing  airloads 
to  increase  the  stresses  in  the  outer  wing.  These  vibratory  loads  also  in¬ 
crease  the  inner  wing  stresses,  but  by  a  ouch  smaller  percentage  because  of 
die  relatively  small  contribution  to  inner  wing  stresses  from  outer  wing 
loads.  The  unfavorable  comparison  between  predicted  and  actual  lives  for 
the  outer  wing  lower  skins  in  Figures  77  and  78  is  considered  to  be  caused 
by  outer  wing  buffeting  which  was  not  simulated  in  the  laboratory  testing. 
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6.  SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 


The  conclusions  which  may  be  formulated  from  this  four  phase  program  on 
the  evaluation  of  structural  reliability  analysis  procedures  as  applied  to  a 
fighter  aircraft  are  as  follows: 

(1)  Scatter  in  spectrum  fatigue  tests  is  considerably  less  than  that  in 
constant  amplitude  fatigue  testa.  This  was  found  to  be  true  for 
both  element  and  full  scale  fatigue  teat  specimens.  For  the  full 
scale  test  article  results,  the  comparison  is  as  follows: 

Spectrum  Constant  Amplitude 


Number  of  Specimens  243 
Number  of  Groups  82 
Average  Standard  Deviation  ,0985 
Average  Weibull  Shape  Parameter  5.44 


491 

143 

.1486 

3.70 


(2)  The  scatter  in  spectrum  testing  exhibited  by  7075  and  2024  aluminum 
in  both  element  and  full  scale  specimens  is  generally  about  the 
same.  Pooling  all  of  the  data  gives  the  following  averages: 

Average  Spectrum  Fatigue  Tear  Data 


Number  of  Specimens  1060 
Number  of  Groups  260 
Average  Standard  Deviation  («)  .0994 
Average  Weibull  Shape  Parameter  (a)  5.27 


(3)  Comparisons  of  the  Weibull  and  the  log-normal  probability  distribu¬ 
tions  to  the  actual  spectrum  fatigue  test  data  indicate  that  the 
Weibull  distribution  (shape  parameter  a  «  5,27)  provides  a  better 
fit  of  the  data  than  the  log-normal  distribution  (standard  deviation 
o  -  .0994) 

(4)  In  addition  to  studies  of  experimental  data,  theoretical  analyses 
were  performed  yielding  the  mathematical  probability  distribution 

for  a  Weibull basad  scatter  factor.  A  laboratory  test  article  can  be 
thought  of  os  one  airplane  selected  at  random  from  the  total 
fleet.  The  scatter  between  the  laboratory  life  and  the  service 
life  of  another  airplane  picked  at  random  from  the  fleet  is 
then  the  ratio  of  two  randomly  selected  variates  free  the  same 
population.  A  scatter  factor  so  defined  as  the  ratio  of  two 
statistically  independent  random  vari  hies  i«s  itself  a  random 
variable  and  its  probability  distribution  can  thus  be  derived  from 
the  distribution  oC  the  parent  population.  The  relationship  between 
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the  reliability  R  (probability  of  no  failure)  and  the  scatter  factor 
S  tor  a  parent  population  described  by  the  Weibull  distribution  is 
_1 

R  a 

3  "  <T=r> 

Using  this  formula  with  a  =»  5.27  gives  a  fatigue  scatter  factor  of 
2.39  for  992  reliability. 

(5)  A  total  of  8200  hours  of  F-4  VfiH  data  were  analyzed.  The  trends 
from  this  data  are  as  follows i 

(a)  The  majority  of  maneuvers  in  both  combat  and  training  operations 
are  executed  in  a  limited  Mach  number /altitude  regime.  In 
both  types  of  operations,  the  majority  of  maneuvers  are  pulled 
at  between  350  and  550  knots  and  at  below  10,000  feet  (the 
combat  average  altitudes  being  slightly  higher  than  these  for 
training  due  to  ground  fire  avoidance) . 

)  The  gross  weights  for  maneuvers  pulled  in  combat  are  higher 
than  thc.se  for  maneuvers  pulled  in  training  operations  (due 
to  the  higher  weapon  payload  required  in  actual  combat  service), 
(e)  In  both  training  and  euohat,  the  average  speed  at  which 

maneuvers  are  pulled  increases  as  the  load  factor  increases 
(higher  airspeeds  required  in  order  to  pull  high  load  factor 
maneuvers) . 

(6)  More  than  2,000,000  flight  hours  of  F~4  counting  accelerometer  data 
ware  studied  to  determine  usage  severity  scatter  trends.  The  data 
indicate  marked  reduction  in  scatter  with  increasing  flight  hours. 
This  trend  reflects  the  fact  that  the  longer  aircraft  are  in  service, 
the  more  likely  they  will  be  subjected  to  a  variety  of  usages 

and  their  repeated  loading  histories  will  "average  out". 

(7)  The  usage  severity  scatter  exhibited  in  the  counting  accelerometer 
data  was  evaluated  for  "goodness  of  fit"  with  the  negative  !  ir.vmi 
distribution.  Comparisons  for  airplanes  having  accumulated  100, 

500,  and  1500  hours  show  definite  correlation  between  the  theoreti¬ 
cal  distribution  and  the  data. 
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(8)  A  study  was  also  conducted  on  the  total  hour  accumulation  on 
individual  F-4  aircraft  to  obtain  information  which  can  be  used 
to  aid  in  establishing  an  aircraft's  design  life.  This  study 
showed  a  trend  similar  to  that  detected  for  usage  severity  scatter, 
i.e.,  flight  hour  usage  also  averages  out  and  aircraft  placed  in 
the  fleet  at  the  same  time  V7ill  tend  to  accumulate  ilar  numbers 
of  hours  over  a  period  of  years.  In  addition  to  this,  however, 
since  F-4  flight  hour  monthly  accumulation  rates  were  shown  to  be 
relatively  invarient,  the  study  also  indicated  that  all  aircraft 

in  a  fleet  can  be  expected  to  accumulate  a  similar  number  of  hours 
after  a  like  number  of  years  in  service.  The  overall  average  for 
the  F-4  is  25  flight  hours  per  month. 

(9)  The  combined  effect  of  fatigue  test  scatter  and  usage  severity 
scatter  was  derived  utilizing  a  joint  scatter  factor  concept.  The 
total  scatter  factor  derived  in  this  manner  is  significantly  smaller 
than  that  determined  by  an  overly  conservative  simple  multiplication 
approach. 

(10)  Probable  minimum  lives  were  computed  for  three  key  fatigue  critical 
areas  on  the  F-4  airplane  based  on  the  reliability  procedures 
presented  in  this  report  and  on  F-4  laboratory  fatigue  and  usage 
data.  The  correlation  with  actual  service  experience  was  excellent 
for  one  of  the  areas,  but  not  for  the  other  two.  tiavy  An  pel 
airplanes  were  inspected  in  the  wing  min  torque  box  lower  skin 
critical  area  using  the  uddy  current  technique.  Two  aircraft  were 
found  to  be  cracked  as  expected.  For  these  aircraft,  the  predicted 
time  to  detection  of  the  first  crack  was  IS 70  hours  as  compared 
to  the  actual  value  of  1830  hours.  Among  462  fleet  aircraft 
inspected  in  the  sacse  area  using  the  dye  penetrant  technique,  there 
were  no  aircraft  found  to  contain  cracks.  This  was  not  unreasonable 
Since  cracks  could  not  be  reliably  detected  in  the  laboratory  by 
this  method  until  2400  hours,  and  the  majority  of  the  fleet  aircraft 
5 tad  accumulated  less  than  1000  hours.  The  correlations  for  he 
critical  areas  in  the  F.S.  303  bulkhead  and  in  the  outer  wing  lower 
skins  were  not  as  favorable.  This  is  attibuted  to  fabrication 
variations  and  to  outer  wing  buffeting. 


The  theoretical  reliability  procedures  outlined  in  Reference  (1)  and 
expanded  in  this  report  for  use  on  fighter  aircraft  provide  satisfactory 
results.  The  formulation  of  systematic  techniques  for  the  incorporation  of 
these  methods  in  designing  for  structural  reliability  in  fighter' aircraft  is 
now  needed.  Two  basic  approaches  could  be  defined;  (1)  method  for  fail  safe 
components,  and  (2)  method  for  non  fail  safe  components.  Method  (1)  would 
determine  scatter  factor  magnitudes  for  fail  safe  components  based  on  main¬ 
tenance  versus  scatter  factor  trade  off  studies  to  give  minimum  total  system 
cost  at  a  given  performance  level.  The  trade  off  studies  would  be  made  dur¬ 
ing  each  airplane's  design  stage.  Method  (2)  would  determine  scatter  factor 
magnitudes  for  non  fail  safe  components  such  that  the  failure  probability 
within  the  design  lifetime  is  extremely  small.  The  exact  magnitude  of  this 
failure  probability  would  be  determined  during  each  airplane's  design  stage. 
It  should  be  noted  that  these  methods  will  yield  different  scatter  factor 
magnitudes  for  different  components  on  the  same  airplane. 

The  method  for  fail  safe  components  would  utilize  the  concept  of  tine 
to  first  failure,  second  failure,  etc.,  to  determine  how  many  aircraft  would 
require  maintenance  action  in  a  given  time  period.  This  type  of  analysis 
would  be  used  in  the  maintenance  versus  scatter  factor  trade  off  studies. 

In  addition,  this  same  type  of  analysis  would  be  used  to  define  inspection 
intervals.  The  method  for  r-n  fail  safe  components  would  be  based  on  re¬ 
quiring  an  extremely  small  probability  of  failure  for  individual  airplanes. 
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APPENDIX  I 


LIST  OF  FATIGUE  TEST  DATA  REFERENCES 

1-  F~15  Fatigue  Design  Analysis  Report ,  MCAIR  MDC  A0928,  1  Jan.  1972. 

2. '  Evaluation  cf  F-4  Blue  Angel  Fatigue  Life,  KCAIR  TM  253.224,  25  June  1970 

(Also  TRJ32-253. 32,300-107). 

3 .  Effect  of  interspersing  negative  load  level  sequence  spectrum  fatigue 
tents ,  MCAIR  Report  604-110  (Also  TR  604-110,10,  -116,  -155,  -180), 

18  August  1965. 

4.  Repeated  cyclic  load  testing  for  ECP  613  elements,  MCAIR  Report  F623, 

30  Dee.  1968  (Algo  T’R  32A-6l6.12.10,  32A-602).  ’ 

5 •  A  comparison  of  the  fatigue  capabilities  of  7075-T 6  a.. A  Titanium  6-4, 

MCAIR  Final  Report  604-124,  21  August  19 66  Ulso  604-124. 10,  -152,  -159, 
-162,  -176). 

6.  A  comparison  of  the  fatigue  capabilities  of  2020-T6,  2024-T81,  7075-To 
aluminum  and  6AL-6VA-23N,  6AL-1FA-1M0,  6AL-4VA  Titanium,  MCAIR  Final 
Report  604-125,  -125,  -157,  16  Jan.  1967. 

7.  Comparison  of  fatigue  life  for  spectra  of  various  slopes,  MCAIR  Report 
604-179,  8  August  1966  (Also  TR  604-179-10,  -194), 

8 .  Effect  on  the  fatigue  life  of  inserting  mandrel  prior  to  installation 
of  Taper-Loka ,  MCAIR  Report  6o4-212,  10  Jan.  1967.  (Also  TR  32-253.32, 
300-197). 

9.  Effect  of  negative  loads  on  F-4  Blue  Angel  Fatigue  Life,  MCAJR  Report 
TM  253.226/25  June  1970  (Also  TR  J32-339-13). 

IQ .  Evaluation  of  Taper-Lok  fasteners  installed  in  titanium  and  aluminum 
combinations ,  MCAIR  Report  701-117,  8  Jan.  1970  (DC-10), 

11 .  Effect  of  fatigue  cracks  on  the  residual  static  strength  of  7075-T651 
aluminum,  MCAIR  Report  604-287  (EMA) ,  1  May  1969. 

12.  Fatigue  and  exfoliation  corrosion  properties  of  7178  and  7075  aluminum 
alloys ,  MCAIR  Report  513-547,  25  Aug.  19 66. 

13.  Fatigue  evaluation  of  coining  fix  for  the  F-4  303  bulkhead,  MCAIR  Report 
TR  J 32- 339. 12,  10  March  1971. 

14.  Effect  of  negative  loads  on  F-4  Th underbird  Fatigue  Life,  MCAIR  "IM253.393, 
2  Dec.  1970  (Also  TR  J32-339.13.10). 

15.  Effect  of  negative  loads  on  spar  element  specimens ,  MCAIR  Report  J32-339* 
13.11  (F-4),  25  Feb.  1971. 

16.  Effect  of  negative  loads  on  F~4  Fatigue  Life,  MCAIR  TM253.333,  7  Aug. 

1970,  (Also  TR  J 32- 339. 17) • 
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17.  Effects  of  countersinking  depth  on  Fatigue  Life,  MCAIR  Final  Report 
J32-3U0.11,  7  July  1970. 

18.  Evaluation  of  crack  growth  characteristics  for  7Q75-T651  aluminum, 

MCAIR  Final  Report  604-339,  19  June  1969. 

19.  Engine  Auxiliary  air  door  cyclic  test  (F-4J),  MCAIR  Final  Report  32A- 
569,  20  June  1967. 

20.  Evaluation  of  the  resistance  to  fatigue  of  the  F-4  wing  main  torque 
box  lower  skin,  MCAIR  Final  Report  604-304,  28  Feb.  1969. 

21.  Investigation  of  various  coatings  as  lubricants  between  stringer  and 
clip  to  reduce  fretting  ( DC-10 ),  MCAIR  TM  253.24,  30  April  1969  (Also 
TR  E0  1A1003) . 

22.  DC-10  stringer- rib  clip  fatigue  test  for  stringer  loading,  MCAIR  Final 
Report  701-110  (DC-10),  13  June  1969. 

23.  Element  fatigue  tests  evaluation  of  short  edge  distance  (F-4),  MCAIR 
Report  TR  J32-339.15.10,  Jan.  1972. 

24.  Fatigue  element  tests  F-4  outer  wing,  lower  skin,  MCAIR  Report  J32-340.10, 
12  Feb.  1971. 

25.  Fatigue  evaluation  of  F-4  outer  wing  fatigue  critical  area,  MCAIR  Report 
J32-340.12,  1  Feb.  1971. 

26.  Fuselage  station  303  bulkhead  element  test  for  5/32"  skin-duct  fastener 
hole,  MCAIR  Final  Report  J32-339.15,  3  Aug.  1971. 

27.  Tapered  .joint  fatigue  life  evaluation  ( DC-10 ),  MCAIR  Report  701-133, 

22  Dec.  1970. 

28.  Stringer  .joint  fatigue  life  evaluation,  MCAIR  Final  Report  701-137, 

29  January  1971. 

29 .  Fracture  surface  examination  of  DC-10  fatigue  test  specimens,  MCAIR 
’Hi  256.285,  18  Feb.  1969  (Also  TR  701-101.10). 

30 .  Fatigue  and  stress  corrosion  tests  of  2024  shot  peened  specimens, 

MCAIR  Final  Report  701-101,  13  Nov.  1969  (Also  701-101.10  and  .11). 

31.  Effect  of  Taper-Loks  on  the  fatigue  life  of  vs  (X)  element  test  speci¬ 
mens  (2024-T351,  7075-TS51  aluminum  alloys) ,  MCAIR  Final  Report  E6610- 
177,  16  March  1967. 

32.  Element  fatigue  test  with  a  reversed-cycle  spectrum  (7075-T651  aluminum 
alloy),  MCAIR  Final  Report  604-219,  24  April  1967. 

33.  Element  test  to  investigate  the  repeated  load  failure  on  main  torque 
box  upper  skin  on  blocks  1,  6,  and  8  test  articles,  MCAIR  Report  32- 
540. 04  (P4H-1),  18  Dec.  1962. 


34.  Fatigue  test-airplane  wing  splice,  MCAIR  Report  TM  TR  4477.6.1*  12,  13, 
April  1956. 

35.  Improvement  in  fatigue  life  of  7075-T651  plate  by  various  means  of 
surface  working,  MCAIR  Final  Report  32-433  (F4H-1) ,  23  Aug.  i960. 

36.  Element  fatigue  test  to  evaluate  short  edge  distance,  MCAIR  Final 
Report  032-722  (F4H-1),  16  Dec.  i960. 

37 •  Fatigue  tests  to  determine  possible  increase  in  fatigue  life  of  alumi¬ 
num  through  special  processing  techniques,  MCAIR  TRM  32-132.1,  -132.2, 

14  Jan.  1953. 

38.  Fatigue  tests  to  determine  possible  increase  in  fatigue  life  of  aluminum 
through  special  processing  techniques,  MCAIR  TRM  32-132,  PTM2,  13  Jan. 
1959. 

39 .  Fatigue  test  of  specimens  to  evaluate  BUWEPS  spectra  and  stresses  typi¬ 
cal  of  F4H-1  airplane  wings ,  MCAIR  Report  32-650.08,  .08.01,  10  August 
1961. 

40.  Evaluation  of  a  fatigue  damage  indicator,  MCAIR  Final  Report  604-298.10, 

10  April  1969. 

41.  Relationship  between  size  of  fatigue  crack  and  fatigue  life  remaining, 
MCAIR  Final  Report  604-266,  7  Feb.  1968. 

42.  Effects  of  local  pre-stressing  on  fatigue  properties  of  7075-T6  aluminum 
plate ,  Briles  Metal  Laboratory  Test  Report  #151,  26  Sept.  1963 • 
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17Q9* 

1.3669 

^ui 

16 

084101 

3 

3101 

.0057 

3116 

7*. 1351 

<JC3 

16 

084101 

2 

2101 

1 2038 

~  24?F 

2.0993 

2C* 

16 

0S4201 

2 

56917 

.1629 

.  75401 

.  1*1732 

<J05 

16 

084801 

2 

165C01 

.0430 

170775 

9:8951 

2Q6 

17 

130U10 

4 

321716 

.5013 

56*575 

.7652 

23’ 

17 

130010 

2 

37000 

.0000 

37000 

VERT  w|OH 

209 

17 

130010 

2 

64722 

.0569 

67731 

7**879 

23  9 

1  T 

130010 

2 

lTWo’™""' 

.0228 

13738' 

16.7063 

J  « 

17 

130010 

2 

109496 

.0028 

110000 

VERY  hjuh 

^  .  » 
w  •  * 

17 

130010 

C 

14*91 

.0212 

1*733 

20.Q937 

J  *  O 

17 

130010 

3 

278*9 

.1307 

3J166 

3.7*69 

d'.'i 

17 

130010 

2 

75299 

.0448 

780*5 

9* *980 

b;* 

17 

130010 

4 

_ _ _i.745B3.__ 

.3037 

2*3*58 

.1*5669 

2:5 

17 

130010 

4 

751208 

.521* 

1232377 

*975* 

2:t 

1’ 

1300 1C 

3 

7030 

.0135 

7132 

56.8011 

-f «  ? 

17 

130010 

3 

35397 

.0647 

3773? 

6.1132 

2  Id 

17 

133010 

3 

79837 

.133* 

8986* 

3.3610 

1’ 

13301C 

3 

138553 

.0*2* 

14*208 

9.6142 

.;?... 

, . j3aoic.„ 

,.  6 

452962. ._ 

782?  __ 

912913 

.6370 

£‘4  * 

1 

»  9 

13301c 

2 

30463 

,0302 

31208 

14.0826 

s!-_ 

i? 

133013 

3 

18757 

.0866 

2Q197 

5,9565 

»  * 

*  ■ 

130.013 

.2 

5*77 

.0560 

572? 

7.6036 

«  4  •■* 

Id 

043001 

2 

66017 

.0360 

67830 

12*1670 

IS 

13- 1.3C 

3 

S32? 

.  tS»6 

962* 

2*6978 

2  it 

i*  . 

_ J3>»3<0 

rj 

. £269g . . 

s'Sn  „ 

2***0 

2.31?? 

«?2? 

iz 

083801 

? 

3672 

.0436 

3803 

9*7620 

*?; 

$8  *  lii  W 

w 

158833 

.0890 

173261 

*.8198 

<.'♦ 

OP 3103 

2 

201*08 

.0615 

£11553 

6*92*1 

232 

‘2  t 

*  « 

060100 

£ 

137625 

,0886 

i*?725 

4*8255 

231 

22 

CsOlO'O 

3 

320723 

<  »S2l 

*92**6 

.896* 

•  *  ;, 

.  22 

..-.  35-200  ^ 

. . 3S63&  . 

. ,1182. . . 

...,,, .JD5037. 

,,->6ce*  , ..  . 

2S.t 

22 

CP’i'vC 

* 

92847 

•  I*?! 

106013 

3**333 

■»’  *  * 

X  v 

26  220  C 

C 

69£i6 

,C3*i 

71132 

12**735 

•  -•  i 

CS  0£CC 

3 

V 

iss^y* 

,  059* 

16ei?S 

4,7600 

c  it- 

£  £ 

28*222 

*4 

c 

33*873 

•  12*1 

369536 

3,4*6* 

;.  .*  * 

*.  € 

2? 2202 

** 

■£ 

i«il37* 

.35U 

;*2as8 

37*5305 

£  }  * 

t  Jf 

,  .;?a<O0i  .. 

......  3. 

-i  7*'3£. 

.10*40 

,  .£*2269 

v  T  '* 

■£  J1 

2*0005 

5 

5**? 

.0719 

3  v*62 

6**130 

• 

ij 

28303; 

3 

97*6 

*1037 

10*33 

4.5338 

i  •  ; 

ii 

CP  2021 

3 

5  5*57 

•  rate 

205SS 

*•5982 

C  *  ' 

is 

Cft-oo; 

Ti 

£’v5; 

♦  C’27 

28798 

6.7409 

T 

US 

5*3021 

3 

522289 

•  031* 

123770 

21.3285 

i«» 

£i 

.  1. 

-  ,0436. ,  . 

.. ..  t??e*e 

S.£?eS 

;•  - 

<,) 

553251 

? 

62822 

•  335* 

86572 

7* 6*52 

2«  5, 

C  1 

'v5322i 

c 

480111 

«C*2& 

0285S 

10*1042 

■>  * 

2J 

35323: 

? 

*$295 

.6161 

90379 

28.1929 

2  3 

r-fc  -*-»r  * 

V»v6V* 

to 

irrsvf 

.0510 

178200 

VE.H*  MjiJA. 

i  * 

£  § 

^  ^  *1 

3 

26*597 

•  CibS 

£77960 

U»989* 

The  UN«i*SEU  P8PT  ESTIMATES  8?  L8U.N93MAL  ANU 
"HBULL  POPULATION  PASAH^TFHS  FgK  L>ATA  GKOUPS. 


U9U.V8HMAL 

L80.N8HMAL 

■m 

It£M 

HE? 

OESCHI^TISN 

SAMPLt 

SUL 

SCALE 

(Myi 

SHAPE 

(SIGMA) 

sca*e 

(tiETfci 

SHAPE 

(ALPHA! 

25q 

a-* 

080101 

3 

89? 

♦  Q922 

98* 

*.3119 

H51 

e* 

080101 

a 

33070 

.1*52 

371*2 

2.9310 

a<» 

3 

....  ,Q**9 

8*9659 

253 

ab 

130101 

5 

19330 

.1150 

213*2 

5.3*** 

25-* 

26 

080101 

.  2 

?53Q 

.032* 

7727 

13.15*2 

255 

ae 

090101 

2 

U?5S 

.0275 

16086 

15.*6?7 

256 

26 

08010.1 

3 

6190 

.0560 

6*53 

11,960? 

25? 

26 

080101 

5 

66-** 

.0*01 

6917 

10.7307 

25S 

36 

090101 

5 

.  1BD8 

.0952 

2063 

7.1379 

259 

26 

090101 

3 

2570 

.0613 

2?i0 

7,7829 

26Q 

26 

Q9Q101 

5 

27*0 

.0897 

2988 

5.1687 

£6l 

26 

caoici 

3 

13709 

.0659 

1*631 

5*998* 

26a 

36 

Q8Q101 

3 

13313 

.029* 

1*310 

13.6892 

263 

26 

0.80101 

3 

10163 

.1783 

11900 

2»52?0 

_56,T 

26  ...  .  083.101  ... 

3 

laa?i 

,jt£U£. 

J2181 

2.0587 

265 

26 

0801C1 

2 

9988 

.1079 

10885 

3»9*51 

266 

26 

CS01Q1 

.  3 

.  .27*3 

.05*3 

395Q 

7.2756 

26? 

26 

090201 

3 

5356 

.2073 

6**7 

3*2382 

268 

26 

090201 

3 

.  6^62 

.  .2321 

9321 

g.8273 

269 

26 

090201 

e 

6600 

.0*10 

6820 

10,3818 

.  .  2?3 

M  . 

,.,.£222511  ... 

,3856£.., 

. .  .15*6.  . 

4*J56 

*n 

28 

020130 

3 

209710 

.06** 

221607 

9*0182 

272 

28 

022120 

3 

3*3*3? 

.  .3356 

*1702? 

1*9630 

2?3 

39 

080131 

5 

.0323 

6SS3 

**60c9 

.2!?* 

39 

.  easui 

. .  .  t  ..... 

.  %M5 

...  -*igpo 

.  5*85 

0.0700 

2?5 

39 

0*0131 

6 

23898 

•  I8U0 

28521 

2*6516 

.. 

.Jl_ 

_ SS0185  __ 

. . . 2 _ 

_ _ 6?££0* _ 

_ *C5&?  _ 

_ £2218?. 

i,Q»27 

277 

30 

0*0100 

6 

1 7755-9 

.1792 

207067 

3.5*12 

2?8 

30  . 

333100 

....  ,5  . 

..—£97**8  .. 

. .11*7 

290*58 

..  .3.9085 

279 

32 

033-lLC 

2 

£3559? 

.1161 

32*333 

3,6682 

283 

30 

0*01.50. 

£ 

. 56127? 

,07*3 

595616 

..a*?3p# 

32 

0*0  It’S 

c 

239333 

.1919 

2**109 

2.2183 

„4»2 

.>£  . . 

. 3.  ... 

_ .23?? 

— - 
31*850  6* S3** 

283 

3: 

0*2 100 

C 

u 

_ __ 

290005 

.088* 

28m 

30 

C  *010.0 

4 

£86*35 

.  #05*1 

303*62 

7*  kin 

285 

3s 

C*3lCQ 

4 

1C5.*0S 

.0382 

409*56 

ic.«3o; 

286 

33 

0*0100 

s 

.  ii5S28  . 

*0?5£ 

12*3*1 

e*im 

287 

0*0200 

3 

15217S 

*  0*3? 

.151359 

9.0791 

288 

_  33_. 

_£»£??* , 

...3  ... 

_U&8*8 

*19*0 

_ 3Pi*« 

2*055*. 

289 

32 

032400 

4 

83555 

.018* 

85069 

23.0*08 

290 

33 

.  020*00 

3 

.17133? 

.56?* 

112*26 

6»203? 

291 

3^ 

\»C*i  + 

3 

13279S 

«m* 

15*939 

2*5*10 

29£ 

33 

4 

36*1?* 

•  18*2 

**1922 

£.3*7? 

293 

30 

020100 

£ 

335812 

.2*2-5 

*276*7 

5,7555 

.  ....... 29*  .32  .. 

_ _£$£l£tL_ 

s  . 

. 1*3325 

.  .  ,208? 

_ _ ;?66i5. 

.  2*112*. . 

295 

3^ 

>*• 

C0S1O0 

3 

102*53 

«C6?1 

109586 

6*063? 

236 

31 

C3S402 

3 

3839 

.,C5*2 

9*19 

*171*7 

297 

31 

030002 

3 

*632 

.5501 

*861 

7,966V 

£36 

31 

■  032002 

£ 

6563 

.0*52 

6526 

*1650* 

299 

31 

030002 

2 

9*65 

*  llv-S 

1033* 

3.8850 

136 


THC  UNBIASED  PBlNTTiTlMATE*  8*  UB^.NBNMAt  ANtf 
WE.L8UUW  POPULATION  PARAMETERS  POP  UATA  (1R6UPS. 


item  REF 

DESCR1WTJ0N 

SAMPLE 

SCALE 

SHAPE 

SCALE 

_ g aatvm: _ 

share 

Hit.. 

t«U) 

ihtaj 

t ALPHA!  . 

.  .  3.0  Q  31- 

030.002 

..  -.  A 

_ 68*8 

.1  C|5i _ 

.  9133 

6*?16{L  . 

301  31 

080002 

3 

5*16 

.0699 

2601 

5.65*8 

. 308  31 

.....  080002 

2 

17380 

.0**9 

16016 

9.67*0 

303  31 

080002 

2 

2388 

.0686 

2*86 

6.1099 

. .  *Q4  32 

.  080002 

. 2  . 

*82  ... 

_ iMM 

669 

_ it  1081.  .... 

30S  32 

080002 

2 

1117 

•  2129 

1326 

2.0000 

306  32 

080002 

£ 

.  xi.483 

. 3510 

—.MIX?  .. 

307  33 

08*221 

2 

8371 

.0000 

5393 

VERY  HUM 

...  .308  3% 

.08*121 

* 

163B 

.1103 

-..  1315 

_ 3 f 1*53 _ 

309  3# 

08*121 

3 

5215 

*  1177 

5629 

3**536 

.....  310-3* 

.-£8*121 

3. 

5185 

-i2?*6  . 

_ _fc71l 

.  .1*467? . . 

311  3* 

08*121 

3 

2356 

*1161 

2618 

3.7476 

312  3* 

.  .  08*121 

2 

_ *£55  . 

...  x  1273 

*710 

...  3.348*  . 

313  33 

0**121 

* 

1676 

.0316 

1727 

14.1638 

31*  35 

_ Q5U£1 _ 

3 

21*8 

.0*25 

2231 

315  33 

08*121 

3 

2056 

.0722 

2192 

6.2363 

316  35  , 

_ 08*121  . 

...  . .  3  ... 

_ _ . 

.  .0798  - 

. 2U5P 

6*1648 

3i  7  33 

0**121 

3 

2200 

.0297 

2258 

IT  ‘"♦13 

.  315  35 

06*121 

.  .  .3  . 

...  ..  ..  *S33 

*0878 

631? 

6*6825 

31S  35 

08*521 

3 

less 

.03*1 

19*1 

14.3343 

,._.3£o.  .25  . 

£•5*121  .. 

3 

*138* 

5089 

.„._£*  3DSt..  _ 

32 1  39 

08*121 

3 

686? 

.0787 

73*S 

8.7817 

322  35 

. ,  08*121 

. .  a 

,, .  11216  ... 

.0363 

am?, 

323  35 

08*121 

3 

6*0* 

.  1220 

93.82 

3.5S67 

.32.*  as 

.  G&*1£1 

.  ,..,3  ,. 

.  ...  .  9133  _ 

.....  *1126 

1020.4 

..  3.S17* 

325  35 

0**121 

3 

3776 

.P2&* 

3662 

£0«25c6 

326  35 

.  08*121  . 

■ 

„ _ 12 1C  5.. 

*0361 . 

._i£3.*6 

..  10 i, 93,77., 

m  n 

08*121 

3 

18**1 

.0785 

13223 

5*2703 

. 326  .  35 

,  06*121 

.....  .3  ... 

.12372 

.,  .*  107*  . 

..  ...  13350 

,5*1018  ...  ..... 

329  35 

08*121 

3 

11098 

.  13*0 

12518 

3»if69 

-  ....  330-  33  . 

.  £3*121 

_ 3.  . 

5556 

,*2260  ... 

.  .  6169 

%$52  . . 

335  35 

cm  2i 

* 

52336 

.0*79 

5513? 

1**2332 

_  _J3£  .  £6  . 

_ Cjcoai . 

£ 

. .  .781  ... 

..*0502  .  ... 

_J3£ 

._.S.4lS?7 . . 

333  36 

030001 

6*6 

.0272 

17.756* 

33*  38 

.  C830U1 

.  .  * 

...  *39 

....  -,*069* 

*87 

10*4500 

335  36 

0800 31 

* 

MO 

.0330 

*96 

13»«0B7 

336  36 

CBS021 

* 

iiUi 

.  -4.17S3 

an 

.  3* £51?  . 

337  36 

CSCOO’ 

* 

637 

.0583 

?nS6t 

338  36 

__C2205i. 

. *  .. 

706 

.  ..*0513.. 

_ JJ*.. 

lO.ixsS?  .  . 

335  36 

C83U01 

69C 

.03*3 

715 

'  11.4013 

.3*3  37. 

.  .  .06*121 

_  *.  ... 

.  _ ^3*33  . 

..  .-  *18%S 

*88.8 

2  s-3223 

3*i  37 

05*121 

2 

.1316 

3*2S6S 

3*2  37 

0?*«£i  . 

.  S.  . 

..._.  *C3« 

. *0733 

*3?4 

3*3  37 

cm  20 

V 

11*04 

.1150 

U76* 

S«  3789 

3**  37 

.  5 

. 38*fi  . 

_ H8?. 

3*5  35 

03*121 

* 

6205. 

.§89* 

676 f 

4*578® 

3*4  38. 

. ..  CB*12l.  . 

% 

.. 

«SSS?  . 

sm 

3*7  36 

05*125 

4 

rasg 

*3667 

•  7708 

?  «2-\26 

3*6  36 

o.s*m 

* 

.  Ul?t  .. 

.O?** 

uw 

£*4fe?6 

3*9  3S 

03*121 

i  £ 

-;*ao 

.0502 

1589 

S.H*¥ 

r  >  <Jr  ..  .  ' 

■  V,  !  ‘if  ;  > 


THE  UNBIASED  P5I*T  ESTIMATES  SF  LOb.NBRMAW  ANV 
WE16UUL  PCPULATI gN  PARAMETERS  FSR  DATA  QR6UPS. 


ITEM  REF  DESCRIPTION  SAMPLE 


361  38 


363  41 

3M._JL§. 

355  43 

356  42 

357  42 


359  4? 


361  42 

M2__4£_ 
363  42 


369  42 


371  42 


080210 

J58jQSiQ 

080210 


osoeio 


0 


385 

45 

Q800CQ 

366 

45 

060000 

388  45 

189  45 

39J2  43 

391  46 

392  45 


396  43 
396  45 
397~~45 
39?  45 

399  45 


080000 
080000 
_o aopsa 
080000 


080000 

_QSl0U01L 

Q80000 

-OAQOOQ. 

030000 


«®N8RMA 


E 


UQ»«MgKMAL  WEJBULU  WEJBUU 


SIQMA) 


42i9730 


5*1510 


2 

75Q4S 

4.3632 

?. 

121407 

•  0606 

_ J27434 

7*3242 

£ 

87988 

i  0662 

92768 

6.4324 

2 

104923 

.0234 

106905 

1B» 1864 

2 

75973 

.Oleg 

7696J 

£6.3359 

.  .1163 


£.0914 


Be  8773 
5.888a 


35.5347 
VERY  HI MH 


00*6473 


4.2407 


6.1147 


'"F  UNCASED  ►*9I~T  ESTIMATES  BR  L8u. NORMAL  ANU 
aEIBULL  MB**LL  AT  I  ®N  PARAMETERS  F0R  UATA  GROUPS. 


L8O.N0KMAL 

L8U.N0NMAL 

AEIBOLL 

WEI  BULL 

ITEM'' 

HER 

DESCRIPTION 

SAMPLt 

SCALE 

SWAPS 

SCALE 

SHARE 

SUE 

(M|J) 

(SIGMA) 

(BETA) 

(ALPHA) 

“CO 

“5 

083000 

E 

394b? 

.2087 

43569 

1  *  8729 

“Cl 

“5 

C80000 

5 

27131 

.3586 

39571 

1.1927 

402 

“b 

08000C 

5 

65D“4 

.2179 

78639 

3* 1874 

“C3 

“5 

C8C0CC 

5 

206215 

.0364 

215754“ 

10* 3721 

“C“  '. 

“5 

080000 

5 

2“24bc 

.  1633 

293476 

2.2655 

405 

“5 

C8000C 

5 

5032 

.1273 

5737 

3.5138 

“06 

45 

C8000C 

5 

3456 

•  2325 

4422 

2*0725 

“07 

45 

08000C 

5 

2437 

•  1341 

2752 

4*4628 

“03 

“5 

C8CO0C 

5 

246106 

.5469 

436271 

*7761 

“03 

4b 

C800CC 

e 

133&0 

.0652 

14074 

6.5315 

“1C 

4b 

csoooc 

5 

2357 

9  2246 

2943 

1*9445 

“11 

4b 

ceoooc 

5 

31444 

.0675 

33716 

6*7466 

“12 

<*b 

ceoooc 

5 

12528 

.  19“4 

15722 

1*9278 

“13 

4b 

080000 

5 

3105 

.  1554 

3736 

2.4120 

“l«k 

4b 

CSOOOC 

5 

79474 

.5264 

143691 

.7340 

“15 

08000C 

5 

7339 

.0670 

7831 

7*6045 

“16 

4b 

CBCOQO 

5 

22114 

.  1696 

26419 

2*7699 

“17 

45 

080000 

5 

349369 

.5196 

611483 

*7738 

“IB 

4b 

C8C00C 

5 

5644 

.1659 

6486 

4*5940 

“19 

4b 

08000C 

5 

141335 

.4800 

251510 

*8037 

“5.1 

43 

08 0000 

5 

30/3 

.2076 

3728 

2.5598 

“21 

4b 

030000 

5 

38344 

.1880 

46781 

2*3798 

“22 

40 

080U0C 

5 

4733 

.5041 

7805 

.8744 

“23 

45 

oaoooc 

6 

37744 

.1005 

41462 

6*8535 

45 

08C00C 

<j 

3730 

.0487 

3924 

8*6911 

“25 

45 

OSOUOO 

6 

12420 

.0733 

13389 

7*2345 

“26 

4/ 

ORQlOC 

5 

15890 

.0681 

17184 

5*7672 

“27 

50 

08010C 

3 

4915 

.0777 

5260 

8*0635 

“2S 

50 

C80-10C 

2 

116*6 

.0153 

11844 

27*0269 

“23 

54 

080100 

2 

10723 

.0429 

11098 

9*9190 

“30 

50 

080100 

2 

13934 

.0286 

14257 

14*8656 

“31 

W  **■ 

V 

oeouco 

2 

264g 

.0116 

2674 

36*7316 

“32 

50 

080000 

3 

37124 

.0569 

38969 

9*4652 

“33 

50 

080000 

2 

42533 

.0454 

44106 

9*3679 

“3“ 

50 

080000 

2 

107“/C 

.1624 

122367 

2*6213 

“35 

50 

080000 

2 

10599 

.  0058 

10648 

73*4721 

“36 

50 

080000 

2 

8831 

•  0763 

9387 

5*5795 

“37 

53 

080100 

2 

17435 

.0528 

18186 

8*0669 

“35 

53 

080100 

2 

51497  . 

.  1707 

59027 

2*4936 

“39 

53 

080100 

C 

93096 

.0855 

99680 

4.9799 

53 

C801C-Q 

?■ 

201988 

.2429 

245282 

1  *  ’.'523 

““i 

53 

ceoioc 

2 

9049 

.0034 

9100 

VE HY  HlUH 

““2 

53 

0*0100 

? 

52411 

.1508 

59126 

2*8228 

““3 

53 

08010C 

2 

93465 

.0164 

94701 

25*9187 

“4“ 

53 

C8010C 

.  _  2 

398618 

.3623 

532499 

_ W75A  . 

“45 

53 

csoioc 

? 

21424 

.1421 

24000 

2*9966 

“46 

53 

080100 

2 

47486 

.0902 

51037 

4*7191 

447 

53 

csoioc 

8. 

195712 

.0938 

21094s 

4.5396 

44  A 

5  3 

08 0 10C 

2 

R4«5 

•  0362 

8734 

11*7698 

449 

53 

0801CC 

2 

102468 

.1221 

112971 

3*4875 

THE  UNSJASEU  ESTIMATES  OF  ueu-NBRUM.  ARL> 


‘Eiauu.  population 

PARAMETERS  F0K 

OAT*  GROUPS. 

Lgu-s’tJKiAL 

LBO.NBKMAL 

WEIBWLL 

WE  I  BULL 

1  T  EM. 

HEf 

DESCRIPTION 

SCALE 

SHAPE 

SCALE 

SHAPE 

...SIZE 

IRU) 

(SIGMA) 

(BETA) 

(ALPHA) 

*»53 

53 

083100 

2 

_1496fc 

«Q41Q 

15465 

10*3818 

451 

53 

080100 

2 

39999 

.2387 

50381 

1*4748 

452 

53 

C80100 

2 

57271 

.2204 

66306 

1*9312 

453 

59 

080100 

2 

127983 

.0551 

13374B 

7*7262 

<*54 

53 

080100 

£ 

19442 

.047  3 

20191 

8*9931 

455 

53 

G8010C 

2 

39949 

o  0307 

40943 

13*8514 

456 

53 

080100 

2 

447317 

.1721 

513291 

2*4735 

467 

53 

oeoioo 

2 

25396 

.3234 

32888 

1*3163 

458 

53 

080100 

2 

34409 

.0446 

35657 

9*5467 

459 

53 

080100 

2 

90464 

•  1937 

105611 

2*1983 

460 

53 

08C100 

2 

16970 

.0362 

17468 

11*7701 

“61 

53 

080100 

2 

47738 

.0642 

50252 

6*6300 

*62 

53 

CSQlOC 

2 

164561 

.0447 

170554 

9*  513B 

463 

64 

080100 

6 

39430 

.0520 

41588 

8*6847 

464 

52 

oaoioc 

4 

3266 

.0502 

3425 

11*4351 

465 

53 

080100 

4 

4377 

.0446 

4555 

12*  009B 

466 

52 

080100 

4 

6083 

.0464 

6365 

9*5666 

46T 

52 

080100 

4 

11650 

.  0689 

12452 

7*3098 

468 

62 

080100 

4 

29850 

.074J 

31997 

6*2628 

*69 

52 

08010C 

4 

5821 

.0367 

6323 

7*3799 

470 

52 

08C10C 

4 

10393 

.0502 

10998 

7*6495 

471 

52 

080100 

4 

13339 

.0673 

14250 

7*1003 

472 

52 

08Q100 

4 

28127 

.0540 

29697 

10*0999 

473 

52 

QS0100 

5 

68800 

.0506 

72571 

9*7482 

474 

52 

08Q100 

4 

10504 

.0720 

11278 

5*7658 

475 

52 

08C10G 

4 

15619 

.1308 

17863 

3*0444 

476 

52 

08010C 

4 

30874 

.0489 

32382 

9*1562 

477 

52 

080100 

5 

71419 

.0727 

76841 

5*8101 

478 

52 

080100 

4 

36313 

.0755 

38909 

6*3518 

479 

52 

C8C100 

4 

71684 

.1588 

82462 

3*7241 

480 

52 

C8010C 

4 

206074 

.  1936 

245468 

3*1839 

481 

52 

OfiOlOC 

4 

135103 

.095u 

153010 

5*  2£36 

482 

52 

08010C 

4 

377988 

.2075 

457997 

2*8959 

483 

53 

080100 

£ 

315434 

.  1428 

353573 

2*9813 

484 

53 

08Q 100 

£ 

12049 

.0025 

12100 

VERY  HJUH 

485 

53 

080100 

£ 

27964 

.3057 

35704 

1*  3927 

486 

53 

080100 

£ 

45166 

•  2513 

55214 

1*6941 

487 

53 

080100 

£ 

215434 

.  1675 

246305 

2*5410 

488 

53 

080100 

£ 

205790 

.0935 

222833 

4*2767 

489 

53 

080100 

£ 

12489 

.0246 

12  nr 

17*3194 

49c 

53 

Q8Q10Q 

. 2  . 

49799 

.1346 

55455 

3*1632 

491 

53 

080100 

3 

£05462 

.6149 

379335 

•  6441 

492 

53 

. 080100 

.3 

. 376181 

.1871 

445024 

2*3363 

493 

53 

C801C0 

£ 

144gg 

.0638 

15176 

6*6764 

494 

53 

08Q10C 

2 

26049 

.0012 

26100 

VERY  HJ  UN 

495 

53 

080100 

£ 

68992 

.0089 

69486 

47*8171 

496 

53 

...  0801QC 

_  3 

. .456196 

.9761 

1 148027 

*4165 

497 

53 

080100 

£ 

236154 

.3351 

311314 

1*2703 

498 

53 

...  0« Cl 00  ... 

.2  . 

35874  . 

.0513 

37376 

8*2985 

499 

53 

080100 

2 

180685 

.0982 

195434 

4*3371 

140 


the  UNBIASED  P0  JNT  EsTlHATES  0F  UBU.NBRMAU  AND 
wgjBWLL  HBPUUATION  PARAMETERS  f&R  DAT*  QH9UPS* 


LftU.NPHWAL  L8U»NBHHAi.  *EI8QLL  HEI9UUL 


ITEM 

«EF 

DESCRIPTION 

SAMPLt 

SCALE 

shape 

SCALE 

SHAPE 

... 

. - 

...SUE. 

IHU) 

(SIQHA) 

(ALPHA) 

_.  ..  S0Q 

...  .082100  . 

_ 2  . 

....  . ...24453  .... 

.0377 

£5201 

..1L»3Q?3 

£01 

53 

080100 

£ 

89860 

.0342 

92347 

12*4641 

£>02 

S3 

080100 

£ 

124474 

.4361 

176663 

*9718 

SQ3 

53 

080100 

£ 

11489 

.0267 

11737 

15.9324 

SQ4 

53 

.  080.1  OQ  ..  _ 

_ 2_ 

32726 

_ *2725  .... 

40689 

1.54,24 

£05 

53 

080100 

£ 

208560 

.0837 

222986 

5*0882 

SQ6 

53 

.  .  ..  08.01. CO  ... 

....  2 

£3366  . 

e  065o 

24624 

6*4909 

£07 

53 

080100 

£ 

56160 

.£339 

66032 

1*7745 

£08 

53 

OBOlOO 

2 

£2360 

.0685 

23619 

6.2129 

£>09 

53 

080100 

£ 

79371 

.0772 

84423 

5*5163 

...  no 

53 

...  .0.8010.0  ... 

_ 2 

27549 

...  ..  1 1 109 

. 30102 

3*6400 

£11 

53 

080100 

£ 

115242 

.4324 

162821 

.9846 

£12 

53 

..  .080  ISO  ... 

....  2 

004.80 

t.02?3 

10736 

.14.5451 

£13 

53 

osoloo 

2 

38884 

p  2017 

4568fc 

2.1107 

£14 

53 

C8010C 

2 

121588 

.2271 

145796 

1.8742 

£15 

4b 

080100 

5 

1414 

•  3835 

2100 

1.0645 

£16 

.  46  . 

. .  caoioo . _ 

_ 5  .. 

..  _ 1436.  . 

.1873 

._  .  175Q 

2*1265 

£17 

45 

080100 

5 

1703 

.2113 

2071 

2*7234 

£18 

46 

.  QSQlQQ 

.  ...5 

.  _ _ 9J& 

*0687 

1004 

6*8909 

£19 

45 

080100 

5 

1890 

.1491 

£187 

3*0814 

£2n 

45 

oanloo 

5 

. 7.05 

.1024 

777 

5.5917 

£21 

45 

080100 

5 

403 

•  0865 

44Q 

4.6907 

£22. 

45 

. oanloo  ... 

. _S_ 

_ _ _2£6  . 

.144Q 

292 

0*8952 

£23 

46 

oaoioo 

5 

204 

.1108 

232 

3.3890 

624 

.45 

. 080100  .... 

_ 5... 

.....  ..  ....  .  439. _ 

.  ..  *1111 

£55 

4*4207 

£25 

46 

080100 

5 

113 

.2886 

145 

2.4184 

£26 

45 

080100 

5 

232 

.1493 

268 

3.0407 

£27 

45 

080100 

5 

8£2 

.0787 

919 

6*7464 

£28 

45 

_ 080100  ..... 

. . 5. 

. .  .  J>79  .. 

....  .1229 

. 771 

. 3*3664 

£29 

45 

08010C 

5 

781 

•  1726 

924 

2.8477 

.  .  £30. 

.45.. 

_ caoioo _ 

_ J5_ . 

966 

.12246 

. 1192  ... 

.1*9680 

631 

45 

0831OC 

5 

£48 

o  0765 

268 

6.7526 

£32 

45 

080100 

K 

662 

.1465 

640 

2*8585 

£33 

45 

080100 

5 

683 

.0884 

738 

7 • 0780 

£34 

45 

-  .080100 

. .6 

....  .  .  ..  I860 

.1178 

£069 

4*2850 

535 

45 

080100 

5 

193 

.2233 

240 

2*1624 

£36 

45 

080100 

5 

162 

.1543 

186 

3*6731 

537 

45 

080100 

5 

306 

.  8066 

370 

2*5212 

£38 

45 

caoioo 

•  5 

492 

.  ...2563  . 

925 

1*4700 

639 

45 

080100 

5 

175 

.0396 

182 

10*6181 

54Q 

45 

08Q105 

5 

. 186 

-1061 

204 

5*0414 

£n 

45 

080100 

5 

508 

.1223 

578 

3*7244 

642 

.45 

. .  .080100  .  .. 

. 5 

_ .612  . 

.1768 

.  63g 

2*1541 

£43 

45 

080100 

5 

745 

.3819 

1170 

•  9891 

544 

45 

oeniac 

5 

1500 

...2591  .  ... 

1956  . 

1*6438 

545 

45 

080100 

5 

74Q 

.1139 

825 

4.2573 

546 

45 

. . -08C1QC  .  .. 

_ b 

_ 6.6.0 

ji1073 

720 

4*5395 

547 

45 

080100 

5 

1211 

.2565 

1641 

2*2329 

548 

45 

C8Q1QQ 

5 

..  1330 

.1377 

1496 

4*5099 

649 

45 

C8010C 

c 

511 

.0188 

521 

23*5434 

■ML  UNBIASED  ESTIMATES  gF'UN-MWWAL  AND 

*£ I3U.LL  >*ePULAT10N  PARAMETERS  FRH  UATA  QKBUfS.  _ 


L9j-NBKMAL 

L8H.NBKMAL 

aeibull 

WE  I  BOLL 

ITEM 

DESCRMTI9M 

SAMPLE 

SUE 

SCALE 
(MU  I 

SHAPE 

(SIGMA) 

SCALE 

(BETA) 

SHAPE 

(ALPHA) 

55q 

45 

.  08.0  WO . 

....  .2  . 

1343 

•  QQ65 

0 

77.6248 

SSI 

43 

C8010C 

5 

4722896 

.2267 

5812934 

2.3128 

S52 

52 

4 

1279314 

.0321 

1321623 

12.5096 

SS3 

32 

080100 

*i 

3769038 

.2678 

4916528 

1.8734 

354 

45 

_  .c&aioc 

5 

200298 

.2403 

253753 

J..9132 

355 

46 

caoioc 

5 

179198 

.7641 

448034 

.5021 

356 

45 

.  .  080100 

5 

958331 

.3312 

1352169 

1*2539 

357 

45 

080100 

5 

20634 

.3303 

29178 

1.3490 

358 

45 

5 

66682 

.4338 

101338 

1.4008 

359 

43 

08Q10C 

5 

381419 

.2523 

492335 

36C 

43 

08Q10C 

5 

3436 

.2835 

4622 

1*4165 

361 

45 

08010C 

5 

53987 

.1734 

63719 

3.6098 

362 

45 

.080 WO-  ... 

5 

42715 

.1868 

53012 

2*0374 

363 

45 

080100 

5 

722476 

.4884 

1177914 

.8771 

364 

43 

5 

2l 

365 

32 

080100 

4 

417206 

.0933 

454976 

5*0708 

366 

43 

.  ....  080 WO. _ 

...  _J5... 

_ 7.63.496  .. 

.6123 

1414067 

.  8048 

367 

58 

080120 

3 

282 

.0399 

293 

13*3692 

568 

58 

Q9Q12C 

£ 

430 

.0479 

472 

9*1528 

369 

58 

080120 

5 

667 

.0505 

702 

10*2380 

379 

38 

7 

935 

.0966 

1023 

5.B447 

*>71 

38 

080120 

5 

1534 

.  0618 

9*4851 

372 

58 

080.120  ..  .. 

8 

_  2973.. 

.0519 

3156 

7*4867 

373 

58 

080120 

5 

4393 

.1080 

4837 

6*0625 

374 

38 

080120 

.. .  . . 5... 

_ 9.413 

9753 

15*3070 

375 

38 

08012C 

5 

20048 

.0332 

20821 

11*1036 

376 

5? 

<■%  •■Suloni#™!! 

5 

31B29 

.0949 

34563 

6*2431 

577 

38 

080120 

6 

71673 

.2179 

93144 

1*7542 

378 

58 

...  080120 

_ 5... 

_  1852.37  . 

.4979 

312633 

.9116 

379 

38 

083120 

4 

3674221 

.7465 

7263820 

>6542 

380 

58 

Q8Q12Q. 

4 

_ 468.39.434 

.  .2970 

61584375 

1*5986 

381 

37 

080110  . 

2 

37 

,0324 

38 

13.1576 

382 

57 

080110 

2 

.  ..  W2 

0 

30.7252 

583 

57 

08CH0 

3 

802 

SRl 

858 

7*0335 

384 

57 

.  0801.10.  .... 

. _e . 

60l>o 

6000 

Vt«r  HJUH 

385 

57 

080110 

3 

26637 

30302 

3.4650 

386 

37 

. .080.110  ... 

. .3.  . 

.  . . . 122541 

.1002 

133572 

4.9131 

387 

37 

080110 

5 

579561 

.2299 

739292 

1*7466 

388 

57 

3 

3.411687 

4*2138  , 

389 

57 

080110 

3 

95 

,0660 

wi 

6*1079 

390 

37 

080UO 

..  2 

405 

.0499 

422 

8*5301 

391 

57 

080110 

4 

33408 

.1731 

39987 

2*2867 

592 

37 

080H0  ... 

5 

. 111981 

.6030 

220879 

•  6264 

593 

57 

C80HC 

3 

1726/50 

.1726 

2004414 

3*5037 

394 

57 

3 

.2338 . 

1*8324 

395 

57 

080110 

£ 

334 

.0495 

348 

8*6073 

396 

57 

oeciiQ 

....  3 

10811 

.0799 

11621 

5*3992 

397 

57 

080110 

5 

194614 

.3722 

288433 

1*0833 

398 

57 

080110 

.  .  4 

...  63316 

.1273 

70895 

4*4191 

399 

57 

OftOllC 

3 

2152700 

.9479 

6863032 

•  6145 

142 


the  UN«TASEO  P0JnT  £  is  T  lMATES  8F  LOP.mOPMAL  ANO 
"E18ULL  POPULATI0’’-  PASA«tTEKB  F6«  OAT*  GPBUPS. 


L90.MBKMAL 

L0U.M8HHAL 

nEIBUU, 

ME10JLL 

lit" 

WEF 

3£SC°  !KT  18s- 

SA^FOt 

SCALE 

SHAPE 

scale 

SHAPE 

SUE 

<hu) 

tSIQMA) 

ttt£TA> 

(ALPHA) 

*03 

57 

0*0110 

3 

19g58375 

.3639 

26562038 

1.1741 

6CJ 

60 

08C1HC 

2 

16 

.0560 

17 

7*6036 

602 

fcC 

08012C 

2 

123 

.0891 

13B 

4.7800 

6Q3 

60 

0801 1C 

5 

25568 

.1367 

29268 

3.5451 

4)04 

60 

OflOUC 

4 

65529 

.1698 

76786 

2.6114 

605 

60 

C80110 

4 

97650 

.7101 

182138 

.7758 

606 

60 

C8011C 

8 

392408 

.1990 

491683 

2.1971 

607 

60 

CPQ12C 

2 

6 

.0473 

6 

8.9932 

608 

60 

C8012C 

3 

8 

.0783 

8 

,  5 ‘242.8 

609 

60 

080120 

2 

89 

•  0892 

51 

0*6459 

610 

60 

080120 

3 

153 

.0086 

129 

8.1373 

611 

6  0 

C80120 

3 

295 

.0541 

309 

8.7028 

612 

60 

CR'OlSO 

4 

1336 

.  0555 

1*06 

10*0199 

613 

60 

crqu: 

3 

4915 

.0777 

5260 

8.0633 

_ 61  it. 

6Q  ... 

CfiOUC 

_ £_ 

7QW, _ _ 

_ O.00-U0._„. 

.  .2P0Q 

..„Vi.HY_,  HJJH 

615 

60 

C8C1U 

2 

11696 

.0158 

118*4 

27.0269 

616 

6: 

080110 

3 

680218 

.3323 

931531 

1*2225 

617 

60 

080120 

2 

6 

.3374 

5 

. *2619 

618 

6  0 

080120 

3 

-22 

.0*9? 

443 

7.9500 

619 

6  0 

or  mo 

2 

7169600 

.3378 

8659499 

1.7904 

,  •  6<U„ 

„«w„ 

...  - - mO  0 12* 

_ 2... 

.  _ *  .  „ 

. .0685  ... 

_ 4 

6*2126 

691 

60 

ohoiao 

? 

15 

•  03B4 

16 

11.0756 

622 

6  0 

080120 

2 

24 

.0000 

2* 

VtW¥  H l  ON 

683 

to 

OH3120 

r 

96 

•  092R 

106 

4*5861 

6g  4 

60 

080120 

«■» 

c 

3*6 

*0319 

335 

13.3503 

625 

60 

080110 

2 

2793212 

.  3225 

3614*77 

1*3202 

„...__6'26 

,„6.*«, 

.wii 31 20 

g 

.0362 

-.^1*7699 

687 

6  0 

080120 

?  " 

1 1 

.0267 

0 

1&.932U 

62  o 

6  „ 

3»0120 

2 

13 

.0528 

0 

16.7068 

6gss 

60 

OR  3120 

2 

702 

.0*5* 

728 

9*3672 

63" 

60 

080 110 

3 

99*6 

.0*36 

1036* 

J0*0»36 

631 

A  3 

O80UC 

2 

13599 

.0338 

106*8 

73.4721 

6.32 

.  6y..-. 

- . -3.. 

,._...  128666 
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Q 
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2 
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.4  755 

6 

47 
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6 
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t 
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* 

7113 
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46 
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4  1 

023101 

6 

7  7*2 

.0661 

' '"a28Q 

7.113* 

*.  <3  * 

43 

JcOlOl 

t 

78U9 

•  0*53 

el2s 

12*51*6 
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“8 

0201C1 

6 
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7733 

6*7623 
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47 

023101 

6 

1*16) 
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1310* 
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021101 

< 
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6*6622 
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M 
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t 
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»H 

022101 

6 

146*5 

.  0*9 1 

19738 

7.83*2 

6»7 

4  n 

0221CS 

0 

2*16? 

.0.499 

23*68 

9*05*8 
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44 

022101 

c 
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.C4U8 
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44 
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6 

iU2fe 
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THE  UN*JAS£0  k'9I"»T  ESTIMATES  OP  U8W»SCWMAL  AND 
^ElbULL  ►,8P,UI(lATJON  PARAMETERS  PQK  UATA  OHtJU^S* 


I  Tt'i  MfcF 


^  _ 

603  '*<1 

65n  4h 

655  *8 

656  <*a 

657  <»3 

658  OH 

659  «*8 
66*  4J 

661  **4 

662  48 

663  45 


’  3E5V.R:^Tl3>i  SA-out 
.  .  SUE 


aeoioi 

C20101 


020101 


0?CUU1 

C201C1 

ceoioi 

c?oiot 

0201(51 

QgQlQl 

020101 

020101 

020101 

C20101 

020101 


L80»\lSW“,A 


LQb.MBR^AL 


4E1UUI.U  RE I 6ULU 


?»5?98 


25B21 

30369 

20786 

29036 

42625 

.0256 
.0365 
.0413 
.  0576 
.0319 

26701 

31587 

21820 

30765 

44134 

14.7574 

10*8894 

;0«89B5 

7.9822 

12*2961 

"i>050f} 

•  0224 

41442 

'  19.9740 

6917  i, 

.0599 

7J547 

7.4791 

41565 

.0248 

42630 

18*3195 

14699 

.  0436 

15245 

13*4921 

46516 
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49944 

6.5359 

43945 
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46U69 

0*1757 

1 

•  0*34 

31763 

13*8216 

Sr.  J"4*  V 
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3 1 19 1 
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TM£  UNdtASEO  eStJmateS  ftP  l8p.k»RHAL  AND 

-he  i  bull  **®p*ju  at  ion  . parameters  fow  uat*  qkbups* 


LBb»\DKMAL 

L8U.N8XMAL 

4EI0ULL 

HEIBULL 

ITEM  KEF 

0EStRl"Tl8N 

SAMPLt 

SCALE 

SHAPE 

SCALE 

SHAPE 

. SUi. 

(My) 

ISjqha) 

tttET*) 

(ALPHA) 

7qq  57 

Q20122 

3 

2*05 

•  0?3q 

256? 

6*0519 

701  67 

020122 

3 

9266 

»Q569 

9752 

7.7010 

702  67 

020122 

3 

7.3359 

703  67 

020182 

3 

6135 

7.12*7 

70*  67 

020122 

* 

*88* 

.0902 

5303 

6.9390 

70S  67 

020122 

3 

Sl*% 

.0709 

8693 

6*3707 

706  67 

Q2Q1E2 

.  6 

..  . . 1628 

.  .11*3 

18*5 

3.4796 

707  57 

020122 

3 

1237 

.0702 

1313 

9.5281 

70B  58 

6 

15373 

_ .094? 

169*9 

5*4888 

>09  58 

C80121 

6 

21991 

.0716 

TTiTTl 

"  7.3*44 . 

/IQ  58 

C3Q1H1 

.  7 

21129 

.0954 

2301* 

6.0705 

711  63 

080121 

6 

23214 

.0502 

2*6*1 

9*0132 

712  S3 

Q8Q121 

6 

.  .  ...  1336$ 
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13828 

15*7739 

713  58 

080121 

7 

2C351 

.0848 

21852 

9.5653 

71V  58 

..  .  080121  . 

6 

5399 

.0376 

5595 

16.0513 

715  58 

080121 

6 

5613 

.0549 

5959 

6i7irP 

?16  56 

080121 

8 

. 13000 

.0*72 

3,3666 

8.9689 

717  58 

080121 

6 

572? 

.  0560 

6069 

9.4368 

718  59 

0BQ121 

6 

2*0f> 

.07*3 

2981 

8.3738 

719  59 

060121 

6 

3&1* 

.0630 

278? 

6*9712 

„  «a  .t* 

rj 

1 

1 

! 

......  6  . 

...  1903 

_ _ ii956 _ 

18*1 

2*1669 

7&|  59 

080122 

6 
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.  1976 
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—  )n  m*m*  m — f  •  — 

2.96  '2 

nz  59 

.080122 

.  .  6 

-  7695 
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12.9297 

723  59 

080122 

6 
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7.5*0* 

n*  59 

warn  . 

.  ,6 
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.1056 

S?50 

5.**?3 

?f5  59 
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6 
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.0983 
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*.7790 
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6 
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_ ,_iC3*9  , 

50* 

^>0*6264 

727  *7 

080102 

ft' 

195 
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030102 

-,7 

. . .  192 

•  0679 
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729  *7 

080102 

0 

1«6 
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6 
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6 
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9*0857 
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060102 

6 
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6 

80 
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8* 
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6 
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.0588 
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080105 

6 
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1.1*5997 

739  *7 

C8Q102 

6 

2G18 

.0383  ' 
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080102 

6 
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t 
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5 
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.0617 
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6*2835 
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6 
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6? 
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70 
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6 

61 
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63 
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6 

66 
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72 
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6 
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T*E  UNBIASED  PgjNT 

'ESTIMATES  '8F  LBU-NSRPAL  AW 

"El BULL 

88POLATION 

PARAMETERS  F0R 

DAT*  CiNBURS* 

Lftb.MBKMAL 

LSU. NORMAL 

WEI8ULL 

WEIBULL 

ITEM 

«EF 

DESCRIPTION  SAMPLE. 

SCALE 

SHAPE 

SCALE 

SHAPE 

- 

-  ..  - 

_ SUE 

<*U> 

(SJQMA) 

_  <“ET\> 

( ALPHA  1 

_  ...750 

49 

080101 

6 

13632 

.0480 

14192 

10*9576 
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49 

080101 

7 
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.0873 
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4.4681 

_  7,5  £ 

49 

080101 

6 
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.0417 

4488 

11*8348 
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49 

080101 

6 
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8.9674 
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49 
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8 

34845 

.0372 
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755 

49 
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6 
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7.1431 
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49 
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6 
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.0773 
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6*3479 
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49 
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6 
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42533 

3.3233 
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49 
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6 
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45 
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5 
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6 
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767 

45 

08^002 

6 

4269 

.0766 

4568 

7*9628 

768 

45 

08QU02 

.......  1.0  . 

819J 

.  Q6U5 

8789 

6.4660 

769 

45 

080U02 

5 

16515 

.1114 

18322 

5*0487 

_  '1 3L. 

45 

080002 

5 

36515 

.0733 

39541 

5.4376 

771 

45 

080002 

6 

880 

.0845 

~STi 

4.520? 

772 

4b 

0PQOO2 

6 

...  ...51731 

.1039 

57741 
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4b 
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7 
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45 
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5 
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the  UNtMASEO  PB1NT  ESTIMATES  BF  LBB,N8RMAL  ANg 
*E1SULL  HBHg^ATION  PARAMETEHS  F0H  DATA  CiRBUPS* 


LSUoMBHMAL 

L05.NBKMAL 

weibou. 

WEI  BOLL 

ITEM 

ref 

OESCRIHTleN 

SAMPLt 

SCALE 

shame 

SCALE 

SHAPE 

SUE 

(MCI) 

(SJGMA) 

(BETA) 

(ALPHA) 

»oo 

63 
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2 
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63 
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2 
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14102 
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3 
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55 
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4 
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55 
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4 
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5*6485 
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55 

084140 

3 
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.0798 

11564 

4*9604 
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55 

08414Q 

4 

14019 

.  0538 

14752 

8*2858 

813 

55 

084140 

7 

124/4 

.0870 

13596 

5*5644 

814 

55 

084140 

4 

5438 

5915 

5*3202 

81b 

55 

04414C 

2 

994 

.0616 

1045 

6*9083 

816 

63 

084140 

2 

1212334 

.1982 

1450*17 

2*1482 

«17 

63 

C8414C 

2 

1507968 

.1312 

1674664 

3*2457 

aia 

63 

C84l‘»Q 

2 

32*03 

0473 

33653 

6*9931 

819 

63 

084140 

2 

141904 

.4795  . 

208191 

» 8873 

.  .  82063 

.  084140 

2 

.  3.9599 

•  1232 

43699 

3*4541 

821 

63 

U84140 

3 

1291677 

.2706 

1641169 

2*4103 

822 

4b 

C841C1 

3 

32689 

•  Q4U5 

33848 

14.7203 

823 

45 

084101 

2 

71386 

•  Q524 

74439 

8*1263 

824 

45 

084101 

a 

52881 

.1808 

61032 

2*3552 

825 

66 

084111 

3 

14317 

.0709 

15351 

5*5826 

. 826 

66 

_ OMJLU 

3 . 

15626 

.0314 

_ 16.1 .15 

_ 12*6343  . . 

827 

66 

084111 

2 

17087 

.0194 

17354 

£1*9365 

828 

66 

084111 

2 

15482 

.1301 

172*8 

3*151“ 

«£9 

66 

084  111 

3 

102bc 

.1207 

11362 

4*3069 

830 

66 

084111 

3 

14071 

*  11*3 

155*3 

4*H46 

831 

66 

084111 

3 

4083 

*09*4 

8779 

6*7064 

.  83g 

.66, 

u... 

3 

6622 

.  ..0816  _ 

7.10.4 

_6«Q140  .. 

833 

66 

C84 1 1 1 

3 

12083 

.  0524 

12**4 

9*1142 

834 

66 

C84i i; 

3 

20911 

.0*93 

2181Q 

10*4648 

835 

65 

024132 

6 

105615 

.0727 

114189 

5*8878 

836 

61 

0201U 

10 

95 

.0979 

105 

4*5053 

8 '3? 

67 

020111 

1  c 

111 

.1022 

128 

3.7763 

833 

6/ 

...JtfQl.U  . 

- 1-0  „. 

98 

.08*0 

107 

5.8752 

839 

6/ 

0201U 

5 

88 

•  0*8* 

95 

10*6947 

84  0 

67 

020111 

1C 

86 

.1020 

96 

4*4184 

1 

6  7 

C  2  0 1 1 1 

K 

124 

.0894 

134 

7*0640 

842 

67 

Ct'jiu 

10 

94 

.1001 

105 

4*6070 

84  3 

67 

020  m 

5 

84 

.0284 

86 

20*9827 

b  i*  •* 

6  7 

£2Cll,L 

. . 1C  ... 

122 

....  1Q6I? 

_ i31 

...  .7*3312 

67 

020  111 

13 

89 

.1223 

101 

3*9104 

846 

67 

020 1 1 1 

1C 

79 

.1033 

91 

3*5667 

847 

67 

CCclll 

5 

83 

.0520 

88 

7*3013 

?  W  * 

6? 

020111 

U 

107 

.1379 

121 

3*7752 

bw-a 

67 

C20U1 

5 

107 

.0887 

117 

4*7768 

THE  unbiased  hqjnt  ESTIMATES  0F  LBu-NBRMAL  anu 
<EiaULU  >J8HUi.AT10N  PARAMETERS  F0R  DATA  GRBUPS. 


LPpoNBHMAL 

LBU-NBKMAL 

REI8ULL 

4EI0ULL 

I  ffcM 

REF 

UE SCRIPT  1 

SAMP^t 

SCALE 

SHAPE 

SCALE 

SHAPE 

SUE 

(Hy) 

(SltiHA) 

(BETA) 

(ALPHA) 

B50 

67 

QPUlil 

10 

*2 

•  06*8 

98 

8.4619 

051 

67 

020111 

5 

07 

.0829 

9S 

5*3466 

052 

67 

Q20111 

1C 

108 

.0991 

142 

4.649* 

053 

67 

C20U1 

10 

104 

.1160 

118 

3.7*34 

ab% 

67 

04Qll 1 

1C 

119 

.0853 

130 

5*9465 

055 

67 

G20111 

1C 

164 

.0804 

180 

S.C-B? 

056 

7  5 

Cfi4l%l 

3 

£S76a 

•  un 

25455 

3.k 737 

057 

70 

C84141 

3 

13804 

.0?36 

14014 

5.4513 

058 

70 

084141 

2 

33008 

.0279 

3465* 

15*2388 

059 

70 

084141 

2 

310/4 

.  2093 

-3673£ 

£.03*3 

060 

7C 

084141 

3 

£4764 

.  1*07 

27994 

3*3333 

861 

70 

084141 

3 

18900 

.0899 

20657 

*>3960 

06£ 

71 

084131 

5 

739139 

.1921 

896545 

£.*167 

863 

71 

084131 

5 

233012 

•  0983 

£54909 

5*0117 

064  71 

p64l31 

3 

212906 

»  Q373 

_JX99J0. 

JL£»3*2S 

065 

71 

084131 

3 

251651 

.0521 

263099 

~  10*7002 

066 

76 

064121 

£ 

.  1?02. 

.02*2 

1706 

17.6113 

867 
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APPENDIX  III 


DERIVATION  OF  PROBABILITY  DISTRIBUTION  FOR  THE 
RATIO  OF  TWO  WEIBULL  VARIATES 


This  appendix  discusses  the  derivation  of  the  probability  distribution 
and  density  function  for  the  ratio  of  the  sample  characteristic  life  of  n 
laboratory  test  articles  to  the  life  of  a  randomly  selected  aircraft  in 
service.  The  scatter  factor  is  written  as 


S  «  0/y  (1) 

where  &  is  the  point  estimator  of  the  Weibull  characteristic  life  (scale 
parameter)  for  a  sample  of  n  failure  times  given  by 

r  n 


0  = 


±  y ...» 

n  L-d  1 

i=l 


(2) 


The  x^'s  in  Equation  (2)  are  the  lives  of  the  individual  laboratory  speci¬ 
mens  and  a  is  the  known  Weibull  shape  parameter.  In  Equation  (l),  y  is  the 
life  of  a  randomly  selected  aircraft  in  service  that  is  a  Weibull  variate. 

To  find  the  probability  distribution  of  the  scatter  factor  defined  in 
Equation  (l),  the  probability  density  functions  of  both  the  sample  charac¬ 
teristic  life  0  as  well  as  the  Weibull  variate  y  must  be  known.  The  prob¬ 
ability  density  function  of  the  Weibull  variate  is  given  by 


(3) 


for  y  >  0  where  a  and  0  are  the  Weibull  shape  and  scale  parameters,  respec¬ 
tively.  The  probability  density  function  of  0  can  be  determined  from  the 
fact  proven  in  Reference  (l)  that  the  statistic 

2nW  =  2n  (M 

has  a  chi-squared  distribution  with  2n  degrees  of  freedom  where  n  is  the 
number  of  failure  observations  in  the  test  sample.  This  gives  that  the 
probability  density  function  for  the  statistic  0  is  given  by 


f(0)  » 


TTiT:<8 


«r-  b®” 


for  0  >  0. 


(5) 


The  probability  density  function  of  the  ratio  S  =  0/y  is 
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=  /  f2 


y  =  |)  dg  for  S  >  0 


where  f(g,  y)  =  f(g)f(y)  since  B  and  y  are  statistically  independent,  and 
f(y)  and  f(g)  are  given  by  Equations  (3)  and  (5).  Putting  f(3,  y)  into 
Equation  (6)  and  collecting  terms  gives 


-fe«r 
[-■#]  -- 
A  =  n  +  (s) 


where 


■/© 


'-r\a(n+l)-l 


If  u  is  substituted  for  Al  —  I  then  the  integral  reduces  to 


ar\s; 


00 

— f  une~U  du. 

i  Y\n+1  J 


The  remaining  integral  is  in  the  form  of  the  Gamma  function  which  gives 


J*  une~udu  =  nr(n) 


This  leads  to  the  final  expression  for  the  probability  density  function  of 


the  scatter  factor 


n+l0an-l 

f(s)  =  -  ,  (12 

(nSa+l)n+1 

Integration  of  Equation  (12}  with  respect  to  S  from  0  to  S  gives  the  cumu¬ 
lative  probability  distribution  of  S  which  is 

n_an 

F(S)  =  — (13 
(1+nS  )n 

If  reliability  H  is  defined  as  in  Reference  (l)  as  the  probability  of  no 
failure  (probability  that  the  ratio  of  life  in  the  laboratory  to  lift  in 
service  is  less  than  S),  the  relationship  between  R  and  the  scatter  factor 
S  is  n  ctn 

»  •  d* 

(l*nS°)n 
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APPENDIX  IV 

LISTS  OF  SERVICE  INSPECTION  RESULTS 
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Bu.  No. 

Flight  Hour*  *t# 
Tim*  of  Retrofit 

F-4C 

637411 

328 

3085.7 

944 

F-4C 

637412 

332 

2079.0 

1340 

F-4C 

637413 

335 

2126.3 

1509 

r-uc 

637414 

339 

2158.3 

722 

F-4C 

637415 

342 

2256.9 

759 

F-4C 

637417 

349 

2253.6 

1226 

F-40 

637418 

352 

2965.9 

1081 

F-4C 

637419 

355 

2947.3 

634 

F-4C 

637420 

357 

2202.6 

8$0 

F-40 

637422 

361 

1572.6 

821 

F-40 

637423 

364 

I668.3 

as 

F-4C 

637426 

371 

1747.4 

783 

F-40 

637428 

376 

1650.0 

1132 

F-40 

637430 

381 

2959.7 

728 

F-40 

637431 

383 

1635.5 

680 

F-40 

637432 

385 

2008.3 

972 

F-40 

637433 

387 

1809.7 

985 

F-4C 

637434 

389 

2039.3 

I486 

F-4C 

637436 

393 

1716.1 

942 

F-4C 

637437 

395 

2008.0 

913 

F-4C 

637439 

399 

2609.3 

755 

F-4C 

637440 

400. 

2C17.9 

896 

F-40 

637442 

403 

2045.9 

1120 

F-40 

637446 

413 

1639.6 

790 

F-4C 

637447 

415 

1159.8 

339 

F-4C 

637448 

41? 

2309.5 

1952 

F-4C 

637449 

419 

1821.3 

1164 

F-4G 

637450 

420 

1864.7 

90S 

F-4C 

637452 

424 

2177.1 

940 

F-40 

637453 

426 

2205.9 

F-40 

637454 

427 

2105.2 

US4 

F-40 

637455 

429 

2038.3 

704 

F-40 

6J7457 

434 

2059.8 

863 

F-40 

—TOT11M 

437 

1894.2 

949 

F-40 

637440 

439 

1992.9 

794 

F-40 

637464 

442 

1602.3 

nS? 

F-40 

637463 

443 

2197.8 

888 
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Figure  BZ-1 

Fleet  Aircraft  With  no  Cracks  Indicated  in  tlia  Key  Area  in 
the  Wing  Main  Torque  Box  Lower  Skin  at  Time  of  Inspection 
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1?3 

warn 

1 

Flight  Hours  at* 
Tlae  of  Retrofit 

637465 

447 

1951.1 

914 

F -4C 

63746a 

455 

1872.9 

1373 

F-4C 

637470 

459 

'476.0 

46; 

F-i»C 

637471 

161 

1966.1 

1023 

[Mi 

637473 

46) 

2047.3 

1065 

F-4C 

637474 

465 

1803.7 

857 

F-4C 

637475 

466 

2092.8 

1)76 

F-4C 

637476 

468 

1760,2 

1278 

F-4C 

637477 

469 

1753.9 

564 

F-4C 

637478 

471 

206). 2 

761 

F-4C 

637479 

472 

1673.7 

925 

F-4C 

637482 

476 

2016.9 

896 

F-4C 

637434 

479 

1784.4 

1660 

F-4C. 

637487 

485 

18)6.) 

1054 

F-4C 

1 

489 

2291.5 

1460 

F-40 

637491 

491 

2201.6 

1841 

F-4C 

637492 

492 

1677.1 

509 

F-4G 

637495 

498 

187) .6 

910 

F-r.C 

637497 

501 

2CTJ2.4 

817 

F-40 

617W 

505 

20)9.7 

1166 

F-fcS 

637501 

50o 

2217.4 

1199 

r-W. . 

637505 

513 

2011.9 

1606 

ESI 

617608 

518 

20)1.0 

764 

F-4G 

637510 

522 

2146.7 

8)0. 

F-4G 

6375U 

525 

212),0 

1)52 

r-4  r 

637515 

529 

20*6.7 

93) 

617416 

5)0 

17)1.1 

629 

6)7520 

5)7 

1525.6 

681 

isa 

637521 

541 

2)42.1 

19)5 

r-4C 

6)7510 

55) 

18)1.1 

867 

r-to 

6)75)2 

556 

2197.0 

1375 

F-40 

6)75)4 

559 

1995.7 

1350 

F-40 

6)75)6 

542 

16)9.4 

428 

F-'G 

■EBB 

544 

16)6.) 

559 

■BQQH 

563 

1632.0 

497 

r-44 

6)7541 

570 

1950.7 

642 

F-4C 

6)7542 

572 

24))>1 

1240 
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Figure  DM  (Continued) 

Fleet  Aircraft  With  no  Cracks  Indicated  in  the  Key  Area  in 
the  Wing  Main  Torque  Cox  Lower  Skin  at  Time  of  Inspection 
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m 

■SHI 

Plight.  Hours  »t  * 
Tise  of  Retrofit 

HBP®EE3i 

P-4C 

637543 

573 

1&12.3 

792 

P-4C 

637545 

576 

1492.1 

603 

P-4C 

637550 

594 

1539.6 

302 

P-4C 

637552 

589 

2553-8 

1512 

P-4C 

637553 

589 

1485.4 

40 

P-4C 

637555 

592 

1716.9 

827 

P-4C 

637556 

59) 

2135.2 

1515 

P-4C 

637562 

603 

1314.9 

948 

P-4C 

637565 

607 

2178.0 

1447 

P-4C 

637566 

609 

1940.2 

871 

P-4C 

637569 

612 

1699.C 

759 

P-4C 

637569 

613 

1261. S 

617 

63757C 

615 

1436.4 

596 

637572 

618 

1473.0 

451 

637574 

621 

1642.3 

899 

637575 

622 

1689.6 

651 

637579 

627 

1954.5 

1170 

637591 

632 

2128*9 

591 

6375s; 

6)6 

1791.7 

603 

F-4C 

63753} 

6)9 

1770.9 

VC96 

P-4Q 

6)9599 

644 

170) .7 

9«a 

637592 

649 

1783.6 

901 

637594 

*5) 

1977.9 

D27 

637595 

654 

1967.1 

992 

F-4C 

637601 

673 

16)0.6 

990 

r~4c 

637602 

67$ 

13)3.9 

6)9 

FHO 

637*07 

604 

217J.7 

125# 

P-4<2 

637609 

696 

2147.0 

991 

p-4C 

637611 

699 

1)73.9 

731 

r-so 

63761? 

699 

SOfe.5 

W# 

P-iC 

63761» 

TCil 

a«*.9 

ItlO 

r-i-3 

6)7622 

701 

M9*.t 

UOO 

F-40 

6)7625 

709 

2007.4 

996 

F-4C 

6J7CA 

710 

1599.9 

715 

F-iC 

6)7625 

713 

1642.) 

#07 

F-6C 

6)76i6 

714 

14*2.9 

67# 

F^G 

6)?*W 

717 

1914.3 

6*3 
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Figure  1V-1  (Continued} 

Fleet  Aircraft  With  no  Cracks  Indicated  in  tire  Key  Area  in 
the  Wing  Mein  Torque  Box  Lower  Skin  at  Time  of  Inspection 
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mm 

UMl 

Plifht  Hours  ft 
of  Retrofit* 

617629 

718 

1967.6 

63763-'' 

720 

19C3.3 

617631 

732 

2290.0 

617612 

723 

1666.0 

617633 

724 

1579. C 

63761  * 

729 

1823.8 

617637 

731 

1990.8 

637618 

734 

1530.1 

61764* 

74$ 

1734.0 

6j76i>6 

748 

1376.1 

637647 

710 

1865.0 

617610 

755 

1898.5 

617664 

762 

931  .«> 

617651 

76J 

1944.1 

617617 

767 

1811.2 

617661 

774 

1932. 8 

617462 

775 

1480.6 

617666 

7*2 

I865.4 

517671 

791 

1917.0 

617672 

791 

20644 

61767} 

794 

1*92.4 

617679 

1697.* 

61?4<H 

20?J  .7 

517W» 

*17 

1*17.7 

41**» 

*20 

16*4.2 

61*491 

*  «2<.“h 

“  u  ,r  1  ,n 

U97.9 

617696 

** 

1»M 

. 

W?S2 

*44 

16VJ.* 

4177171 

<46 

1671.7 

617K4 

to* 

itoM 

61770$ 

to) 

1VH.9 

tms? 

*62 

1427,9 

»j»i 

to* 

1174.9 

IKOI 
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VJ7.7 

6i<*19 

*76 

1644.0 

<46461 

*7* 

1197.4 

MO664 

*fc 

327.9 

Figure  15Z-1  (Continued} 

Fl**t  Awctalr  With  no  Credit  Indicated  in  the  Key  Area  bn 
the  Wmg  Main  Torque  Bon  Lower  Skin  at  Time  o!  Inspection 
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Figure  iiM  {Continued) 

Fleet  Aircraft  With  no  Cracks  Indicated  <n  the  Key  Area  in 
the  Wing  Main  Torque  Box  Lower  Shin  at  Time  of  Inspection 
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tys? 

Bu«  Ko« 

K6C  Cua  No. 

Flight  Hour*  ot# 
71b*  of  Retrofit 

Equivaltnt 
Ub  Hour* 

p-ac 

6A07SA 

1092 

1627.5 

738 

SAC 792 

1105 

2369.7 

1540 

r-4C 

&AOT93 

110? 

2052.3 

1409 

P-AC 

64C796 

1112 

1049.3 

363 

P-AC 

640802 

1122 

1747.2 

1790 

?-AC 

6a0«6 

1128 

1999.0 

1483 

6AOSU 

1137 

1683.1 

&4 

ran 

640*13 

1141 

1897.9 

1288 

csa 

64CS17 

1147 

700.9 

55 

d 

4408*2 

1155 

2050.4 

1748 

m 

640825 

1162 

aso.7 

1693 

P-AC 

640623 

1167 

883.3 

445 

f-ac 

640629 

1169 

2014.0 

1452 

p-ac 

640931 

1173 

2083.7 

1707 

p-ac 

640036 

1183 

1952.1 

996 

p-ac 

640639 

1187 

2370.0 

1881 

f-ac 

640841 

1191 

2166.6 

1964 

F-AC 

640844 

1197 

990.4 

453 

r-Ac 

640647 

1204 

1475.2 

1361 

r-Ac 

640854 

1216 

1169.4 

691 

P-AC 

640655 

1220 

1110.4 

810 

7-AC 

6AOS56 

1221 

1300.0 

786 

P-AC 

640357 

1224 

1217.4 

705 

P-AC 

6A0958 

1225 

1255.3 

909 

640859 

1229 

1344.5 

1C86 

P-AC 

640961 

1233 

1279.5 

610 

P-4C 

640662 

1235 

1219.7 

489 

P-AC 

64O863 

1238 

14OT.C 

805 

P-AC 

640864 

1240 

1497.1 

695 

P-AC 

640669 

1250 

623.3 

132 

P-AC 

640672 

1257 

1162.8 

589 

F-AC 

640875 

1265 

772.8 

85 

P-AC 

6408?v 

1270 

13,32.8 

808 

P-AC 

640879 

1272 

1297.6 

501 

P-AC 

640879 

1276 

3158.1 

668 

P-AC 

640690 

1277 

1166.1 

966 

P-AC 

640691 

1231 

1168.1 

S33 
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Figure  CM  (Continued) 

Fleet  Aircraft  With  no  Cracks  Indicated  in  the  Key  Area  in 
the  Wing  Main  Torque  Box  Lower  Skin  at  Time  of  Inspection 
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wm 

Bu.  No. 

MAC  Cum  No. 

Plight  hours  at  * 
Tiaso  of  Retrofit 

Blulvslsnt 
Ub  Hours 

F-oC 

640883 

1286 

1155.3 

681 

F-4C 

640884 

1239 

1438.8 

977 

P-iC 

640387 

1296 

128C.2 

1030 

F-4C 

6403% 

1304 

1483.8 

741 

F-/.C 

640391 

1307 

1554.7 

750 

640392 

1308 

1266.8 

888 

?-LC 

640693 

1311 

1C45.5 

643 

P-4C 

640396 

1320 

1282.0 

725 

F-iiC 

640399 

1328 

1209.2 

682 

' 

FHjC 

6409CC 

1331 

1574.7 

59C 

?-uc 

6409 ca 

1337 

1326.1 

666 

p-4c 

640903 

1339 

1554.9 

912 

F-4C 

640305 

1346 

1429.7 

529 

FW.C 

640906 

1349 

1607.0 

ion 

F-J*C 

640907 

1353 

1213.0 

583 

F-I»C 

640911 

1365 

1328.0 

492 

F-4C 

640912 

1368 

1529-6 

1384 

7-LC 

640913 

1372 

1001.0 

770 

F-1»C 

640914 

1376 

1118.9 

662 

F-4C 

640915 

1378 

1854.0 

910 

P-4C 

640917 

1385 

1799.0 
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Figure  EZ-1  (Continued) 
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the  Wing  Main  Torque  Box  Lower  Skin  at  Time  of  Inspection 


181 


Figure  12-1  (Continued) 

Fleet  Aircraft  With  no  Cracks  Indicated  in  the  Key  Area  in 
the  Wing  Main  Torque  Box  Lower  Skin  at  Time  of  Inspection 
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Figure  ISC-1  (Continued) 

Fleet  Aircraft  With  no  Cracks  Indicated  in  the  Key  Area  in 
the  Wing  Main  Torque  Box  Lower  Skin  at  Time  of  Inspection 
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Figure  DZ-1  (Continued) 

Fteet  Aircraft  With  no  Cracks  Indicated  in  the  Key  Area  in 
the  Wing  Main  Torque  Box  Lower  Skin  at  Time  of  Inpsection 
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Figure  12*1  (Continued) 

Fleet  Aircraft  With  no  Cracks  Indicated  in  the  Key  Area  in 
the  Wing  Main  Torque  Box  Lower  Skin  at  Tima  of  Inspection 
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Figaro  IV-4  (Continued) 

Ffaet  Aircraft  with  no  Cracks  indicated  in  the  Key 
Area  in  the  Outer  Wing  Lower  Surface  at  Time  of  inspection 
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Figure  DZ-4  (Continued) 

Fleet  Aircraft  with  no  Cracks  Indicated  in  the  Key 
Area  in  the  Outer  Wing  Lower  Surface  at  Time  of  Inspection 
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Figure  Iff-4  (Continued) 

Fleet  Aircraft  with  no  Crack*  Indicated  in  the  Key 
Area  in  die  Outer  VIMng  Lower  Surface  at  Time  of  I  nspection 
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Figure  IE-4  (Continued) 

Fleet  Aircraft  with  no  Cracks  Indicated  in  the  Key 
Area  in  the  Outer  Wing  Lower  Surface  at  Time  of  Inspection 
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Figure  E£~5 

FJeet  AL  jaft  in  Which  Cracks  Were  Detected  in  tha  Key 
Area  in  the  Outer  Wing  Lower  Surface  at  Tima  of  Inspection 
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Figure  TSL5  (Continued) 

Fleet  Aircraft  in  Which  Cracks  Were  Detected  in  the  Key 
tee  tn  the  Outer  Wing  Lower  Surface  at  Time  of  Inspection 
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Figure  iV-7 

Fleet  Aircraft  With  No  Cracks  Indicated  in  tire  Key  Area 
in  the  F.S.  303  Bulkhead  at  Time  of  Inspection 
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Figure  12-7  (Continued) 

Fleet  Aircraft  with  No  Cracks  Indicated  in  the  Key  Area 
in  the  F.S,  303  Bulkhead  at  Time  of  Inspection 
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Figure  12-11 

Cracks  Detected  in  the  Wing  Main  Torque  Box  Upper  Skin  in  Service  Operations 
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Cracks  Detected  in  the  Wing  Main  Torque  Box  Upper  Skin  in  Service  Operations 
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Figure  IE-12 

Cracks  Detected  in  the  Lower  Longeron  Dog  Bond  Fitting  in  Service  Operations 
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where  S  is  the  scatter  factor  and  R  is  reliability  or  the  probability  of  no 
failure.  VGH  and  load  factor  counting  accelerometer  data  from  the  F-4  fighter 
airplane  were  utilized  to  correlate  that  airplane’s  laboratory  and  service 
fatigue  experience.  Probable  minimum  service  lives  considering  the  F-4  fleet 
size  and  individual  airplane  usage  were  computed  based  on  the  Weibull  based 
scatter  factor  and  order  statistics.  The  combined  effect  of  fatigue  test 
scatter  and  usage  severity  acattay  was  derived  utilizing  a  joint  scatter 
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